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When probed at very high energies or small Bjorken a;bj, QCD degrees of freedom manifest themselves 
' as a medium of dense gluon matter called the Color Glass Condensate. Its key property is the presence 

, of a density induced correlation length or inverse saturation scale Rs — 1/Qs- Energy dependence 

' of observables in this regime is calculable through evolution equations, the JIMWLK equations, and 

OA , characterized by scaling behavior in terms of Qa. These evolution equations share strong parallels 

f—j ' with specific counterparts in jet physics. Experimental relevance ranges from lepton proton and 

, lepton nucleus collisions to heavy ion collisions and cross correlates physics at virtually all modern 

collider experiments. 
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1 Introduction 

1.1 QCD in collider experiments at high energies 

With the advent of modern coUiders, from the Tcvatron and HERA to RHIC and planned 
faciUties Uke LHC and EIC, the high energy asymptotics of QCD has gained new prominence 
and importance. Needless to say that these facilities have been built with very different research 
goals in mind. On the particle physics side one finds most prominently the search for the 
Higgs particle as the last missing ingredient of the standard model and hope for clear evidence of 
physics beyond it. Here high energies are needed mainly to exceed particle production thresholds 
to facilitate the generation of telltale signals with sufficiently large cross section. On the nuclear 
or heavy ion physics side it is the search for the quark gluon plasma and its properties that 
requires high energies to create the energy densities and temperatures needed to cross the QCD 
deconfinement phase transition and create a medium of deconfined quarks and gluons whose 
exact nature is still hotly under debate. 

In all these experiments one is talking about enormous center of mass energies per participant 
(i.e. nucleon): RHIC operates at center of mass energies of 200GeV/nucleon, LHC in its heavy ion 
mode is aiming at 2TeV/nucleon while the pure particle physics experiments will be conducted 
via proton proton collisions at 7TeV. This implies enormous boost factors between 200 and 
several thousands that will strongly affect the QCD aspects of these experiments. In such an 
environment one knows that soft gluon emission is enhanced by logarithms that arise from dE/E 
factors in the phase space measure. 
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In QED the analogous process of photon emission can be treated by resummation, which 
gives rise to Sudakov form factors. (These provide the resolution of the well known infrared 
"problem" in QED.) For the non-Abelian case the situation is more involved. The (Abelian) 
Sudakov type of exponentiation or eikonalization is still present, but captures only a part of 
the effects. In QCD the problem of multiple soft emission is intrinsically nonlinear - giuons 
carry color charge and thus in contradistinction to photons act as sources for further emission. 
This provides a growth mechanism for cross sections absent in QED. Complementarily, after 
a sufficient number of soft emission steps, one has to start to consider recombination effects: 
further emission into a region that is already populated by a large number of other color charges 
will be modified by recombination and absorption, a mechanism that will invariably slow down 
further growth in gluon numbers and lead to saturation. Both these effects are irrelevant in QED 
(suppressed by an additional factor of a) and only relevant in the non-Abelian theory. There 
they are usually associated with the successes and failures of the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL-) equation [1-5], which is the prototype of an evolution equation meant to resum soft 
gluon emissions at high energies. This equation carefully takes into account the multiple emission 
aspect and in this sense addresses a specifically non-Abelian issue that is absent in QED, but 
excludes recombination effects that arise at large gluon densities and thus is no longer valid where 
these become relevant. 

The BFKL equation, if taken literally and used beyond this limit of applicability leads to 
an untamed, characteristic power-like growth of cross sections with cm energy, for instance in 
deep inelastic scattering (DIS) as conducted in electron proton collisions at HERA. This in itself 
would lead to a violation of the unitarity bound on cross sections.^ 

On the other side, at vanishing (or just very small) momentum transfer the BFKL equation is 
plagued by a diffusion into the infrared: giuons are emitted at smaller and smaller momenta and 
will eventually reach nonperturbative scales where the assumption of small coupling underlying 
the derivation of the equation is no longer valid. 

Both of these main problems of the BFKL equation can be cured by taking into account 
nonlinear effects in the gluon generation mechanism at high energies, at least for central colli- 
sions, where we expect gluon densities to become large. This is the domain of a highly dense 
gluonic medium and the topic of this review. The main theoretical tool is an associated evolu- 
tion (or Wilson type renormalization group) equation that describes the creation and change of 
the medium with increasing energy, the Jalilian-Marian-|-Iancu-|-McLerran-|-Weigert-|-Leonidov- 
-l-Kovner equation [7-15] [the order of names was chosen by Al Mueller to give rise to the acronym 
JIMWLK, pronounced "gym walk"]. The quite distinct combination of features, the generation 
of large gluon densities which can arise only in non-Abelian theories, the time scale differences 
between soft and hard modes that enter the dynamics and the onset of recombination and sat- 
uration effects have led McLerran and Venugopalan to coin the term "Color Glass Condensate" 
(CGC) which subsequently has been generally adopted as a convenient label for this situation in 
particular with regards to the HERA (electron proton), RHIC (heavy ion) EIC (lepton nucleus) 
and LHC (both in proton proton and heavy ion mode) experiments. 

It should be kept in mind, however, that the central underlying concept, the nonlinear effects 
in soft gluon emission feature prominently also in other areas. An example that provides oppor- 
tunities for fruitful comparison can be found in the study of so called non-local jet observables. 
This is a class of observables in which soft emission from secondary particles -relatively " hard" 
jet constituents other than the original leading partons- plays an important role. Accordingly 
jet evolution can be described by evolution equations that have a, at first sight surprising, strong 
structural and conceptual analogy with the JIMWLK equation. For so called global observables, 

^Unitarity in QCD requires that cross sections asymptotically rise not faster with invariant energy s than 
(Ins)^. This is usually called Froissart's theorem [6]. 
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where the origin of secondary gluon radiation can be ascribed to the original leading partons. the 
equations reduce to liner evolution equations that to a large extent arc dominated by Sudakov 
type physics - to a degree that one can even derive evolution equations without any overt refer- 
ence to medium effects. In the generic case, however, the medium is the crucial ingredient. This 
situation will be used throughout this review to highlight the generic nature of concepts and the 
mechanisms appearing in the generation of a dense gluonic medium at high energy, although 
phenomenology and applications in both cases are quite different. 

Returning to the CGC, one has to highlight its generic feature, the creation of a density 
induced correlation length Rs = l/Qs- The associated momentum scale Qs is called the satura- 
tion scale. They simply characterize the onset of the saturation effects induced by the medium. 
Beyond 1/Qs color charges start to screen each other and since the underlying gluon densities 
grow with energy, Qs will, by necessity, follow this trend. Qs is the characteristic scale of the 
medium and responsible for all the "good" features that go beyond the physics of the BFKL 
equation: Below Qs the growth of the gluon density and thus of certain cross sections is checked 
on the one hand, and on the other hand the infrared problems are cured: a shrinking correlation 
length Rs prevents the appearance of non-perturbative modes near and below the QCD scale as 
long as the dense medium is present, in particular for central collisions. Generically, one would 
expect Qs, as the characteristic scale of the medium, also to be the relevant scale for the coupling 
constant in this situation. This holds true, although it has turned out that the issue is a bit 
more subtle than initially expected. 

As the list of experiments given at the outset already indicates, the situation is quite generic. 
Given the high energy involved in all these experiments, it is quite obvious that soft gluon 
emission and thus JIMWLK like evolution should affect all of them, be it a relatively simple 
experiment such as deep inelastic lepton proton (henceforth referred to by ep or 7*p, as done at 
HERA) or lepton nucleus scattering {eA or 7* A as planned with the EIC experiment) or proton 
proton (pp), proton nucleus (pA) or nucleus nucleus {AA') experiments. [A refers to the atomic 
number of the nucleus involved.] In all of them one expects CGC physics and thus the presence 
of an energy dependent saturation scale Qs to affect particle production rates and cross sections. 
This should allow to cross correlate many of these experiments in this respect and gain a clear 
understanding of the relevance of the CGC in high energy experiments. It should be noted, 
however, that the tools are most developed for the simple case of deep inelastic scattering of 
leptons on protons and nuclei. This example then will also be my starting point to introduce the 
physical ideas -with a minimal amount of formulae- in the following section. 

1.2 The CGC in DIS: high energies, gluon densities and the saturation scale 

To acquaint the reader with the main features of the CGC, it is perhaps best to refer to the 
example in which the physics ingredients are best understood and the technical tools that were 
developed to incorporate them are most directly applicable, deep inelastic scattering of leptons 
on (preferably large) nuclei. I will attempt to keep the technical details in this section to a 
minimum. All of the relevant physics ideas collected here will have their theoretical foundations 
reviewed in later sections. 

To start with, consider a collision of a virtual photon that imparts a (large) spacelike mo- 
mentum q {Q^ := —q^ is large) on a nucleus of momentum p as shown in Fig^ Being spacelike, 
Q^ sets the transverse resolution scale in the problem. At large energies the only other relevant 
invariant is Bjorken-x (or Xbj), defined as Xbj := as in Fig.^ ^h] has the interpretation of 
a (total) momentum fraction carried by a struck parton. At small Xbj this has the interpreta- 
tion of a light cone momentum fraction fc+ = a^bjP^ where k is the momentum of the probed 
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constituent and p the target momentum.^ x^j directly corresponds to the large boost factor 
separating (photon) projectile and target (nucleus) via = as shown in Fig|2'' At high 




Figure 1: Kinematics in DIS of leptons off protons 
or nuclei 



Figure 2: Kinematic variables in a Minkowski di- 
agram. Projectile and target are separated by a 
large boost factor — — . := -k=(x'^ ± x'^) 

are the lightlike directions in the plane often 

referred to as "+" and "— " directions. 



transverse resolutions the situation is simplest. Many quantities, amongst them the total cross 
section, can be addressed in the framework of an expansion in inverse powers of the transverse 
resolution Q^. This is known as twist expansion and Operator Product Expansion (OPE). Tech- 
nicalities aside, this expansion has the character of a density expansion: at high resolution a 
nucleon in the target (which may be a proton in the simplest case) appears to be a dilute sys- 
tem of quarks and gluons. In such a dilute system and in lowest order of the expansion it is 
sufficient to express the cross section entirely in terms of two point functions, the quark and 
gluon distributions. Only higher orders of the expansion would involve particle correlations - 
which are naturally suppressed when densities a small. Nevertheless, even the individual terms 
in this "density expansion" receive quantum corrections which are important. At large the 
dominant ones carry factors of (a^lnQ^)" where a large logarithm (InQ^) compensates for the 
smallness of the coupling. These contributions need to be resummed to all orders to understand 
the dependence of quark and gluon distributions at high resolutions. This can be done dia- 
grammatically -known under the keyword 'resummation of ladder diagrams'- or by deriving a 
renormalization group equation, a differential equation for the distribution functions with respect 
to InQ^, the Dokshitzer-Gribov-Altarelli-Parisi (DGLAP) equations. By solving these equations 
one resums the logarithmically enhanced corrections. In fact, the differential equations are the 
only information that can be extracted perturbatively: their initial conditions, the quark and 
gluon distributions at some resolution Qq do contain nonperturbative information about the tar- 
get. As a consequence one needs experimental measurements over a range of to extract the 
gluon distributions at Qq by a fit consistent with further evolution. One gains quark and gluon 



^This is not to be confused with Feynman x or xp, defined as the longitudinal momentum fraction of a 
constituent. The two agree to lowest order in the parton model of DIS, but start do differ when pcrturbative 
emission into the inclusive final state are considered. Less inclusive measurements bring xp into its full right: 
longitudinal momentum fractions of measured particles come into the game and relate in specific ways to the 
longitudinal momentum fractions of incoming constituents giving rise to process specific definitions of xp. 

more precise definition would employ to the worldlines of constituent quarks and gluons in the projectile 
and target wavcfunctions which whill appear below. 
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distributions that can be used also in other experiments on the proton. A quite striking result of 
this procedure is a strong growth of the gluon distributions at small Xbj as sketched qualitatively 
in Fig. 13 This rise is driven by soft gluon radiation, (sea) quark distributions simply follow their 



. Q2^20GeV2 


i Q'- = 200GeV^ 


q2 = 5Qpv2 





IQ--^ W-^ IQ-- 10-' 



F = lnl/Xbj 

Figure 3: Qualitative growth of gluon distributions as from typical DGLAP fits of DIS at HERA (eA) 
for difi'erent Q^. 

rise. 

This procedure is self consistent and perturbatively under control as long as one stays within 
the large region. Starting from a large enough Qq and going to larger values, one finds that 
the objects one counts with the quark and gluon distributions increase in number but they stay 
dilute: as pointlike excitations their apparent sizes follow the transverse resolution scale 1/Q. 
At small Xbj, however, the growth of the gluon distribution is particularly pronounced - a clear 
sign that corrections of the form (a^ ln(l/a;bj))") start to become important. Any attempt to 
track the Xbj -dependence of the cross sections by summing these is immediately faced with an 
additional complication: moving towards small a;bj at fixed increases the number of gluons of 
apparent size 1/Q, so that the objects resolved will necessarily start to overlap. This is depicted 
in Fig.0J At this point one has long since left the region in which a treatment at leading order in a 
expansion (OPE at leading twist) and the concept of such an expansion itself is meaningful. 
Instead of single particle properties like distribution functions one needs to take into account all 
(higher twist) effects that are related to the high density situation in which the gluons in the 
target overlap. At small .Tbj one is thus faced with a situation in which one clearly needs to go 
beyond the standard tools of perturbative QCD usually based on a twist expansion. 

There are different ways to isolate the leading contributions to the cross section at small a;bj 
in a background of many or dense gluons and the history of the field is long [7-31]. I will use 
a physically quite intuitive picture first utilized in the McLerran-Venugopalan model which was 
originally given in the infinite momentum frame of the nuclear target. Knowing already that the 
interaction will be dominated by gluonic configurations, one tries to isolate the boost enhanced 
ones among them. Looking at the gluon field strength of the target for that purpose, one finds 
that only the components are boost enhanced, all others can be neglected. At the same time 
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Figure 4: Left: Partons in the Y,lnQ^ plane represented as dots of size Numbers grow both 

with and ln(l/a;bj)- Going to sufficiently small Xhj at fixed partons start to overlap. Right: 
corresponding regions of phase space with qualitatively different behavior. The lower boudary of the 
CGC region is determined by Qa{xhj)- The extendeded scaling region will be explaind in Sec. Ii-i.71 



one finds a strong Lorentz contraction in direction, and a time dilation, correspondingly, in 
direction. With such a field strength tensor, one is left with only one important degree of 
freedom, which, by choice of gauge, can be taken to be the + component of a gauge field. Taking 
into account time dilation and Lorentz contraction, the gauge field can then be written as 

A^b + SA with b'-^O, b+ ^ f]{x)S{x^) . (1) 

where /? is an unspecified function of the transverse components x := (a;^,a;^) of the four vector 
x. This exhibits a leading contribution 6+ which is x^ independent and Lorentz contracted to 
a 5-function in x~ . The 6 A symbolizes possible corrections that are kinematically suppressed. 
They will start to play a role as quantum corrections in the derivation of evolution equations. 
The boost and Lorentz contraction arguments for the leading contribution b are Xhj dependent. 
In the interpretation of Xbj as a lightcone momentum fraction it determines a resolution in 
x~ : the 5{x~) in is "localized" only with the resolution available with a probe providing the 
corresponding momentum fraction. A similar cautionary remark applies to the x^ independence. 

Note that mathematically one can always trade a component of a gauge field for a path- 
ordered exponential along the direction that picks up this component: 

,.*^.i3*U)U^ f/..Poxp{- = /..-»*(,-,x,0)}^ (2) 

If multiple eikonal interactions are relevant as the high energy nature of the process would 
suggest, one expects these path ordered exponentials to be the natural degrees of freedom, since 
this is what they encode: Diagrammatically 
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where the path ordered integrals along z~ are represented by the fermion line. Ignoring the small 
fluctuations bA for a moment, one has a picture in which the 7* nucleus cross section arises from 
a diagram in which the photon splits into a qq pair which then interacts with the background 
field. Due to the 8 like support of the qq are not deflected in the transverse direction during 
that interaction. This just reflects the largeness of the longitudinal momentum component of 
these partons at small Xbj. 



This is shown in Fig. 5(a) 



To be sure, the physics content is not 
frame dependent although it is encoded differently in, say the rest frame of the target. There, 
one encounters neither Lorentz contraction nor time dilation. However the scale relations are 
preserved: the photon splits into a qq pair far outside the target and its jT is so large that typical 
variations of the target are negligible during the interaction. As a consequence the probe 
is not deflected in the transverse direction, picking up any multiple interactions with (gluonic) 



scattering centers as it punches though the target. This is shown in Fig. 5(b) 



Lorentz contracted to )" ^ — : huge ~ A}'!'^' finite 

■''l>j 




(a) (b) 

Figure 5: qq pairs interacting witfi a nuclear target at small a^bj. (a) shows the situation in the targets 
infinite momentum frame with fully Lorentz-contracted target fields, (b) shows the situation in the target 
rest frame in which the qq vertex is far outside the target at a distance proportional to the relative 
boost factor l/xbj = e^. 

That multiple interactions are of relevance immediately becomes obvious, once one tries to 
calculate such diagrams with a background field method. To this end one evaluates the diagram 
shown in Fig. [S] in the background of a field of the type shown in Eq. ignoring the small 
fluctuations. Such a calculation will prove somewhat involved irrespective of the method chosen 
-one of them will be sketched in Sec. Eh but the resulting expression* 

aDis(:^bj,0')=Im™<CQ;>— = Jd\\^p'\{r^Q') j d% {^^^^—^^) (4) 

has a clearcut interpretation. Beginning with notation, r = x y corresponds to the transverse 
size of the qq dipole and b = {x + y)/2 to its impact parameter relative to the target.^ This 
result shows a clear separation into a wave function part (the factor \ijj'^\{r^Q-^)) and the dipole 
cross section part 

f tr(l-J7^J7t) 

(Jdipoic{Y = ln(l/xbj),r j d b ( — )Y=inii/x^,) ■ (5) 

The former has the interpretation as the absolute value squared of the qq part of the photon wave 
function, contains the ^*qq vertices and can, alternatively to the background field method, be 

^An integral over longitudinal momentum fractions of the q and q has been absorbed into the wave function 
part of this expression for clarity - the full expression is given in Eqns. 1181 and 1221 

defined like this is the Fourier conjugate variable to the transverse momentum transfer in non-forward 
matrix elements. Sec also IllOi . 
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calculated entirely within QED. The wave functions themselves turn out to be Bessel functions 
Kq^i depending on the polarization of the virtual photon. This factor carries all of the direct 
Q^-dependence of the cross section. 

The dipole cross section part embodies all the interaction with the background gluon fields in 
terms of the path ordered exponentials. In particular, in the absence of a gluon field representing 
the interaction with the target, both U factors reduce to unit matrices and the cross section 
vanishes. The averaging procedure (. . .)y contains all the information about the QCD action 
and the target wave functions that are relevant at small .Tbj . 

Even without exploring any of these in detail, there is one feature about this interaction 
which is very characteristic for DIS at small .Tbj which is already visible in the above expression: 
although momentum is exchanged between projectile and target, there is no net color exchange: 
the qq pair enters and leaves the interaction region as a color singlet. The same applies to other 
intrinsic quantum numbers such as spin. In a diction dating back to pre-QCD times this is often 
referred to as "pomeron exchange," although it is now acknowledged that the underlying theory 
of the interaction does not exhibit a particle cut that would allow an interpretation in terms of 
a simple particle exchange. 

It is fairly clear from the above discussion, that the isolation of the leading 6+ contribution 
which defines this average is a resolution dependent idea: as one lowers Xbj, additional modes 
-up to now contained in SA of Eq. (pj- will take on the features of 6+. They will Lorentz contract 
and their .T+-dependence will freeze. Accordingly the averaging procedure will have to change. 
If one writes the average as 



the weights WV[6"'"] or Zy[U] will, by necessity, be F = ln(l/a;bj) dependent. 

Besides the perspective adopted above, there are equivalent ways to understand and derive 
the Xbj dependence of such cross sections in which the path ordered exponentials do not appear 
in closed form in the first calculational step but instead are built up perturbativcly. Among 
those is the most straightforward procedure in which one allows the incoming qq (into which the 
incoming 7* splits initially) to branch off additional gluons before it hits the target. Looking 
at the corresponding amplitudes this would amount to consider perturbative corrections to the 
7* ~* QQ amplitude contained in Fig. of the form 



All but the last of these diagrams would be present already in QED. It is this last type of 
contribution that is responsible for the strong growth and the saturation features of the CGC. 
That all these methods lead to the same result will be thoroughly explained in the main text. 

Fig. El highlights another feature of the dominant contributions: as the distance between the 
7* — > qq vertex and the target grows with increasing energy like the corresponding boost factor, 
the probability of multiple gluon emission grows accordingly. This leads to the iterative scheme 
indicated in Fig. |H| (with only part of the above diagrams displayed for brevity). Note that the 
leading contributions come, as usual, from ordered emission,^ so that the lines depicted there all 

^Logarithmically enhanced contributions that lead to small x evolution equations arise from a region in phase 
space in which the x^j values in an n-gluon amplitude are strictly ordered x^^j-^ > Xf~,j-^ > . . .aj^j^. This is in 
complete analogy with transverse momentum ordering in the DGLAP evolution case. Ordering in x^j also implies 
a hierarchy of longitudinal distance scales as indicated in Fig. |H] 




or, cquivalently (• • ■)y 




(6) 




+ ... . (7) 
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target 




Figure 6: RG corrections to tlie average over background field configurations. Tfie sliaded areas contain 
tiie contributions subsumed into the averaging procedure at successive values of y = ln(l/a;bj)- 



have large momentum components in the photon direction. They themselves wiW then interact 
with the target in a similar way as the original qq component: These punch straight through 
the target at fixed transverse locations x and y, giving rise to the U^U^ in the expression for 
the cross section in Eq. I^J. The additional gluons shown in Eq. (0) and Fig. leave behind 
new gluonic Wilson lines U (the tilde indicating that they are in the adjoint representation) at 
new locations Zi. It is clear from gauge invariance that the second diagram in Fig. will involve 
the operator f/°^ 2tr(t°J7^t^[/^) where [7°^ represents the added gluon. From this perspective 
one would expect that calculating additional corrections would lead to an infinite hierarchy of 
coupled equations for more and more complicated correlators of Wilson lines. This indeed is the 
perspective taken by Balitsky in his derivation of this hierarchy of evolution equations [26] , the 
Balitsky hierarchy. The perspective taken above is different in the sense that these same diagrams 
have been interpreted to contribute a change in the averaging procedure that describes the dipole 
cross section of Eq. (O . The qq dipole at smaller Xbj is taken to interact with additional gluons 
and thus the cross section and averaging procedure changes. In Fig. this redefinition process 
is indicated by the shaded areas that include the target, and subsequently the target and all 
the gluons softer than the initial qq pair. Both approaches are indeed equivalent as will become 
apparent from the mathematical treatment. 

As these descriptions indicate, it is possible to calculate the change of the weights in Eq. ((HJ 
with F to a certain perturbative accuracy - at the moment complete calculations of these evo- 
lution equations based on independent techniques exist at leading log accuracy, i.e. to accuracy 
ofg ln(l/a::bj)- while Zy[U] at a given Yq remains incalculable without nonperturbative input. 
This is again the situation one has to face in any Renormalization Group (RG) setting in which 
the initial condition for evolution remains outside the scope of the calculation and must be 
determined from experiment. 

The renormalization group character of the calculation can be made explicit by assuming one 
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knew, say Zy[U] at an initial Yq. This defines the ensemble of initial background fields and 
one may then integrate over the fluctuations SA around between the old and new cutoffs Yq 
and Yi . Taking the limit 5Y = Yi — Yq ^ Q one then gets a renormalization group equation for 
Zy[U]. 

Because one integrates over SA in the background of arbitrary b of the form J^l, the equation 
treats the background field exactly to all orders and captures all nonlinear effects also in its 
interaction with the target. The resulting RG equation is a functional equation for Zy[U] that is 
nonlinear in [/. It sums corrections to the leading diagram shown in Eq. in which additional 
gluons are radiated off the initial qq pair as shown schematically in Fig. All multiple eikonal 
scatterings inside the target (the shaded areas) are accounted for. 

The result of this procedure is the JIMWLK equation alluded to above. This equation and 
its limiting cases will be discussed in detail in Sec^l 

Let me close here with a few phenomenological expectations for the key ingredient of the 
present phenomenological discussion, the dipole cross section. Clearly, one expects the underlying 
expectation value of the dipole operator, 

tr(l - UxUl) 

Ny xy := ( >F=ln(l/2;bj) : (8) 

to vanish outside the target, so that the impact parameter integral essentially samples the trans- 
verse target size and thus scales roughly with A^/'^. For central collisions, i.e. deep inside the 
target, NY,xy is expected to interpolate between zero at small r = x — y and one at large r. Note 
that Ny^xy is necessarily limited to lie inside these boundaries as an immediate consequence of 
the resummation of the gluon fields into eikonal factors U^^\ Saturation of these limits at small 
and large distances is the idea of color transparency for small dipoles supplemented with satura- 
tion for large dipoles and is sketched schematically in FigEI Focusing on the rightmost curve in 



1 
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Figure 7: Generic evolution trend for a one scale dipole correlator N^y plotted against |a; — y | as 2: — > 0: 
the curves move towards the left as Y = ln(l/xbj) increases, the asymptotic form being a step function 
at the origin. 

FigO which corresponds to Yq one has a single scale, Qs{Yo) — l/Rs{Yo), that characterizes the 
transition between transparent and saturated, hence its name, saturation scale. Looking back 
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at the definition of Ny xy in terms of the relevant fields U''^\ it is clear that the region where 
it differs from unity is simply the region where the fields are correlated: Rg = 1/Qs has the 
interpretation of a correlation length. From this it is apparent that evolution towards smaller 
Xhj at fixed (along vertical lines in FigEJ, which will add more gluons, will necessarily lead 
to shorter correlation length and thus larger saturation scales: The curve will move toward the 
left as indicated, with smaller and smaller Rs{y) characterizing the changeover between the two 
regimes. Since small Xbj evolution equations respect the generic saturation features shown in 
Fig. 13 Qs is also the scale that characterizes the onset of color recombination effects that slow 
down and eventually stop the growth of the dipole cross section for large dipole sizes. 

From here one can immediately understand that going from small to large hadronic targets 
(i.e. by increasing the target's atomic number A) one also increases the apparent gluon density 
in this kinematic domain: The starting point again is that the energy is assumed to be high 
enough that the projectile would punch straight through the target so that multiple scattering 
contributions sum into the eikonal factors entering Ny xy In large nuclei, one scatters off gluons 
from many independent nucleons along the qq trajectory as sketched in Fig. |H1 Following the 




Figure 8: The origin of growth of Qa with A in dilute targets with decorrelated scattering centers in 
longitudinal direction for a small and a large nuclear target. The lines indicate the qq pair, the dots the 
uncorrelated scattering centers. 

Lorcntz contraction argument from above, the whole depth of the nucleus contributes to the 
charge density seen at a point in the transverse plane. Since, due to confinement, gluons from 
different nucleons -encountered at different longitudinal positions- are necessarily decorrelated, 
this naturally leads to a rescaling of the transverse correlation length by the inverse nuclear 
radius. Since the nuclear volume scales with the atomic number A, this implies scaling of 
with A^/^: 

QfiYo)' ^ g?(>o)' ^1/3 . (9) 

This behavior will reappear in a model that was developed to describe qualitative features of 
saturation effects at small Xbj, the McLerran-Venugopalan (or MV) model, and explains the 
interest in studying DIS at small Xbj off large hadronic targets. The argument leading to Eq. ^ 
is simplistic and a more careful treatment will lead to a modification of the A dependence - 
although an enhancement typically persists. The most obvious source of such modifications are 
correlations induced by evolution towards smaller Xbj . This is why the scaling relation is written 
for some initial (not too small) Yo at which such effects can be expected to be small. 

Returning to the Xhj dependence sketched in Fig.[7| it is crucial to note that such behavior has 
far reaching consequences as it leads to perturbative consistence of the RG approach described 
above. If one takes the Y dependence of N sketched in Fig. as given, any RG equation 
that leads to this behavior will predict the largest changes at momentum scales of the order 
of Qs{Y) - this is where the infinitesimal change between neighboring curves and thus the 
r.h.s. of the corresponding evolution equation peaks. Contributions at momenta below Qs are 
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strongly suppressed. As Qs grows with Y, evolution moves away from Aqcd and thus remains 
perturbativc. This is illustrated in Fig.|51 which shows the curves of Fig.|7|in a 3d plot on the left 
and the scales active in the evolution equation i.e. dyNy on the right. Fig.lHlanticipates a further 




Figure 9: Dipole function (left) and its rapidity derivative (right) for central collisions plotted against 
dipole size. The latter shows that the change induced by evolution happens at larger and larger momen- 
tum scales. This implies IR safety. 

characteristic feature of evolution towards small Xhj that, at this point, can not be predicted, but 
was observed in DIS data by Golec-Biernat and Wiisthoff (G-B+W) [32-34], namely scaling with 
Qs{Y), a phenomenon called geometric scaling [34] by Stasto, Golec-Biernat, and Kwiecinski. 
They have found that for Xhj values below 10~^ all HERA data show beautiful scaling, if in 
Eq.ljHl) one assumes the scale in the r dependence carries at the same time all the Y — ln(l/j:bj)- 
dependence according to 

Ny ^y = N{{x-yfQAY)^) . (10) 

Such scaling takes all the curves represented in Fig. and maps them onto a single curve as 
shown in Fig.^l The same carries through to the DIS cross section via Eq. (@J) and leads to the 




rQs{Y) rQs{Y) 

Figure 10: Qa scaled version of Fig|^ assuming that geometric scaling holds. 

successful fit of the HERA data just alluded to. 

There are many more physics ideas associated with the label CGC, such as the existence of 
an extended scaling region at resolutions Q above Qs indicated already in Fig. 01 which is quite 
important for the conceptual consistency of the scaling fit to HERA data, the relevance of the 
CGC to the initial conditions of heavy ion collisions, consequences for particle multiplicities and 
the Cronin effect, as well as an intimate relationship with other QCD phenomena that involve soft 
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gluon emission such as jet physics. Among the core observations is the fact that the JIMWLK 
evolution equation intrinsically leads to the scaling behavior discussed in Figs. 1^ and 1101 

Many phenomenological applications are based on this scaling behavior and I strongly rec- 
ommend to have a first glance at Sec.|Hlto get an idea of the broad range of affected phenomena 
before delving into the main text. This should provide additional a strong motivation to under- 
stand the underlying theoretical tools and evolution equations. 

1.3 How to read this review 

Having given a brief overview over the core physics ideas in the previous subsection, it is left to 
the body of the paper to substantiate these claims and flesh out the context from which they 
arise. 

This review contains two interwoven threads of arguments, one of them centered around the 
derivation of evolution equations that govern soft gluon emission from QCD at small coupling 
but large gluonic densities, the other aimed at exploring their consequences and interpretation 
using models and numerical work as well as phenomenological applications. They can of course 
not be fully separated. Clearly the theoretical grounding is necessary for any application, but 
also the other direction is vital: for instance the numerical work in [35] has had consequences on 
the understanding of the necessity of running coupling corrections, not only in a quantitative, but 
also a conceptual sense. Similarly, the discovery of geometric scaling in HERA data has prompted 
a look for scaling features in the evolution equations. This in mind, the paper develops both 
threads in stages. It presents a different perspective on the derivation and interpretation of the 
underlying evolution equations than other existing reviews [36,37], not the least by highlighting 
the parallels with nonlinear effects in jet physics in Sec. |31 I have also been aiming at exposing 
the close connection between the McLerran Venugopalan model with the BK equation in Sec. 13. 51 
These formulations have not been shown elsewhere, but are hidden under the surface of many 
discussions of the Glauber-Mueller type nature of the multiple scattering phenomena omnipresent 
in small Xbj physics. I have attempted to present an up to date description of the role of 
BFKL physics as the driving force of small Xbj evolution and have included results of numerical 
simulations of the full JIMWLK equations as well as a thorough discussion of running coupling 
in the BK case. 

Sec. 12 is devoted to the first important type of nonlinearity encountered: eikonalization 
of soft gluon fields into path ordered exponentials often also called eikonal factors^ or Wilson 
lines. Sec l2.1l aims at explaining their appearance in cross sections and propagators and Scc l2.2l 
discusses an immediate consequence, boundedness of the underlying correlators, in terms of the 
MV model and its generalizations. The link to gluon distributions, density effects and A scaling 
also appear in this context. 

SecOlthen gives an overview over the formulation of the JIMWLK equation and the associated 
physics ideas. It starts off firmly in the formal thread with the full derivation in SeclHl- I consider 
this quite informative, but the reminder can be followed without having read this in a first go. 
Sec. 13.21 explains the JIMWLK Hamiltonian, its properties and the nature of the equation as a 
Fokkcr-Planck equation. Sec. 13.41 formulates the idea of scaling in the context of these evolution 
equations and Sec. 13.51 discusses the relation of JIMWLK, BK and BFKL (in both forward 
and nonforward cases) equations from the perspective of Nc and density expansions. How the 
underlying growth pattern is imprinted by BFKL physics is discussed in Sec. 13.61 This gives rise 
to the idea of a scaling window as emphasized in Sec. 13.71 and already shown in Fig. ^ as the 
extended scaling region. 

^This term has found various applications: besides this usage, also exponentiation of low order contributions 
in correlators like JSJ is often labelled as eikonalization. 
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Sec.0]performs a translation of the JIMWLK equation to a Langevin type formalism: Sec 14. II 
gives a toy example to illustrate the idea while Sec l4.2l givcs the result which is the basis for numer- 
ical simulations. Scc l4 . 31 introduces shower operators which will, in Sec.Elfind an interpretation 
as operators that generate the soft gluon cloud in the photon wavefunction and counts among 
the conceptual developments. 

Sec. pleads up to this interpretation and helps to understand that JIMWLK evolution shows a 
structure which is generic for soft gluon emission physics. This is illustrated with an example from 
the context of non-global jet observables for which Sec. l5.1l and l5.2l give evolution equations with 
clear structural analogies to BK and JIMWLK. Next is a derivation of these equations which may 
be skipped in a first reading. In this derivation one first constructs the leading part of generic soft 
n-gluon amplitudes (corresponding to the photon wave functions) in Sec. 15.41 As a second step 
one utilizes the iterative nature of the construction process to extract a corresponding evolution 
equation in the next subsections. Sec. l5.7l highlights the two types of nonlinearities encountered in 
the derivation -eikonalization and nonlinear evolution- and explores Nc or correlator factorization 
once more in the jet situation. 

Sec. 15.81 briefly discusses the lessons to be learned for small Xbj evolution from the preceding 
section and provides a first glance at what is to be expected if jets are created inside a dense 
medium. 

Sec. El and 13 discuss numerical results from solving JIMWLK and BK evolution equations. 
Discussed are scaling, 16.11 correlator factorization violations, 16.21 and the problems of fixed 
coupling simulations in the UV, 16.31 With running coupling UV phase space is under control 
and allows quantitative studies of Qs and the scaling features. 17.11 This subsection contains the 
core results on scaling behavior. Sec. 17.21 returns to A scaling already introduced in the MV 
model and discusses the changes imposed by evolution. 

Sec. IHl tries to give a flavor of the phenomenology associated with the previous theoretical 
results. Highlighted are geometric scaling, initial conditions of heavy ion collisions, consequences 
for the Cronin effect and a few more. 

Experimental efforts in connection with RHIC, LHC and EIC increase the pressure to come 
up with more observable consequences and better tools to make predictions - a few ideas are 
mentioned in the conclusions. Sec. El 

2 Eikonalization at high energies 

2.1 Propagators and cross sections 

There are two distinct types of nonlinearities present in the formulation of cross sections and 
their Xhj dependence as presented in Sec ll.2l the first is the resummation of the gluon field into 
path-ordered exponentials U^ ^ collinear to the projectile, the second is the nonlincarity of the 
evolution in correlators of this type of field. The first aspect is present already in the Abelian 
theory, the second appears only in the non- Abelian case but can not be separated from the first. 
Without the first step, the identification of the relevant variables, all the following steps have 
no proper foundation and formulae like Eq.Q remain mysterious. Let me, therefore, start with 
identifying the variables and a brief derivation of this formula for the DIS cross section. 

The starting point is the kinematics. A collision at large energies singles out a longitudinal 
direction, say x^, as the collision axis. In such a frame the large energy is characterized by 
a big rapidity separation between target and projectile. To judge how a projectile, say the 
qq dipole emerging at lowest order in the photon wave function of Eq. Q), sees the target, 
it is best to consider the color Field generated by the target. Remembering that light cone 
components a;* :— -75(2^" ± x^) are eigendirections of a boost in 3 direction {x^ x^ 
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with 7 = = . ^ ) it becomes clear that the components of the field strength tensor are 

enhanced by a factor in the plus components and suppressed by a factor e~^' in the minus 
components. At the same time Lorentz contraction "shrinks" the support of the color field to 
a (near) delta function in x" and time dilatation renders it (nearly) independent on . The 
leading components are 

F'+{x) =6{x-)d'l3{x) . (11) 

On the r.h.s. I have already made use of the gauge freedom to express via a single scalar 
function which can be associated with the 6+ of Eq. and (jSJ . 

To calculate cross sections like one needs to consider current correlators of the form 
{J^{u)J'^{v)) in the background field representing the target and subtract noninteracting coun- 
terpart (to go from an S-matrix to a T- matrix contribution). To lowest order, the cross section (@J 
then takes the form 

<7Dis(a;bj, Q^) = Im — ^^I^^^^^^J^^^' — ^ 

(12) 

where the interaction with the target is now explicitly trough the fcrmion propagators in a 
background field — ®^ — and the average (...) is over background fields b as mandated by the 
interaction with the target. The simplification lies in the fact that the types of field allowed to 
enter the interaction is restricted to the type in (|1HI . With it one can find results for cross sections 
and also include perturbative corrections that will modify (|12|l once quantum fiuctuations around 
b become large. 

The latter will be discussed in detail in later sections, here I am interested in the physics aspect 
of the calculation leading to (0J. The core clement here is the background field propagator which 
encodes all the interaction with the target. It is sufficient to study this propagator to understand 
the origin of the eikonal factors in the dipole cross section as well as the general form of the result. 
The key physics insight will be that the interaction is solely mediated through an exchange of 
momentum. 

Besides these physics results, gauge invariance is a somewhat delicate issue: While parametric 
enhancement and Lorentz contraction of components in the Field strength tensor can be argued 
in a gauge independent manner, the extraction of path-ordered exponentials that lead to formulae 
like Eq. appears to be gauge dependent: What happens in a gauge in which one forces the 
"-|-" component of the gauge field to zero? This will be bricfiy touched upon towards the end of 
this section. 

To explore the situation, consider the propagator of a scalar field in the background of some 
gauge field. ^ This quantity has also been used as the core ingredient to the gluon propagator 
in [13] in the derivation of the final form of the JIMWLK equation. Using first a Schwinger 
parameter representation and then standard (nonrelativistic) Feynman-pathintcgral techniques 
in which the Schwinger parameter takes the role of "time," one writes its propagator as 

^^{x,y) = J^ dse-^'^'[^l = dsj [dz]e-^od-^ Pe^p[-^gJ dz^b,{z)} (13) 

where the path-integral is over trajectories z, parametrized by k, that connect x (at k = 0) 
to y (at K ~ s). The expression at 6 = is the free propagator. Based on the kinematical 
arguments above, the background field is of the form given in |^: b^{z) — g^^S{z~)/3{z). This 

*This is in essence the simplified version of an argument used first in [26] from which the spin degrees of 
freedom are exorcized for clarity. 
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implies that only the "— " component of the trajectories are probed and only a dependence on 
the transverse coordinates z remains. Due to the shape one finds free propagation unless 

X and y are on different sides of the z~ = hyperplane, in which case a nontrivial contribution 
from the path-ordered exponential emerges with free propagation before and after. Since the field 
is concentrated at z~ = 0, one may deform the cikonal factors outside the z~ = hyperplane 
into straight line eikonals according to 



Pexpi^-igJ dz''6^(z)| ^ Pcxp| -i5 J dz^ b+ {0, z, z^)^ 



(14) 



and one has to distinguish 4 cases depending on which side of one places x~ and . For 
x~ > > y~ it is possible to extend the path to cover the whole line from, z — con to 2 + oon 
(n a lightlike vector in "— " direction). This is what has been denoted Uz throughout: 



z+oon 



J/^ := Pcxpj -i.g j dz- b+{Q,z,z-)^ 



z — oon 



(15) 



For x~ <Q < y~ one obtains U^, and for same side propagation -with x~ ,y~ > or x^ , y" < 0- 
one simply gets the unit matrix and thus the free propagator. Following the lines of [26] one 
eventually arrives at an explicit expression for the scalar propagator of the form 



D^[b] 



- [eix- - y-)e{p-) - e{y- - x-)e{-p-)\ jd^p'd\' 



X 

e-*(f'-''> Pexp{-zg J dz' b+ {0, z, z-)y'^' y 



(16) 



where some of the components of these momenta are interrelated according top'^ — ^ 2p- 



p 

2p- 



Q' ' = 2^ and p' 



q' — p ■ The massive version simply has on-shell "+" momenta shifted 
by rr? Ilp~ . An alternative calculation of this propagator using a spectral method can be found 
for example in [22]. 

This background field propagator can now be used in the diagrams of ()12|l: 



Im ( / d^xe"'"-^ / d^ye-^^-yq" 



D^[b] 



D^[b] 



{y, x) — (same at 6 = 0) ) (17) 



where the momentum flowing through the diagram is q, the vertices at x and y carry factors of the 

incoming (outgoing) momentum (q^) [For scalar "quarks" the current is J^^{u) = (j)^ (u) 9^'/'(u) 
instead of the J^{u) = ip(u)j^'4>{u) for fermions]. The remaining calculation is straightforward 
but tedious. I will only report the key features needed to understand the result: 



• The minus component of the momenta is conserved. It is easy to convince oneself by 
inspection of the theta functions that the corresponding loop momentum will be constrained 
to lie between and q~ . Parametrizing it as a dimensionless fraction of q~ one sees a 
fraction a flowing trough one of the lines while a fraction I — a enters the other. This 
momentum fraction integral remains present to the very end. 
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• Since all propagation is free unless x and y lie on different sides of the z =0 hyperplane, 
all "same side" contributions cancel between the two terms of H17() . 

• In the remaining contributions one encounters eikonal factors of the form H15|l in the com- 
bination trUl_^Uz2/Nc. The second term at 6 = will simply subtract a 1 since all other 
ingredients are equal. The trace originates from the contraction of color indices at the 
j*qq vertex. This is the origin of the dipole operator tr(l — UziUl^)/Nc and allows the key 
observation that no color is exchanged with the target. 

• All integrals besides those over the transverse coordinates Zi, Z2 can be performed. The 
momentum exchange with the target is parametrized by this coordinate dependence. 

• The remaining x~ and integrals have the form dx~ dy~ or /^^^ dx~ dy~ . 
These cover regions outside the z~ = hyperplane where the qq pair propagates freely. 
Each x~ and y~ integral can be identified as an integral representation of a modified Bessel 
or McDonald function. They describe the amplitudes for the 7* qq transitions before 
and after the interaction with the target - the "virtual photon wavefunctions." [See [38] 
for a direct calculation of these amplitudes.] 



The final result, reads 



tr(l-f/.,C/t )^ 



<JmsiY,Q') = / d'r j da [V-K r^ Q2)|2 / ^'b C ^^'~'^'' ) ■ (18) 







For longitudinal polarizations (these are chosen for closer similarity to the fermion results with 
physical transverse polarizations) and vanishing masses one finds 



lVL(a, r2, Q2)|2 ^ ^-^a{l - a)Q^Kl{^a{l - a)Q^) . (19) 

The physically relevant case of spin 1/2 quarks can be derived following the same lines and 
leads to a very similar result - it is essentially only the Bessel functions that are modified in the 
end (see Eq. (|22|l 'l. However, starting the calculation with the fermion propagator in the small x 
background field adds technical complexity to intermediate stages of the calculation [26] . Since 
they do not influence the core physics content I will only briefly sketch what has to be done in 
addition the steps outlined above. 

One first has to derive the fermion propagator starting from the equivalent to the l.h.s. of (|13|l 

S[A]{x,y) = [i$ + m],,{-D^'+ga^'''F^,-m'y\x,v) (20) 
[a^'^ (X [7'',7''] generates spinor rotations) and ends up with the expression 

S[A]{x,y)=~i[ill)^m], j ^^—[e{x- - y-)0{p-) - 6{y- - x-)6{-p-)\ jd^p'd\' 

X~ 

_0_e-(pW). Pexp{-z,|d.(^A+ + ^)(0,.,^-(.))}e-- 



(21) 

with exactly the same restrictions on momenta as in the scalar case, Eq. (|16|l . In fact the uF 
contribution simplifies considerably, if one takes into account that {ctFY ~ due to cr^^a^^ = 
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(for arbitrary transverse indices i and j). The remaining contribution does not enter the DIS cross 
section at all. The final result is again of the form H18() although polynomials in the numerator 
lead to different types of Bessel functions. The wavefunction contribution now is given by 



6ao 




|^T.L(a,r,02)P ^ ^ _ (22) 

/ [4Q^a^l-a)^K^iQr) . 



J^k (23) 



where subscripts T and L refer to transverse and longitudinal polarizations of the photon and 
= a{l-a)Q^ + mj. 

There exist many alternative derivations in the literature. [39-41] for example build up the 
wavefunctions perturbatively and then calculate the expression for the cross section. In these 
approaches eikonalization -the appearance of the path ordered exponentials- has to be obtained 
perturbatively from diagrams like the ones in Eq. The key observation there is that in the 
kinematic situation at small Xbj , in which, say, a quark, with high energy of momentum p that 
emits a soft gluon of momentum fc, the product of the propagator of the emitting particle and 
the emission vertex, simplifies drastically and takes the form of what is often called an eikonal 
current 

p.k 

to which the soft gluon couples as JpfJ3^{k). With p pointing in "—"-direction this object involves 
a fc"*" and a only, and after Fourier transformation of the expression, one begins to recognize 
that the denominator induces the ^-function structure needed to build up the path-ordering in 
the expressions above. Iterating the procedure carefully leads to Eq. ((3J. This correspondence 
becomes important again when one starts to write down functional expressions that implement 
soft gluon clouds in both small Xbj and jet physics, the so called shower operators first encountered 
here in connection with a Langevin rewrite of the JIMWLK equation. 

All of these approaches share a common technical difficulty: as soon as one tries to choose 
a gauge in which the component (the component entering the path ordered exponentials 
Uz) vanishes, the above arguments appear to fail pathologically. Since such a choice has been 
used in the original formulation of the McLerran-Venugopalan model and in early versions of the 
JIMWLK equation, I need to explain how this confusion is resolved. Such a calculation shifts the 
leading contributions into transverse components of the gauge field bi{x) where the bar serves 
to characterize a gauge in which 6+ = 0. One then proceeds to show that there is only a single 
degree of freedom contained in the two components of hi by showing that they can always be 
written in terms of a group valued field U^- .x via 

This expression is of the form 6{x~)f3i{x), it has support over "half" of space time instead of 
only on a hyperplane at x~ =0. Only by comparing to the calculation presented above does one 
realize that C/^j-.x is nothing but the gauge transform that relates the two (barred and unbarred) 
gauges 

= Pexpj J dz^ b+{Q,x,z^)^ , (25) 

— oo 

although within the barred formulation, such an interpretation is not available. As was to be 
expected, one merely ends up reshuffling the same degree of freedom by changing the gauge. 
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Nevertheless, the calculation of the propagators via path-integral representations l|13() becomes 
quite nontrivial due to the altered support of the background field and a method based on 
spectral representations and wave functions in the presence of these background fields becomes 
more efficient. 

2.2 The McLerran-Venugopalan model and its generalizations 

The McLerran-Venugopalan (or MV) model was formulated as a means to describe cross sections 
at small but fixed Xbj, like the DIS cross section in Eq. Q), by parametrizing the averaging 
procedure on the r.h.s. with a suitable ansatz. In essence, one assumes that the dominant 
configurations b of Eq. ^ are governed by a Gaussian weight. This was argued to be reliable 
in the case of large nuclei and formulated in terms of color sources in the large hadronic target. 
The large field b of the introduction can be thought of as generated from these sources via the 
Yang-Mills equations. Just as the background field in the gauge chosen in the introduction, the 
corresponding color current is of the form 

J''(x) = 6''+Six-)p{x) (26) 

where p{x) = p'^{x)t°' is a matrix in color space that describes the color charge density in the 
target. As explained in the introduction, at large enough to be in the perturbativc domain, 
the probe can transvcrsally resolve colored structures inside individual constituent nucleons. In 
longitudinal direction the situation is different. Just as with the color field b created from it, the 
S{x~) structure of the current is to be taken with respect to a longitudinal resolution imposed by 
the value of Xhj in the experiment. In going beyond that resolution, or by viewing the situation 
in the target rest frame, one should think of p{x) as the integral in x~ of color charges at a given 
transverse position x: p{x) :~ J dx^p{x^ ,x). A probe with a;bj so small that it can not resolve 
longitudinal internal structures will couple to this integral directly instead of the individual 
color charges inside the constituent nucleons. Since the latter arc color neutral on their own, 
the target sees a incoherent superposition of color charges. Let me denote the incoherent sum 
of color charges in a tube extending the full length of the nucleus with transverse area l/Q^ 
(according to transverse resolution) Q° (x) . This may be written in terms of p as 



Q"(a;) := 




where (/)q2 represents a coarse graining function adapted to the transverse resolution scale 1 / 
(by its normalization J (Py (f)Q2{y) ~ l/Q^) whose precise nature will not be important. If 
we ignore geometrical complications and assume uniform longitudinal thickness ("cylindrical 
nuclei" ) should obey (locally, in the coarse grained sense) 

(S") = (Q^'Q") = 7^-^ . (28) 

The average here can be thought of as a configuration average. The first of these equalities states 
color neutrality and the second gives the "typical charge" squared. With Ra ~ -Rproton^^^'^ 
both the individual charge of Ea. H27|l and the correlator in Ea. H28() scale with A^^^, so that 
asymptotically [Q", Q'^] = if'^^'^Q''- ^ and the charges can be treated as commuting objects. 
This is the original McLerran-Venugopalan argument that would lead to correlators of color 
charges (Q""^ . . . Q"^) that are fully determined by the two point function and hence by a Gaussian 
functional weight. This weight is characterized by a width pj^, which should be thought of as Xhj 
and dependent. It is common use to not explicitly discuss the coarse graining in transverse 
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space and write this distribution in terms of p{x). At this point one should caution that although 
the x~ integral in H27|l will scale with A^/'^. the actual proportionality constant can only be found 
in a more careful treatment. As indicated, geometry alone leads to an additional factor that is 
gencrically smaller than 1. At a given impact parameter b, this could be taken into account 
by writing nR^SAib) instead of A^^^ [42]. Here SA{b) is the Wood-Saxon formfactor and R is 
the gluon radius of the nucleon, the size of the gluon distribution in the transverse plane. This 
would still imply scaling of Ql with A^^^ at large enough nuclei but would lead to a noticeable 
reduction for small nuclei. The naive approximation of a cylindrical nucleus on the other hand 
would amount to Qf{xhi) = Qf (.Tbj)^A^/^. 

A slightly different treatment of color charges has been adopted in the calculation of gluon 
distributions in [7] and in the derivation of the BK equation by Kovchcgov [30] . Here one would 
resolve the nuclcons in the nucleus and treat the color charges at different x~ as belonging 
to different nuclcons. Again one can ignore all higher order correlators of p and describe the 
p-distribution by a distribution local in transverse and longitudinal space 

Wuvip] ■■= A^exp { - i 1 dx-d'x . (29) 

The link between the two descriptions is provided by the requirement that the widths in the 
distributions be related by /i^ = J dx~ p,^{x~) wherever this integral is large enough for the 
original arguments leading to a Gaussian weight to apply. The first formulation directly uses a 
large surface density of gluons (as parametrized by p.\)., the latter formulation explicitly builds 
it from a locally small "per nucleon" volume-density p^^x"). Instead of the original classicality 
argument it is the dccorrclation in x^ which implies that charge commutators never play a role 
in the evaluation of correlators like the dipole function and thus that the Gaussian weight (|29|) is 
adequate. The latter description is actually capable of interpolating between two opposing limits 
as one increases p^{x~), say as a consequence of soft gluon radiation and coarse graining in x^ 
as one progresses towards smaller Xhy Both versions can be found under the name McLerran- 
Venugopalan model. A reliable treatment of the region in between the extremes where individual 
p?'(x~) grow large requires additional information that lies beyond this level of modeling. Infor- 
mation of this type is expected to arise from small Xbj evolution equations such as the JIMWLK 
equation. It is because of the link with one of the original derivations of the BK equation that I 
will further discuss H29|l and its generalizations. 

Note that due to its local nature and the gauge transformation properties of the current 
J, H29(l is a gauge invariant choice. The Gaussian leads to a local correlator for p{x^ , x) both in 
x~ and X where p'^{x~) turns out to be related to the gluon density. 

The above uses the case of large nuclei in which individual nuclcons are used as decorrelated 
sources of a cumulatively large number of gluons, but the resulting description in terms of a 
Gaussian weight may also apply to small nuclei or even a proton in a situation where evolution 
towards small x^,] builds up sizeable gluon fields, provided these fields are not strongly correlated 
over a large Xbj range. [A dependence of this alternative source of gluons would of course be 
different, see also Sec. 17.21 ] In this sense one may substitute small Xbj for large A in the above 
reasoning. 

It is instructive to calculate dipole functions and gluon distributions in the McLerran-Venugo- 
palan model which allows one to formulate a number of (possibly simplistic) phenomenological 
expectations that should be closer to reality the larger A and the smaller Xbj become. The 
simplest objects of interest are the dipole cross section at fixed Sbj, which contains the operator 
trUxU^ and, somewhat more complicated, the gluon distribution, which contains the operator 
tT{UdiU'^)x{UdiW)y. Since most correlators of interest -as well as the evolution equations 
considered below- can be expressed via path ordered exponentials which depend on 5+ instead 
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of p with —didib~^ = p (from the + component of the Yang-Mills equation), it is natural to 
replace p by &^ and write the weight of the McLerran-Venugopalan model as 



/i2(x-,Q2) 

It is clear that the 6-correlator remains local in x~ but becomes nonlocal in the transverse 
direction. Still, gauge invariance is guaranteed by its equivalence to the original H29|l . 

I will explain below that the local nature in x~ -which leads to the absence of commutator 
contributions in simple correlators- is, if not equivalent to, then at least compatible with the large 
Nc limit and is an important structural feature of the model. The specific nature of transverse 
nonlocality will be changed through gluon emission as induced by small x evolution be it in 
the JIMWLK or BK evolution equations. I will anticipate this and allow distributions with 
correlators of the form 

with a more general transverse coordinate dependence. Such a generalization has also been 
argued for by Mahlon and Lam [43] in order to implement overall color neutrality by imposing 
conditions on G^-'xy that were incompatible with the original MV form.^ 

For clarity of interpretation it is useful to change notation still a bit further and replace 
x~ in the above by Y itself. This is possible, since evolution towards smaller Xbj (larger Y) is 
related to a change of resolution in x~ as discussed above. A more detailed discussion of this is 
given in [14, 15], but an early version can already be found in [7]. To summarize, one ends up 
interpreting the U factors as path ordered exponentials 



Ux = 



P exp { - ig jdz-b+_ ^ } = : P cxp { - i JdY ay.^ } (32) 



where a can be taken to be defined via this change of variables. The above implies a Gaussian 
weight for a of the form 

Y 

Wria] exp ( - i jdY' jd\d\ a^,^GyU"?".0 (^3) 

with 

{aYa^ay.y) = 5{Y' - Y")Gy"^...B{Y - Y') (34) 

Simple correlators are calculable, for instance the (singlet) two point functions with fields in the 
fundamental and adjoint representations come out to be 

{i^iU^uD) =N, e-^'^^-" (35a) 
{HU^Ul)) ='2NcC! e-^=e-- (35b) 



^In this generalization gauge invariance still holds provided the x^j -dependence is derived using one of these 
evolution equations - they resum the gauge invariant, leading x^j-dependence and the initial condition is gauge 
invariant. An illustration that some apparent gauge dependence may be spurious at small x^j is provided by 
the dipole cross section itself: This is gauge invariant in the sense that the observable predominantly probes 
configurations b as in 0. Hence, it does not matter along which curves the Wilson lines are closed at infinity, 
any curve provides the same factor 1 which is not displayed in the formulae - this information is not probed 
(suppressed) at small x^j and the result is gauge invariant "to leading order at small x^j." The same reasoning 
applies to the evolution equations themselves. 
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with 



Y 

Qv.xy j dY' (^Y',xy — 



Gyi 



Y',xx 



Gy' 



y',yy 



(36) 



a symmetric function of x and y that vanishes in the local limit. The color factors arc generic: 
the prefactor is simply the dimension of the representation, the factor in the exponent is the first 
Casimir. 

Since such Glauber-Mueller type exponentiations are very characteristic for multiple scatter- 
ing events in the CGC I will, once, outline how to arrive at this result. This will also allow me 
to highlight the role played by color flow in this type of event. I will start with a diagrammatic 
interpretation of the lowest order contribution to, say (|35a|l . One finds 



N,C{ dY' GY',xy - 



G 



Y' ,XX 



G 



Y',yy 



+ 



(37) 



The horizontal lines correspond to the Wilson lines expanded to the order indicated by the 
number of giuon insertions. This implies the relative factor of — ^ of the last two diagrams 
compared to the first as shown on the l.h.s.. The hatched blob stands for G and the lower 
legs indicate that these correlators represent interaction with the target. The giuons hook into 
the Wilson lines at Y' , which is then integrated over up to Y, the value characterizing the 
functional weight l|33() . The external closing lines "(^" and stand for Kroneckcr deltas on 
the external color indices and implement the color trace. The path ordered exponentials furnish 
two Y integrals per factor of G and since the Y structure of the weight is local only one of 
them remains. It turns out that the contributions to all orders can be rearranged into a locally 
subtracted version of the two point function called Q in the above. Explicitly displaying the Y 
integrals, the relevant precursor to Q, which still contains nontrivial color structure is 



Y' 





(38) 



On the r.h.s. the color structure of the first digram is t° (8) t"- while that of the other two is 
simply Cr ■ \ ® 1 with Cr the first Casimir of the representation in question. To arrive at a 
generic expression for higher orders in G, locality in Y is essential. It implies that the n-th order 
contribution is a simple iteration of this structure: 



Y 



Y„, 



dY" 



(39) 



Note that the color exchange with the target is zero in each individual factor. Consequently an 
insertion of H38(l will never mix inequivalent multiplcts in the multiplet decomposition of Wilson 
lines considered. For the correlators in (|35|) this means that it is sufficient to calculate the singlet 
channel contribution of an individual insertion and iterate that to all orders. Since there is 



^"A closer examination of tiic multiplet structures would reveal the following; For the qq case an insertion 
of I38i can be written as a diagonal 2x2 matrix in a space consisting of a singlet and an octet. For the gg case 
the matrix operates on8®8 = lffi8ffi8©10©10©27 (in SU(3)) and will indeed mix the two equivalent adjoint 
representations therein. 
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only one singlet, the multiplet of interest completely decouples from all others. Writing the 
singlet projection operator generically as ^^^)C for ^ 2- Wilson-line example in a representation 
R with dimension d{R) , one finds that the color structures of all three terms in l)39|l in the singlet 
projection are reduced to Cr -j^^C so that one identifies 



:j7mX = - ( Gy,.^y - + ) Cn -—X (40) 



as the relevant building block for iteration in the singlet channel of l|39|l . Extending the inte- 
gration range of all integrals up to Y then restores powers of Q and furnishes the 1/n! factors 
needed to arrive at the simple exponential form shown in (|35|l . The external color trace acts on 
the overall singlet projector and provides the d{R) in front of the exponential factor. 

The fact that there is no color transferred to the target means that in terms of color, an 
individual insertion is planar. Locality in Y ensures that this remains the case for multiple 
insertions: all the diagrams entering the above correlators arc planar in their color structure and 
in this sense correspond to a large Nc limit [44]. 

Quite generically, one expects correlators of Wilson lines to show some exponentiating prop- 
erties in that their logarithm is in some sense generically simpler than the object itself, just as 
in H35|l - but in general without any reference to Gaussian weights or the Nc limit. Such examples 
are provided by general calculations of Sudakov type form factors, Drell-Yan production [45] or 
certain jet observables. The latter case will be discussed in Sec. and the simplification there 
is in terms of logarithmically enhanced contributions as stated in H151|l . In general the natural- 
ness of such a parametrization gives no clue as to whether it should arise from simple averaging 
procedures as employed here. Their use in this case is mainly motivated by the MV model and 
an application to the evolution equations with its resulting simplifications and interpretation to 
be discussed in Sec. 13.51 The comparison with the Sudakov example should provide enough of a 
hint as to how to think about extending such a picture at least in principle. 

Returning to the planar diagrams encountered with MV type models one can start to prepare 
for better phenomenological understanding. To this end consider first the gluon distribution 
which was originally calculated in [7] in the McLerran Venugopalan model in a somewhat more 
pedestrian manner. 

Since we deal with a gauge theory the definition of the gluon distribution is not simply 
formulated as a naive number operator construction. It involves a bilinear in field strength 
operators instead of a simple bilinear of fields. In = gauge, the expression can be rendered 
as 

XbjG(xbj,02) = IJ ^e(Q2-fc2)(F^+Kj,fc)FrKj,-fc)) . (41) 
This is written in a form that exhibits its relationship with the unintegrated gluon density 

' iV2 - 1 ttR^ dYd^k ~ ttR^ 7V2 - 1 ■ ^^^> 

In a general gauge one has to Fourier transform the expressions back to coordinate space and 
insert a Wilson line in x" direction. In the same = gauge, the Fourier transform of the 
FF correlator of the unintegrated Gluon density involves the operator tr{UdiW)x{UdiW)y (c.f. 
Eq.|2il)) which is perfectly calculable in this class of MV like models. One may even consider 

{[UmrjUd,U%) ^ ^ 5^' ^^]^^g^':f U e-^^^-A . (43) 
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It is clearly Q and its gradients that determine the size of the gluon distribution. Many phe- 
nomenological applications exist that fall into this class of parametrization and it is worth to 
give a first flavor of the issues discussed phenomenologically by listing some of them. The 
original MV expression of [7] emerges by setting GY;xy = g'^xO^:Q^)[l{x ^ V) ^ 7(0)]j where 
l{x) = -^{x) = ^£c2ln(a;2A^cD) + 7(0) and x{Y,Q^) ■= ^ dY' ^j?{Y' ,Q^) is the total color 
charge sampled inside the target at a given transverse resolution Q^. The gluon distribution 
proper emerges in either case after traces over color and Lorentz-indices are taken. The MV 
result reads (see [7]) 



mi - 1) 



NcX 



2 



l-ix'K 



2a2 X^X{Y..Q^ 



qcd; 



and leads to a unintegrated gluon distribution of the form 

In this expression x(F, Q^), which parametrizes the color charge in longitudinal direction, has 
been reinterpreted in terms of a Golec-Biernat+Wiisthoff type saturation scale Qg. Equally often 
one finds Qs in this expression identified with the gluon density of individual nucleons in the 
nuclear target [46-49]: 

Ql{y. -^ b) ^ x.,G(.., 1/.^)^ . (45) 

This allows for additional impact parameter (b-) dependence and introduces G{xh],k'^), the 
nucleon gluon distribution (with the identification = 1/x'^) and a nuclear profile function p 
that counts participants at fixed impact parameter. Often the gluon distribution in the nucleon 
is taken to be of the simple perturbative form compatible with the logarithmic structure of 7 in 
the MV model. Another model for Q to be listed in this context is the celebrated Golec-Biernat 
Wiisthoff model used to fit the HERA data. Accordingly this was used not to describe gluon 
densities but dipole correlators via 

i'^^1^3^),^l^,--^Qliry^ (46) 

in the same spirit as H44|l , but with a parametrization of the Y dependence added to cope with 
Xhj dependence in the data. All these interpretations share the same A^^^ scaling present in the 
MV model. Here is a collection of generic features: 

• Dipole correlators of the form tr(l — UxU^)/Nc will naturally be bounded by 1 at large 
distances and will show the color transparency + saturation asymptotics as soon as G grows 
with distance. 

• Any growth of Q with Y will result in the qualitative behavior sketched in Fig. [7| 

• As long as the assumption of uncorrelated scattering centers holds, one would expect that 
G scales like A^^^. This clearly will enhance the importance of the nonlinearities for large 
nuclei. Going to large nuclei in this sense has a similar effect as going to small Xbj . Quantum 
corrections that drive the change towards small .Xbj, however, will induce correlations and 
this will not occur when going to larger A. The consequent change in the naive A^^^ scaling 
will be discussed in Sec. 17.21 



^See [50—52] for studies of impact parameter dependence which remain largely outside the scope of this review. 
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The first two of these are consistency requirements, the last property, scaling of the exponent 
has been used to reinterpret the expression in many different ways as already indicated. This 
aspect connects to the realm of model building. 

While any reference to small Xbj and A scaling arc specific to the problem at hand, the 
exponentiation of leading order contributions in expressions (|35|l and l|43|l is not. Already in 
QED, in the calculation of soft photon bremsstrahlung, occurs exponentiation of this type. The 
interpretation of the object in the exponent there is that of a probability to not emit soft photons 
below a certain experimental resolution. Other calculations involving soft gauge bosons will 
acquire a similar form, in particular in situations in which the second nonlinearity mentioned at 
the beginning of the previous section, the creation of additional "hard" particles, is excluded. 
An example arc contributions to jet obscrvablcs from soft gluons going into the "empty" region 
outside the hard jets (see SecEl). 

The first decisive step beyond this stage is the derivation of an evolution equation that 
determines the Xbj or Y dependence of Q and its generalizations. This is the topic of the next 
section. 



3 JIMWLK evolution and the Balitsky hierarchy 

After having identified the relevant variables for a description of QCD scattering at high energies 
-which constitutes the first resummation, the eikonalization of gluons into Wilson lines U~ one 
now needs to calculate how general correlators of such fields change with Xbj ■ This leads to the 
second resummation necessary, this time in the guise of an RG equation. Sec. 13. l| gives a brief 
overview on how this calculation can be most efficiently organized and can be skipped by those 
only interested in a discussion of the results, which will be given in the subsections following it. 

3.1 A systematic derivation 

From the discussion above it is clear that one needs to understand perturbative, logarithmically 
enhanced corrections to correlators of Wilson lines of the general form (Ux^ . . . C5i Ux})^ with 
C/'^'f) in the fundamental representation.^^ 

One way to do this efficiently is to introduce a generating functional for such correlators and 
perform the calculation directly for this general case. I therefore introduce 

J] (e■5ext[^•^^'^l)^ (47) 

where 

SZ[A,J\J] = Jd'x{tr{iJlYUx[A+])+tiij',ul[A+])} (48) 
is an external source term, correlators are extracted via j'^^-* derivatives in the usual way 

6J ' 6J^ 

or explicitly 

-^■■■^^■■■^Z[JlJ] = {ul®■■■®Ul®U^®...®Urn), . (50) 



^^This is completely general, as one may write any higher representation as a (local) product of U^^K For 
example adjoint links emerge as a combination of two fundamental ones by virtue of U"-'' = 2tT[t"-Ut^U'']. 
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2 UUy 



2 obx oOy 



So far the physics of small Xbj is encoded in the type of correlation functions considered 
-exclusively correlators of link operators U^^^- all the rest is mathematical convenience to help 
summarize the result. To extract the logarithmic corrections, one now expands the gluon field 
A around b (whose correlators are assumed to be known) and keeps fluctuations to the order 
as ln(l/a;bj)- That is to say that one will expand around b to one loop accuracy and select the 
terms carrying a ln(l/xbj) factor. This way one will be able to infer the change of correlation 
functions as one lowers Xbj . This is the second statement about physics or rather what one can 
learn about it through such an approach. 

Now turn back to the actual calculation. At one loop one needs at most second order in 
fluctuations:^'^ 

(51) 

+ l{{SAJAy)sA[b] (2 l^-l^SZib, J\ J]) iJ^flAb. J\ J]) + ^^'^eftl^ J\J\) 
xe^'^t[^'^'''^l), + ... . 

The second term is the calculable perturbative correction to Z[j\ J], the generating functional 
for all generalized distribution functions of the target, an object which can not be calculated as 
such with present tools. 

To understand the individual terms in Eq. H51() in detail requires the calculation of the one 
loop corrections to path ordered exponentials [/[foj^^^-* in the presence of a background field b of 
the form Eq. Before diving into this it may be helpful to specialize once more to the DIS 
example from above to illustrate the physics content of the terms in Eq. 1)51(1 . Here one needs to 
look at 14 

tr(l - JJ^JJ-) (e'5-'[''+^^'^'"^l)^ = (tr(l - U^[b] C/^[6]))^ + quantum corrections . 



(52) 



The quantum corrections are induced by fluctuations 5 A incorporating physics below the Xbj 
values for which the original b were a good approximation. One would like to use these quantum 
corrections to redefine the h average to adapt to the new, lower values of Xbj . If this is possible 
the averaging procedure (. . becomes Xhj dependent. Of course one needs to perform this step 
in general for the whole generating functional, not only the specific correlator (|52|l . This means 
one has to carry out this step for the second term in Eq. ((51|l . This sets the task: // one is able 
to calculate these quantum corrections before taking the b average, that is, for all relevant b, one 
can deduce how the U^^^ correlation functions or the weight Z[U,W] that defines (. . .)(, evolves 
with Xbj- 



^^In writing Eq. 1511 one lias anticipated tliat (^SA^") g^[b] = as in the free case. This is a consequence of the 
structure of the propagator in this background field and fias been used repeatedly [12,26]. 

After integration over impact parameter (simply x + y, if one puts the target at the origin in transverse space) 
this is the dipole cross section that features prominently in calculations of 7*p or 7*^ cross sections at small 
Xbj [38,41,53]. 
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Clearly there are two generic types of corrections corresponding to the two terms in 



-{SA,JAy)sA[b] 2 j-SZ[b, J\ J] TrSllAh, J\J] + TT^S^^itlb. J\ J 



Sb„ 



5bx Shy 



(53) 



The first term represents contributions where the gluon propagator connects two quarks (Us) 
or antiquaries (U^s) as well as a quark to an antiquark. In addition there are pure self energy 
corrections dressing one quark or antiquark line instead of connecting two of them, represented 
by the second term. 

Particular correlators are again selected taking any number of J derivatives and setting J 
to zero. One point functions get contributions from the second term only, while both terms 
contribute to anything with more than two J^^^ derivatives. Clearly it is sufficient to calculate 



1 S 



2 S4'^ 



1 



J(t)=o 



-{SA^6A.,)sA[b] 2 —SZ[b,J^,J] —SZ[b,.P,J] 



Sb. 



Sb^ 



and 



SJ. 



(t) 



,/(t)=o 



t^SAJA,,)sA[b](^^J'jAb,JlJ 



(54) 



(55) 



to completely reconstruct Eq. 1)53(1 as these terms are precisely second respectively first order in 
J^''\ Eq. H53|) then defines the change in all other correlators and one has reached the goal of 
finding the fluctuation induced corrections. 

The task is clear and the terms appearing in the actual calculation are best visualized dia- 
grammatically. For the gluon exchange diagrams corresponding to Eq. (|54|l one defines 



ln(I/a;bj) -X 




ln(l/xbj) • Xl 




(56) 



and analogously for xl" 



and xll- 



vxy Axj/- already seen in [12,26], these split naturally into 
ordered contributions when one combines the structure of the vertices (the b derivatives of S'cxt 
in Eqns. H54() and H55() ) and the gluon propagator in the background field. Take x'' ^ a-n 
example: 




^Representations for Xxy Xxy s-^d Xxy result from reversing the quark lines aceordingly. 
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The quark self energy correction in Eq. (|55l) and its x ordered parts are given by 



ln(l/xbj) al 




(58) 



the antiquark diagrams are completely analogous. Again, the various components of x correspond 
to gluon exchange diagrams while the trs are self energy corrections. As for the time ordered 
ones, there are those with and without interaction with the target (with the gluon line crossing 
the x+ axis or not respectively as shown in Eqns. (|57|) and ^b^\ ). These latter versions determine 
the JJ*-^^ dependence of % and a in detail. They have been calculated several times using different 
methods. The common feature of all of them is that diagrams without target interaction leave 
the number of C/^^^^s invariant, while those with target interaction insert an additional adjoint 
J7^^^ at the point of interaction. This implies a nonlinearity in the evolution: Any correlator 
with a fixed number of fundamental representation C/^^^s will couple to other correlators with 
two additional factors of fundamental representation t/^^^'s in each step in I'bj. 

The diagrams have been calculated in [26] or [12], here I only quote the results for x a. 
For convenience, define the integral kernel 



{x - z).{z - y) 



{x-zY{y-zf 



(59) 



Explicit expressions for the components of tr are given by 



as 
27r2 



(60a) 
(60b) 



while those of x read 

(61a) 

= Jd'^ '<^y^^(lUzUl]kAUlUyU + [UyUlhMu.U - [UyUl]k,Su - SkAUlUyU) (61b) 
[xllhki ■■= fd^z IC^.y ([[/.]«; [f/y(7]f/.]fe, + [UMUy]ulU.]kj - [U^]iilUy]kj - [Uy]alU^]kj 



(61c) 



(61d) 



At this point all the calculational effort pays off: Using the above definitions for x ^-^d f j 
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and remembering that S^^^^[A, j\ J] depends on A via U^^^ only, recasts Eq. H51|l as 
- \{{5A^5Ay)sA[b] (2 (J-S^J^[b,j\j]) (J-SZ[b,j\j]] + ^l-SZ[b,j\j]] e^'^S^^''''''^ 



2" «/"-LJl^ l^^^^-extL '"'J ySby ' ' 7 Sb^Sby 



(62) 



The logarithm comes from an integration over SA with momenta in a finite interval of a;bj values. 
Accordingly the equation has the form of a finite difference equation for Z[J^ , J] with respect to 
Inl/ccbj. 

This gives rise to the final RG equation in differential form, which I write using condensed 
notation with J = {,P , J) and jj = {gjf, jj) as 

-Z[J] = J^J^XuvIyt] + Ju^uIyi] Z[J] (63) 



51nl/xbj ^2"""''^^™^5J^ ' ""■^"^JJJ 

where u and 1; now also stand for q and q in addition to color and transverse coordinates. Eq. (|63|l 
represents a full set of coupled evolution equations for all correlators of J/^^-' also known as the 
Balitsky hierarchy. It is obtained from the above by taking j'^^^ derivatives evaluated at J^'''^ = 
as demonstrated in Eq. 1)52(1 for the dipole cross section. 

Although Eg . 1)63(1 contains all the information we are looking for it is still not written in 
an optimal way. In fact it requires a functional Fourier transform from J to C7 to expose 
its statistical nature as a Fokker-Planck equation. Let me do this, first taking U and U'^ as 
independent variables. 

If one recalls the definition of Z[J\ := (e'^™t['''"^l)b and the initial idea of interpreting (. . .);, 
as a statistical average with associated statistical weight Z according to 



)u= jD[U]{Z[U]...) , (64) 



one may try to rewrite Eq. ((63)1 as an equation for Z. Indeed -remembering that 5gxt[^i J] 
iSgxt[C^: J] is actually given as a functional of U^^^- 



d 



9 In 1 / Xbj 



j D[U]Z[U]e'''^'^^'J''^ = J D[U]Z[U] |ix™[C/] J. + a„[i7] J„| e-^c'^tlt/,-/! . (65) 



Using the fact that if one treats U and W as independent variables one finds, conjugate to 
Eqns. 1(33, that 

± e-^'^'I^^-^l = Jt e-^c'^tlf/../] , _i_ e^Sdu.J] ^ j ^S^du.J] (66) 

Hence, Eq. 1)65(1 can be rewritten as an equation for Z shown in 1)67(1 . Note that by treating U 
and W as independent integration variables [this is necessary for 1)66(1 to hold and to perform 
the translation step] one violates the SU{Nc) condition UW = 1. The variables C/'^-' introduced 
here are more general (invertible) Nc x matrices at this point. This will be rectified shortly 

^^The convention for index contraction again is defined in order to simplify notation. Products arc defined as 
in jto"? := [ [o^'lij with integration over transverse coordinates implied. 
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by using properties of x a and hence the QCD derived nature of evolution for Z. Leaving 
this for later, the evolution equation for Z is given by 

where the J/'-^-' derivatives act on both the kernels x,tT and Z. This equation exhibits exactly 
the same structures as the JKLW equation, an earlier version formulated in terms of the charge 
densities p [7-12], although the variables and details of the kernels are, of course, different. 
Both have the form of a Fokker-Planck (FP) equation, simply due to the fact that the original 
equation as represented in Eq. (|65|l had only terms quadratic and linear in J . 

As with any other FP equation, the fact that the right hand side is a total derivative implies 
that the normalization of Z -its [/, integral- is conserved under Xbj evolution and, if one 
chooses this norm to be 1, allows to interpret Z as a probability distribution at any stage of the 
evolution. 

Clearly this formulation is completely equivalent to the original version involving the generat- 
ing functional Z. Instead of obtaining equations for correlators by functional differentiation, one 
now extracts them, again as in any FP formulation, by multiplying with the desired mononomial 
in U^^"^ followed by integrating the result over U, W . 

Before exploring how to reconcile this Fourier transformation step with the SU{Nc) group 
nature of the physical U fields I use the observation -based on a lengthy but straightforward 
calculation using the explicit form of x a-nd tr- that 

-'?u[«7] =0 (68) 

(for both the q and q components) to rewrite the Fokker-Planck equation in a form that somewhat 
resembles Brownian motion on a curved space 

For consistency one now needs to be able to impose the group constraint on this equation. What 
one needs in fact is, that for physical distributions, which are necessarily of the form 

Zpbys[U, U^] = S{UU^ - l),5(det U - l)Z[U] , (70) 

the constraint factor S{UW — l)S{det U — 1) is left invariant by the evolution operator, i.e. that 
one may commute it though the delta functions. 

In order to get a feeling for what is involved, it is useful to collect some ideas of what to 
expect, if this were actually possible. The main step in rewriting the evolution within the group 
consists of turning ^ and from independent to dependent derivatives via [/C/^ = 1. What 
is to replace them? 

Pondering this question a little (c.f. App. C of [13] for a pedestrian approach), one inevitably 
ends up considering the natural derivatives along SU{Nc), its left and right invariant vector fields. 
The former would be given by-'^^ 

^V?.:=[C/n.4-[-t"t/-].^^^, (71) 



-'^^It is easy to see that this is the correct form as both variants amongst themselves and with each other satisfy 
the SU{Nc) commutation relations at each point. 
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so that as a first step in this direction one would like to reexpress things via [Ut"^]ij-^ and 
[—t^U'^]ij -g^ instead of and jipr^y Indeed, if one now just tries out what this transfor- 
mation will do, one finds remarkable results for the transformed components of x'-^^ 

• in SU (Nc) all four of the transformed components are equal. 

• in SU{Nc) all of them vanish if either A or B is 0. Only octet components survive. All 
matrices involved now are explicitly in the adjoint representation. 

• The unique common form even factorizes naturally into the "square" of a much simpler 
factor. The result, after pulling out a minus sign for later convenience, will be denoted by 
X^y[C/] and reads 

r.l[U]:=-^ • (72) 

The diagonalizing factors of x deserve a name for further reference 

The results of the above transformation of x are so suggestive that it seems obvious what the 
form of the FP operator on physical configurations has to be: — ^iV^^x^^iV^ . The chain of 
argument to prove this is the following. One sets out to show that 

'ifu^^^m:, ^(C/C/t - l)<5(det;7 - 1) = 5{UU^ - l)<5(dctC/- 1) V^X^^V'^^ . (74) 
This amounts to verifying that 

^ -xu.4r u^l ^---^ui^l ^^^^^ = ^?j^xZK, ui-'^ ® . . . ® ui-y (75) 



for arbitrary mononomials U^l ■ . . Uwl, a task that once more is best performed via a 
generating functional. This strategy has the added benefit that it makes it readily apparent how 
the two operators lead to the same hierarchy of equations for the correlators, which provided the 
starting point for all these deliberations. The details of this exercise arc given in App. D of [13]. 

This leads to what is now known as the JIMWLK equation, the RG/FP equation on the 
physical configuration space, where 

ZpY,ys[U, U^] = S{UU'' - l)5(det U - l)Z[U] . (76) 

The constraint factor may now be absorbed into the measure to form a functional Haar measure 
according to 

D[U] ~ D[U]S{UU^ - l)(5(det U - 1) . (77) 

This constitutes the natural measure to use for averaging correlators of the form (j/l,/'' (E) . . . (8) 
ui,J'^)u with the probability distribution Z[U]. The RG/FP equation simplifies to 

^ Z[U]{x^i) = -hvij^r'^K^ mix^i) . (78) 



^*To achieve this technically, one needs to be able to represent all of S/Suf^ components. Accordingly one 

has to replace the index a e {1, . . . , A^^ - 1} by A e {0, . . . , N'^ - 1} and supply t° = l/^/2Nc 1, such that 
tr{t'^t^) = 5'^^/2 for all A, B. Then the transformation is easily accomplished, (see App. C of [13] for details). 
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This formulation is fully equivalent to the initial one, but with a large amount of redundancy 
removed. 

A crucial observation is that the evolution operator on the right hand side of Eq. H78|l 

i/jiMWLK := l^lj^xfy^^u, (79) 

is actually self adjoint and may be called a Fokker-Planck Hamiltonian. It has been shown to 
be positive semidefinite and to posses a (trivial) fixed point at infinite energy with vanishing 
correlation length [13]. 

3.2 The JIMWLK equation and the Balitsky hierarchy 

I will begin this section with a short summary of the results of the previous subsection and then 
turn to their interpretation. The JIMWLK equation describes the evolution of correlators of 
Wilson lines coUinear to the projectile direction which are defined through 

(...)y- = J D[U]...Zy[U] (80) 

where D[U] denotes a functional Haar measure, in keeping with the group valued nature of the 
field variables U. The evolution equation can be compactly written as 

dyZYiU] = -HnMWLK[U]ZY[U] . (81) 

This form highlights its statistical nature as a (functional) diffusion equation. The Fokker-Planck 
Hamiltonian is defined as 

1 -rrb 



i?,iiMWLK ^'^IxZiK (82a) 

^ jdhK.^,y[{l-UlU,){l-UlUy)Y'' (82b) 

{x-z)-{z- y) 
{x - zY{z - y)2 



ab 



and has been shown to be positive semidefinite. In Eqns. (|82|) . both a summation and an 
integration convention is applied to repeated indices and coordinates. 
iV° is a functional derivative operator defined as 

^Vl■.= -\V^e\-^ . (83) 
is the ordinary functional (or variational) derivative w.r.t. the components of the U field: 

= S,kS,iS^^^ (84) 

where S^y := S'^'^^x — y) for compactness. 

This operator in fact corresponds to the left invariant vector field on the group manifold -for 
a physicist, the most familiar interpretation might be that as a vielbein or moving frame. For 
the present purpose, only a few operational facts are needed, which are briefly summarized as 
follows: 
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• Operationally the definition (|83() leads to 

zV?,Uy := - U^t^S^^^ , iVlUl := t^C/t^i^^ . (85a) 

There is, of course, a corresponding definition for the right invariant vector fields : 

i\7lUy ■.=t^U^Si'^ , iVZUl := -t/tr<5(2^ . (85b) 

• Their main properties are the commutation relations (which I display leaving the functional 
nature aside for a second) 

iV''] = if^H^J" [iV", iV''] iP^H^" [iV", iV''] =0 . (86a) 

V and V arc interrelated by 

^'^l = ~[UiT'^y^■. ^yi^-[U^]f^< (86b) 

and "representation conscious:" With the above definitions for the action on U and in 
the q and q representation, it automatically follows from representation theory that acting 
on U or in an arbitrary representation produces analogous formulae with the generators 
appearing on the r.h.s. in that representation. 

• The mathematical concepts of right and left invariant vector fields are used to generate right 
and left translations on the group manifold. For the present purpose this is best expressed 
by writing e"*'^°(*'^°)[/ = t/e"^"*" and e-''^°(*'^°)[/ = e"*'^°*°t/. Functional forms of this 
of course involve an integral in the exponent as in 

e.~^L<i^l)Uy = Uye"^y''^ or e'' ^ ':<^'^'^^^Uy = e-'^*"^ Uy . (87) 

This property has been used in [54] to write an operator that creates soft gluon emission 
in jets. Below, it will be explained how this can be used to summarize efficiently the soft 
gluon content of the projectile wave function entering for example in DIS. 

Using these tools, iJjiMWLK is elegantly written as 

i^jiMWLK = -\^ ^-y [^V^zV^ + + UfiiSJli^l + ^VS^V^)] (88) 

[integration convention for x,z,y\. While the factorized form of Eq. H82() is most useful in a 
derivation of a Langcvin description of the evolution that allows a numerical implementation, 
this second form is more economical in the derivation of evolution equations for given correlators 
that follow as a consequence from Eq. (|81|l . 

This form also clearly distinguishes two types of terms: those that add a new gluon to the 
diagrams at z and those that don't. Diagrammatically this separation corresponds to splitting 
the time ordered diagrams of Sec. 13.11 according to whether the new gluon crosses the a:;~ = 
line and thus interacts with the target or is reabsorbed before or after without interaction and 
thus acts as a virtual correction to the projectile wave function. 

This allows to make contact with Eq. Q used to illustrate the importance of soft gluon 
emission in the virtual photon wavefunction. The contributions shown there directly correspond 
to real diagrams (adding new factors of U into the expressions), the virtual contributions have 
been omitted. It is easy to see now how iterated application of the JIMWLK operator adds new 
gluons, if not on an amplitude level but at least for the DIS cross section as a whole, as shown 
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Figure 11: The JIMWLK operators adds gluons into the photon wavefunction, virtual corrections are 
not shown 



in Fig. ^2 To do so, one needs to write the evolution equation for a general correlator 0[U]. To 
explore the first step shown in FigE] one might choose 0[U] to be 



tr{UiUy)/N, 



(89a) 



which appears in the top left "bare" diagram. For convenience, I also define a shorthand for its 
Y dependent expectation value 

Sy-.xy ■= {Sxy)Y ■ (89b) 

To obtain the evolution equation, one simply multiplies both sides of H81|l with 0[U] and averages 
with (|8U|) . The result then is an evolution equation for {0[U])y 



dY{0[U])Y = -{HnMWLKO[Ul 



(90) 



Because of the nonlinear nature of iJjiMWLK the r.h.s. will in no case be expressible by {0[U])y 
alone, but instead involves other correlators as well. The evolution equation of these new quan- 
tities will then also be needed, and the argument repeats itself, ultimately leading to an infinite 
coupled hierarchy of evolution equations, the Balitsky hierarchy. For Sxy, one immediately 
arrives at 



where 



^xzy {'^^xzy ^xzx ^yzyj 



iVc /y 

(91) 
(92) 

Observe on the r.h.s. the operator [C/2]"^2 tr{t°'Uxt^U^) already anticipated in the corresponding 
discussion of Sec. II. 21 This term corresponds to the real contribution shown in the middle diagram 
of Fig. 111! while the other term is a virtual contribution not explicitly shown. I will return to 
explore this viewpoint in more detail in Sec. [S] 

To compare with the factorizing limit and in particular the BK equation, one applies the 
Fierz identity 



[Uzr'trit'^Uxt'>Ul) = ^{triUxUlMU.UD - ^tr((7. 



to rewrite Eq. as 



xy/Y 



OsNc 

27r2 



Y • 



(93) 



(94) 
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Alternatively in terms of 

NY,a.y - {N^y)Y , N^y = tr(l - UiUy)/N, (95) 

this equation reads 

dY{N^y)Y = jd?z K.^.y[{N^. + N^y - N^y)^ - {N^.N^y)^'^ (96) 

These immediately reduce to a closed equation for Sxy or N^y respectively, if one assumes 
factorization in the nonlinearity on the r.h.s. by replacing 

{SxzSzy)Y Sy;xz SY-zy Or equivalcntly {NxzNzy)Y ^ Ny-xz NY-zy ■ (97) 

The result is known as the BK equation and typically written in terms of NY.xy'- 

dYNY,xy = jd^Z iCxzy [Ny;^z + NY;zy ^ NY,xy - Ny .xzNy .zy] ■ (98) 

To come back to the Balitsky hierarchy, there are several equivalent options on how to write 
the next equation in the hierarchy built on top of Ny-xz if factorization does not hold in 
Eq H96|l . One may consider the evolution equation of {NxzNzy)Y or {SxzSzy)Y as well as 
([1/2] "^2 tT(t°'Uxt^U^))Y or any other combination on the r.h.s. of this "pre-BK equation." In 
any case, iJjiMWLK as given in allows to arrive at these equations most efficiently. 



3.3 General properties of JIMWLK and the fixed point at infinite energy 

There arc a number of quite general properties which have been shown to hold already on the 
JIMWLK level [13]. The most generic of these concerns the interpretation of the JIMWLK 
equation in a probabilistic sense. Due to the fact that the r.h.s. can be viewed as a total 
derivative (c.f. Eq. I|78(l ). or equivalcntly by looking at Eq. H90|l with 0\U] = 1, one finds that 
"probabilities" (the normalization of Zy[U]) are conserved: 

Oy J D[U]Zy[U] = . (99) 

In [13] it has been shown that the Fokker-Planck Hamiltonian in the JIMWLK equation is 
positive semidcfinite. This leads to stable asymptotic behavior: a finite evolution step is given 

by 

Zy[U] = e-"-""^"^'^^-^"^ZYo[U] , (100) 

which, given positivity, means that evolution slowly erases high lying modes. The lowest lying 
mode, the solution with zero eigenvalues, on the other hand, is simply Zy[U] — 1, as the Fokker- 
Planck Hamiltonian Eq. H88() acts with at least one derivative. This solution acts therefore as 
the fixed point at F = cx) and carries the (conserved) normalization of all these "states." This 
state then acts as a fixed point of evolution at infinite energy. In this situation {Zy[U] = 1) all 
there is to define correlators is the Haar measure, which acts locally in each point in transverse 
space. A first implication is that the correlation length in this situation is strictly zero. A second 
one is the fact that, via Eq. (|100() the system moves toward this solution - in RG parlance this 
solution acts as an attractive fixed point. This is fully in line with the physics expectations laid 
out in the introduction: going to higher and higher energy adds more and more gluons of the 
same size, increases density and necessarily leads to a shrinking correlation length. This element 
of the physics picture then is already visible on the abstract level of positivity of the spectrum 
and existence of an attractive fixed point in this evolution equation. 
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3.4 Geometric scaling in BK and JIMWLK - the idea 

At this point one is in a position to make first contact with the scahng idea of Golec-Biernat 
and WiisthofF that provides such an efficient parametrization of the small Xbj HERA data. 
The argument is due to lancu, Itakura, and McLerran [55]. They assume (in the context of 
the fixed coupling BK equation just derived) that Ny xy depends only on the combination 
{x — y)^Qs{Y)^ —: r'^Qs{Y)^. Then y-derivativcs can be traded for derivatives and the l.h.s. 
of the BK equation is easily rewritten as 

dYN{r^Q,(Yf)=r^dr2N{r^QsiYf)dYlnQs{Y)^ r^dr2N{r^Q,{Y)^)2\(Y) , (101) 

where the last equality defines a new quantity, the evolution rate X{Y) = dy ^nQsiY). Taking 
into account the spatial boundary conditions imposed by "color transparency+saturation," i.e. 
the fact that N{r'^ = 0) = and N{r'^ = oo) — 1, allows to take the zcroth moment of the BK 
equation -to integrate with ^ipf~ in order to isolate A on the l.h.s.. The r.h.s then provides an 
integral expression for it in form of its zcroth moment: 

27rA(F) = ^ / ^ (Nin^Ql) + N{v'QI) - N{r'Ql) N{u'QI)N{v'QI)) (102) 

where I have set u = x — z and v = z — y ior compactness. Due to scale invariance of the 
integral, the r.h.s. is indeed independent of Y: the scale common to all N can be chosen at 
will without changing the integral, A is constant. Using the definition of A in Eq. (|101|l . this Y 
independence of A immediately implies that Qs grows like a power in x^y. 

Qs{Y) = e^(^-^")Q,,(ro) = (^)'Q.(^bjo) ■ (103) 

This simple power like growth induced by constant A only holds at fixed coupling. Running 
coupling effects render A Y dependent and thus modify the Y dependence of Qs with the general 
trend to slow down evolution. These effects will be discussed in detail later on. For the moment, 
let me return to the implications of scaling features in small Xhj evolution equations as such: if 
such scaling is to occur, the shape of the solution to the BK equation is strongly constrained 
~ one has to simultaneously satisfy the scaling condition and the evolution equation. These 
combined requirements should be sufficient to completely fix the shape of the scaling solution in 
r, a conclusion that has been confirmed by numerical work again and again. 

The above argument only applies if the scaling situation is ever reached in evolution. Recently, 
Munier and Peschanski [56-58] have argued analytically that this feature should be attractive 
and emerge quite generally as soon as the "color transparency+saturation" condition is met on 
the initial condition. This was done using an approximation that maps the BK equation onto 
equations of the type of the Fisher and Kolmogorov-Petrovsky-Piscounov (KPP) equation [59,60] 
which are known to posses "traveling wave" solutions which in this context appears as scaling. 
While numerical work on the BK equation had demonstrated this earlier on for certain ini- 
tial conditions [61-69] this is a nice theoretical confirmation despite the need for an additional 
approximation. 

This has very important implications: //such a regime is reached, the situation is qualitatively 
different from the situation encountered, say, in the case of DGLAP evolution. There, features 
of the initial conditions and the noncalculable nonperturbative information contained therein 
will be slowly washed out but never completely forgotten. In the scaling regime everything is 
determined by the (perturbativcly calculated) evolution equation - the nonperturbative input of 
the initial condition is completely forgotten in the sense that it is sufficient to state that scaling 
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is reached and that the saturation scale has a given value to fully determine the system at this 
and all further Y . The initial condition will be visible only before one enters the scaling regime 
and will determine the Qs value at which scaling initially emerges, but has no further influence 
on the evolution beyond that point. That this behavior is characteristic for this type of evolution 
equation, also in the more general JIMWLK case, has been explored extensively in numerical 
simulations [35] and will be further discussed in Sec. 

The discussion above applies to the BK equation, but it can be equally applied to the 
JIMWLK equation: With the scaling dimension of the U fields at 1, all the above scaling 
arguments go through, if one assumes all correlators of U fields to carry F-dependence only via 
Qs{Y). In principle this is not impossible and might emerge from evolution in a natural way, but 
it is easy to envision that it is at least violated in a particular physics situation via the initial 
condition to evolution. This should then be seen in experimental data over a certain range of 
Xhj values. 

The experimental implications are also striking. As an example consider the DIS cross section 
of Eq. . If the dipole cross section depends on Y via Qs only, one should expect the DIS cross 
section to scale with Q/QsiY) as discussed in Sec. 18.11 This is the phenomenon originally called 
geometric scaling [34]. 

3.5 From JIMWLK and BK to BFKL, large Nc and density expansions 

Another important result is the statement that the BFKL equation emerges as the small density 
limit of the JIMWLK (and as a corollary, also of the BK) equation. This was first demonstrated 
in [8]. The argument involves an expansion of the dipole correlators {Nxy)Y in powers of the 
gluon field and thus in terms of gluon densities. As a consequence the nonlinear terms in Eq. (|96|l 
become negligible, and it turns out that Fourier-transforming the remainder leads directly to the 
BFKL equation. To state precisely what the expansion parameter for this low density limit may 
be, it is sufficient to consider the ansatz ()33|l and to rephrase the BK equation in terms of Q. 
Since Q is a, two gluon correlator, it is then easy to see that the low density expansion can be 
made precise as an expansion in powers of Q. The name density expansion is justified in the sense 
that Q also fully determines the gluon density as seen from Eq. H43|l . but one has to remember 
that large Q does not literally translate into large gluon densities. 

The evolution I want to consider is in fact the evolution equation for the two point function 
before factorization into the BK form proper, i.e. Eq. (|94|1 . which I write explicitly as 

5y(tr(C/,C/t)) = ^ jd^z iC^.y[-2Ci (tr(C/,C/t)) + {[U,Yhtv{eU^t'ul))) . (104) 

With the two point averages already known from Eq. (|35a|l all that is left to calculate is the three 
point function. This can be calculated along the same lines as the various two point functions. 
One finds 

{[U^f^ 2iY{eU^t^Ul)) = 27VcCfe"3(^<=(ev-.^.+eY-,«.)-%^) _ (^05) 

Instead of an explicit construction let it suffice here that this result satisfies all necessary limiting 
relations. Using, where needed, the Fierz identity H93|l one gets coincidence limits for the fields 
that have to be respected also for their correlators: 

lim \UX'' 2iv{eU^t^Ul) =2Cf iv{U^Ul) (106a) 
lim [[/J"^ 2iv{eU^tW) =tv{ij^Ul) (106b) 
lim [IJ^Y'' 2tY{t''Uj''Ul)=2NcCf . (106c) 
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Obviously Eqns. H35|l and (|105|l arc compatible with these requirements. 

Incidentally this truncates the nonlinear equation for the correlators of [/ to a single, highly 
nonlinear equation for Q. Using the above, Eq. (|104|l reads 

or, equivalently, 

dYgY,^y = ^ Jd^z jc^,y(l - e-^(e-.-+e— e--)) . (107) 

Since one involves only planar diagrams in the construction of the correlators due to locality 
of the weight in y, one might have reached this result also via a different route: One rewrites 
the bracketed expression on the r.h.s of H104|l once more using the Fierz identity H93|l to map it 
onto the form given already in (|96|l . Now one may factorize the nonlinearity in N in the large 
Nc spirit to obtain the BK equation proper. Then one only needs to introduce an exponential 
parametrization of the form (|35a|l . without any reference to the Gaussian weight of Eq. (|33|l at all 
and approximate Cf to arrive again at (|107|l . This highlights the link of the BK equation 

with the large Nc expansion quite clearly and demonstrates that (|107|l is fully equivalent to the 
more common forms of the BK equations shown earlier. 

It is important to note that at this point one has strictly speaking left the realm of the original 
MV model with its assumption of uncorrclated scattering centers in longitudinal direction. Aside 
from the initial condition, t/, as it emerges as a solution to (|107|) . will no longer scale with A^^^ . 
This will be discussed further in Sec. 17.21 

For small Q ()107|l immediately reduces to the coordinate space BFKL equation in the form 

dvQY.xy = ^_^2 j ^xzy{QY.xz + Gv.yz — QY,xy) ■ (108) 

That this equation indeed relates to the better kown momentum space version requires a little 
work, but the correspondence is expected: it has been demonstrated for the forward limit of the 
JIMWLK and BFKL equations already in [8]. Recently, Bartels, Lipatov and Vacca [70] have 
shown explicitly that the lowest order part of (|107|l corresponds to the full non-forward BFKL 
equation (at leading order in ln(l/2;bj)). Since this equation is usually written for a gluon 4-point 
function, this requires an explanation. 

In this low density limit, all that Eq. (|1()7|1 does is to iterate the BFKL variant of what is called 
rcggeized two gluon exchange between projectile and target. This allows to draw both sides more 
symmetrically than in Eqns. (|37|l and (|38|) and to expose the two coordinates (labelled u and 
V below) at which the exchanged pair of gluons couples to the target. This corresponds to the 
step from a the solution of a Bethe Salpeter type equation to the corresponding Greens function. 
This latter object will be denoted Gv; xy.uv: with Fourier-transform S^^''{q — q')^Y{k,k';q'). 
Diagrammatically, this is written as 




GY;xy,uv=^^ ^Y{k,k';q')= . (109) 
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[The definition in terms of underlying gluon correlators now contains local subtractions of the 
type shown in H38() both for top and bottom coordinates.] Explicitly, with this assignment of 
variables 

6^^\q- q')^Y{k,k';q') ^ J d?xSyd?ud?v Gy, a^y^uv e''=^+*('J-'=)2'-'*='"-^('j'-'=')'' (110) 

and to make contact with the standard momentum space formulation reduces to a simple exercise 
in Fourier transformation. 

The connection between the result of these manipulations and the conformal formulation in 
coordinate space is well documented in the literature^^ so that taking this route might be easier 
to the non-expert than the direct calculation of [70]. Only by using conformal symmetry tools 
is it possible to write complete solutions to the nonforward BFKL equation - but this is outside 
the scope of this review. 

While the authors in [70] arrive at the BFKL limit in the usual way, by simply dropping 
the nonlinearity in Eq. and hence refer to N as the object of their evolution equation^" 

their calculations make it clear that they discuss a two gluon object coupling to the photon 
wavefunction or "impact factor." This object has to be identified with G as done above. 

The form of the nonlinearity in Eq. (|107|l is intriguing in its exponential nature. This directly 
corresponds to the eikonal nature of the multiple scattering induced by the presence of complete 
Wilson lines in the definition of correlators Sy-xy or Ny-xy This becomes important as soon 
as one attempts to go beyond leading order in density, as becomes necessary if one tries to step 
beyond the large Nc limit. For example [70] aims at providing a 1/Nc correction to the BK 
equation. They find a contribution which contains the square of ^f-{Qy^xz + QY,yz ~ QY,xy^, 
the factor in the exponent of the (iVc-leading) expression H107|) . If the eikonalization feature of 
JIMWLK and BK equations is taken seriously, this can not be done in a way consistent with the 
density expansion without also including higher order contributions from the nonlinearity on the 
r.h.s. of H107|) . At least for large nuclei this becomes an important consideration. 



3.6 Growth in BFKL and BK, the Mueller- Triantafyllopoulos description 

To identify BFKL as the low density limit of JIMWLK and BK implies that linear BFKL 
evolution is present in the more general framework and should manifest itself in those areas of 
phase space in which densities are not yet large. BFKL is known for fast and unlimited growth 
of its solutions, the C/ or $ above. The role of the nonlinear effects is clearly to tame this growth 
and to ensure that the objects of evolution, Sy-xy or Ny-xy for example, with their interpretation 
as correlators of eikonal factors Ux'' stay bounded between and 1.^^ 

Mueller and Triantafyllopoulos [73, 74] have provided a very elucidating perspective on this 
matter that employs BFKL as the driving force behind all growth at small Xbj which turns out 
to be surprisingly successful. The ideas entering this model naturally incorporate the notion of 
a scaling window [55] - a region in phase space at distances shorter than 1/Qs, i-c. away from 
the high density limit, at which the presence of Qs exerts its infiuence. 

As indicated, one cornerstone of these ideas is the strong growth and diffusive character of 
BFKL evolution which can be most easily understood in the forward case (i.e. the above equation 

^^See for example [71], for a pedagogical presentation. 

■^'^Note that both A'' and Q vanish in the local limit, a fact that is called the Moebius representation in [70]. 

■^^In the case of the BK equation, which can be written without any reference to this interpretation via eikonal 
factors, this boundedness has to be a property of the evolution equation itself. [72] has given an argument in the 
context of the Banfi-Marchesini-Smye (BMS) equation, why this is the case. Their argument directly translates 
to the BK case due to a structural analogy of the BMS and BK equations that will be discussed in more detail 
in Sec. |5] 
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at g =: 0). At g = the above equations can be diagonalized by power functions 

e('^'")(fc) = 27rV2(fc2)-i-''^e-™"^ ; G M, n £ Z (111) 
in momentum- as well as in coordinate space. The eigenvalues are given by 



X{v,n) := 



(112) 

with ■0 defined as the logarithmic derivative of the F-function {iIj{x) = r'(a;)/r(a;)) and x{^) 
a shorthand for 0). The solution then reads 



+00 /•+00 J 

$Y(fc,fc')= y / — e('^'"Hfe) e^'^('''") e('''")*(fc') • (113) 



n— — oo 



Of particular interest is the large Y asymptotics of this equation. There the contributions will 
be dominated by the modes with n = 0, shown here in the coordinate space version to underline 
the analogy [28]: 

= yfr^^-'e^iiu - vf)-^'^'^ = j ^^yxi^)~^^~o)^P~p') (114) 

where p \vl(x — y)^/^^ and p' := ln(tt — t>)^/i^. The integration contour in the v (or Mellin-) 
integral is to be taken parallel to the imaginary axis, to the right of all singularities; in the last 
version I have substituted v := iv with an appropriate change in definition of x- 

Further insight emerges by taking the saddle point approximation to in Eq. (|113|l or (|114|1 . 
This selects z^o = as the extremal value in the v integration. The result displays a factorizing 
structure, 

gi-xli^o) / '^'^ ^yj;y"(i;o)(£'-i>o)" + (i>-£>o)(p-p') ^ gyy(->o) ^ p- x^T^^pjy fii5) 

where the first factor, e^^'^°' = [n^J ' describes the strong, power like growth in l/xbj 
alluded to already in the introduction. The second factor, which originates from the Gaussian 
integral in i), describes diffusion in the transverse plane of dipole sizes: 

i^Yip-p') ■■= , ^. . e-^^^ (116) 

is immediately recognizable as the solution to the diffusion equation 

dy^Yip) = ^^a>y(p) (117) 

with initial condition ipoip) = 6(p). The same picture holds in momentum space as well. This is 
the well known IR problem of the BFKL equation at one loop accuracy: even an initially short 
range configuration spreads to nonperturbatively large distances. Since the main focus here is in 
nonlinear effects I will refrain from a discussion of two loop corrections to BFKL. 



The result does not depend on ^. 
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This diffusion and growth now has to be contrasted with what happens in JIMWLK and BK 
evolution. The nonhncarity tames the growth and, at the same time, prevents diffusion into the 
infrared. At small Q, though, evolution is correctly described by BFKL dynamics. 

Mueller and Triantafyllopoulos therefore have suggested to simulate the presence of the non- 
linear term by an absorptive barrier in the above diffusion process that is to be self consistently 
adjusted such that the solution to the so modified BFKL equation never grows large. The result 
can be interpreted as a change of the dominant contributions to the integral in, say H114() . in 
such a way that one now follows the saddle point but allows diffusion only from small dipoles at 
size ranges where G is not yet large. Creation of new objects in the already saturated domain is 
eliminated via absorption across the boundary. 

In order to explore the first part of this argument and to understand how to choose the 
absorptive boundary, one can look for directions in Y — Q phase space in which BFKL evolu- 
tion remains bounded by simultaneously looking for a saddle point and Y dependent starting 
sizes p'(Y), for which the exponent in Eq. (|114(l vanishes and thus G remains approximately 
constant [55,73]. The conditions to satisfy in this case are 

yx'ii'c) + {p- p'{Y)) = saddle point (118a) 

Yxii>c) + (1 - i^c){p - p'{Y)) = (almost) const. G (118b) 



with solutions 

X{Vc) 



(119a) 



p-p'(y) . (119b) 

1 - I/c 

(The saddle point value Vc remains implicit, but changes slightly from the BFKL value vq.) 
Translating the latter back into distances p =: ln(r^/x^) and p'{Y) =: lni?Q(y)/x^ one finds 

RliY) =r^e-^'^ . (119c) 
Using this as a saddle point for the integral in G one finds 



^g>'|x"(i>c)(S-Se)^ + (5-i>e)ln(^) 

Ro{Y) J J 2txH 

.2 ^- 



A list of observations apply [73] : 



(120) 



• Within the diffusion radius the dominant factor is (r^/i?o(i^))''^ ^"^ ■ The solution has 
scaling form, just as argued above for the solutions of the JIMWLK and BK equations after 
some initial settling down. The scale variable is i?o which has exponential Y dependence, 
again as required in the BK case. Moreover the coefficient is known from Eq. I|119(l . 



Apart from a mild correction through the square root factor, i?o(i^) describes lines of 
constant G ■ G evolved in this manner will not grow exponentially. This is because emission 
from large objects has been prohibited in sharp contrast to the diffusion of the BFKL 
equation. 
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Of course the argument can be equally phrased in momentum space which essentially replaces 
r ^ 1/Q and Rq 1/Qq. 



One may remove the nonscaling square root factor by replacing condition (|118b|l . with the 
requirement of truly constant Q and work with the resulting definition of Qq Qg without 
introducing absorptive boundaries. However, evolution speed A = cJylnQs, will necessarily 
change if one really implements the effect of the nonlinearity in form of the absorptive boundary. 
Mueller and Triantafyllopoulos [73, 74] continue to formalize their argument by also including 
running coupling effects. In their treatment, the scale for the coupling is set by the size of the 
emitting object, the parent "dipole." These refinements do not change the generic conclusions 
given above but they ensure scaling and provide quantitative results for A and Qs- Due to the 
approximations made, their results are valid asymptotically at large Y and contain an unknown 
constant (called Y below) which is related to the (uncontrolled) small Y details. I write their 
result for the evolution as in [35]: 



A 9y In g,(y) 



.90 



0.47 



VyTy (r + r)5/6 



higher inverse powers 



This can be integrated to give Qs at the price of an additional integration constant c'P'^ 



Qs{Y) 



A, 



QCD 



■ cxp[2 • o.9(r + y)5 - 6 • o.47(r + y)s] 



(121) 



(122) 



In particular the coefficient of the second term is affected by the presence of the absorptive 
boundary. Its value will turn out to be important in the comparison with numerical simulations 
shown in Fig. 1231 

Clearly running coupling effects do slow down evolution quite considerably, the leading contri- 
bution goes from an exponent linear in y to a square root in the exponent but they do not destroy 
scaling. The functional form of the diffusive factor in Eq. Ijl20() is modified through the correct 

determination of lines of constant Q but functional form of the leading {fp / Q siXY) '^'^^ is 
imaltered (although Vc changes slightly from the solution to (|119|l ) [73]: 



^Y,Q.,Q,{Y) 



C 



-(i-!>c) r 



In- 



(123) 



Dc (or in pure BFKL) is usually called the anomalous dimension governing the leading large 
k behavior of the unintegrated gluon distribution through (|43|) . This result is expected to hold 
in what is called the scaling window which will be discussed next. 



3.7 Scaling above Qs, the existence of a scaling window 

The scaling form of Eq. H123|l is valid at small Q which immediately implies distances shorter than 
l/Qs- Consequently, one has to ask how far away from this can be trusted. This, in fact, is 
important not only for the approximate solution discussed above but also as a limitation on the 
range of validity of the evolution equations, be it JIMWLK or BK. That this range is limited 
as a matter of principle should be clear from the derivation of these equations, which selects 
contributions of the form ofg ln(l/a;bj) and rejects contributions enhanced by other logarithms, 
such as asln{Q'^ / fY) as subleading. If one looks at the phase space plot of Fig. 01 one would 
naively expect ln(l/a;bj) corrections to be most important above the diagonal, while ln((5^//i^) 

^''Alternatively, one could have directly used the expressions for Qs in [73, 74] which contain a free parameter 
with the same role. 
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corrections should be most important below the diagonal, with a region in between in which 
double logarithmic corrections of the form as ln(l/xbj) ln(Q^//i^) are most important. But where 
can one safely neglect these other logarithms? 

Even before Mueller and Triantafyllopoulos constructed their scaling solutions, [55] have in 
fact argued that the range of validity of the evolution equations and also of scaling expressions of 
the type (|123|l is surprisingly large. It is expected to extend over a large scaling window defined 
by the condition 

Ql<Q'< Qt/A|cD (124) 

or, more stringently 

1 < HQVQ^s) « HQI/A^cd) ■ (125) 

The original argument is based on the diffusion radius of BFKL if one follows evolution along 
lines of (almost) constant Q as in Ea. p2()(l . This sets a natural range for the validity of scaling 
formulae like H123|l . 

Besides the theoretical considerations above, the existence of a scaling window is also impor- 
tant phenomenologically. The key example is the G-B+W idea to fit the HERA data with a 
scaling ansatz. In the HERA range Qs comes out to lie in between 1 — 2GeV^, while the momenta 
entering the fit range extend up to of several hundred GeV^. Thus it is natural to expect 
these data to exhibit scaling features only if a scaling window exists that reaches far beyond Qs- 
In a sense, the success of the G-B+W fit directly poses the question for the existence of such a 
window whose theoretical origin was only understood much later in [55]. 



4 Fokker-Planck and Langevin 

At this point all of the regions marked in the phase space plot of Fig.^have been fully explained 
in terms of a QCD evolution equation although for some of the steps certain approximations have 
been used. To go beyond this stage numerical work is needed. To allow for an efficient numerical 
treatment of the JIMWLK equation a reformulation is required. This takes the form of a trans- 
lation of the Fokker-Planck equation to a Langevin description. This is done in Secs. l4.ll and l4.2l 
Interestingly enough, in this process one finds theoretical tools that also help in interpreting the 
gluon radiation process and directly lead to the diagrammatic interpretation of Eq. ((T)). This 
step is performed in Sec. 14.31 while the interpretation will become clear only in Sec. by way of 
a surprisingly close analogy with the theory of QCD jets. 



4.1 From Fokker-Planck to Langevin: Illustrating the idea 

Although this topic is covered in many textbooks, the presentation is often somewhat old- 
fashioned and the intimate connection with path-integrals is not always mentioned. For this 
reason I demonstrate the relationship of Fokker-Planck and Langevin formulations with the aid 
of a simple toy example, a particle diffusion problem, that, for ease of comparison, shares some 
of the features of our problem. The toy model equation is 

dyPYix) = -Hfp Py{x) (126) 

with a Fokker-Planck Hamiltonian 

HYP^\^^^X^.Ax)i^''^■ (127) 
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The particle coordinates x (in D dimensional space) correspond to the field variables of the 
original problem and the probability distribution Py to Zy ■ It serves to define correlation 
function of some operator 0{x) via {0{x))y = J d^x 0{x)Py{x). 

Eq. H126(l admits a path-integral solution which is, as usual, built up from infinitesimal steps 
in Y: this is a solution P{Y' , x'; Y, x) with initial condition 

P(Y,x';Y,x)^d'^^'>{x' -x) (128) 

that is composed of many, infinitesimally small steps at times Yk := Y + ke. In this solution, the 
product rule 

P{Y', x'; Y,x) = J d^Xn-i . . . d^xi P{Y' , x'; Yn-i,x)P{Yn-i,x„-i; Y„-2, Xn-2) ■ ■ ■ P{Yi,xi\Y, x) 

(129) 

expresses the finite Y' — Y solution P(Y' ,x' \Y,x) in terms of infinitesimal steps P(Yi,Xi\Yi--i_,Xi^\). 
Arbitrary initial conditions Pyo{x) are then recovered via Py{x) = / d^yP{Y,x;Yo,y)PYo{y)- 

The derivation of the infinitesimal step from i^-i to Yi = Yi-i + e is textbook material. It 
reads (the index i on the coordinate also refers to the time step, not to the vector component) 

F(y„a;,;y,_i,a;,_i) -iV fd^p, e-^(^P>?x°''(--i)+*P'-(^"^--(--i))) + 0{^) (130) 



where (J^{x) := ■^id^^Xani^) and A'^ is a coordinate independent normalization factor. 

From here, the step to the Langevin formulation is trivial: as this expression is quadratic in 
the momenta one can trivially rewrite it with the help of an auxiliary variable [i again the 
step label] as 

P{Y,,x,;Y,^uX,^,)=N jd^'i, yditW^e-^^''^"^^')^- {x, - [x,^^ + e{i, + n{x,^{}]) . 

(131) 

The 5^{. . .) arises from the momentum integration and determines Xi in terms of a;,;_i and the 
correlated noise ^^.^^ The equation for a;,, 

X, = x,^i + e(e, + cr(a;,_i)) , (132) 

is called the Langevin equation. It contains both a deterministic term [the a(xi-\) term] and 
a stochastic term [the ^ term]. To fully define the problem without any reference to its path- 
integral nature, one then needs to state the Gaussian nature of the (correlated) noise separately. 
This is the formulation often found in textbooks. The path-integral version would appear to be 
much more elucidating and up to date.^^ 

If, as in the JIMWLK case, x factorizes as 

Xixv{x) = £^a{x)£ua{x) , (133) 

a simple redefinition of the noise variable leads to a version of the Langevin formulation with a 
Gaussian white noise: 

PiY,,xf,Y,^i,x,^i) = N jd^i, e -^^- S'^ix, - [a;,_i + e{£{x,^i)i, + a{x,^^)]) . (134) 

Re-expressing the delta function via a momentum integral and performing the Gaussian integral over 
immediately recovers Eq. ImH 

^^One of the main sources of confusion with stochastic differential equations, the choice of discretization [the 
details of the Xi and Xi—i dependence] and its impact on the form of the equation itself, finds its natural 
explanation there. It corresponds to the choice of discretization in the path integral solution Eq. I13UI . The 
choice here is called the Stratonovich form. One might have as well have given the Ito form by choosing a 
different discretization for Eq. ll3Ui or yet another version with the same physics content. 
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The correlation is now absorbed into the Langevin equation through the appearance of S'^^ which 
reads 

X., = a;,_i + e(£(a;,,_i)|j + a{xi^i)) . (135) 

Wherever this form is available, it will be the most efficient version to use in a numerical simu- 
lation, as the associated white noise can be generated more efficiently than the correlated noise 
of the more general case. 

What is left is to sketch a numerical procedure to implement the solution given in l|129(l . 
H13Q() using the expressions Eq. H131(l or Eq. (|134|l . As a first step one replaces the average with 
the probability distribution Py{x) by an ensemble average 

{0{x))y= [d^xO{x) Py{x)^^ ^ 0{x) (136) 

c/eE[Pv] 

where N is the size of the (large) ensemble which is taken according to the probability distribution 
Py{x). This one can easily do for the initial condition at Yq. The Langevin equation then allows 
to propagate each ensemble member in Y by e, thereby providing a new ensemble at Yi = Y + e. 
Iteration then creates a chain of ensembles for a discrete set of Yi that provide an approximation 
to a set of Py^ (a;) in that they allow one to measure any correlator {0{x))y according to Eq. I|136(l . 

The above clearly shows that a "Langevin system" is a way to write a path-integral solution to 
a differential equation that is particularly easy to implement numerically. Many, often confusing 
issues, such as the time step discretization issues, find their natural resolution in this path integral 
setting. 



4.2 A Langevin formulation for JIMWLK evolution 

Translating Fokker-Planck equations to Langevin equations for numerical work is not only useful 
in the case of a particle system as portrayed in the previous subsection. Also in the case of 
the JIMWLK equation this is the only promising angle of attack. Besides the evolution "time"- 
discretization needed in any example, it is also necessary to discretize transverse space to render 
the number degrees of freedom finite. Similar methods have been discussed in lattice gauge 
theory under the general heading stochastic quantization. 

In the present context, this approach has already been suggested in [13] and a derivation has 
been given in [75] to which I refer the reader for details. To affect the translation in analogy to 
the example above, one may follow [13] and write the JIMWLK equation as 



dYZYiU] = -VI - {\^'yXty) Zy[U] (137) 



and define 

K := = -^{~ Jd'zj^^trimU.)) (138) 

(the indicating adjoint matrices and traces). 

One then abandons the description in terms of weight functionals Zy in favor of one in terms 
of ensembles of fields which are governed by the corresponding Langevin equations. 

Explicitly, to calculate any observable 0[U] of the fields U one writes 

{0[U])y^ Jd[U]0[U]Zy[U]^^ J2 0[U] (139) 
ueE[ZY] 

■^^The dimensions D and n of configuration space and noise now need not be equal as is illustrated by the 
JIMWLK example. Again, write the <5-function as a momentum integral and carry out the Gaussian integral in 
^ to recover Eq. ITi^ from Eq. tT^ . 
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where, separately at each Y, the sum is over an ensemble ^[Zy] of N configurations U whose 
members were created randomly according to the distribution Zy- Clearly, for N —* oo, the 
ensemble and Zy contain the same information. 

The Langcvin equation then schematically^^ reads 



dy [C/.]y- = [ jSy [SZlUOk + (140) 



where 



£t. = ( ^ 



(141) 



TT^ / (x — y 



is the "square root" of x, xSy = ^x'z^z'y^ ^^"^ ^ ^'^'^ independent Gaussian random variables with 
correlators determined according to 

(...).= y'i?K](...)e-^« . (142) 

This clearly represents an equation for an infinitesimal change of an element of SU{Nc). In fact. 



the components of u;^ 



fc 



can be directly interpreted as the "angles" 



parametrizing a local gauge transformation in transverse space. 

It is of particular importance to note that the possibility to formulate the stochastic term via 
a completely decorrelated Gaussian noise that is to say with (C'Cy^ ) = S"'' S'^ S^^y , reduces 
the numerical cost for a simulation like this considerably. 

All that is left to do conceptually is the explicit discretization and the creation of an initial 
condition in form of an ensemble of U fields at Yq. 

A reasonable choice for the initial condition is an ensemble that leads to a shape for {U^Uy) 
that is compatible with the generic physics requirements as shown in Fig. 13 for instance with 

(UlUy) (X exp (^^-^) . (143) 



4i?2 

The scale R here is close to the initial saturation scale Rg (depending on the precise definition 
of R,). 

A possible choice for the discretizations [35] is as follows: one discretizes the Langcvin equa- 
tion in Eq. p4()|l using a regular square lattice of iV^ sites (volume (Na)"^, where a is the lattice 
spacing), with periodic boundary conditions. In dimcnsionless units (a = 1) and defining a 
rescaled, discrete, evolution time 

^Y (144) 



7r2 



the Langcvin equation becomes 



Ux-s+Ss = Ux-s exp[ii"w°.^] (145) 



where 



/^^i^X^ - y)[l - Ul.Uy^sl'^'C;' - SsY,S{x y)lil■{^Wl,Uy^s]■ (146) 
y y 



^'^This is a continuous time version of the equation that strictly speaking is not unique. The path-integral 
derivation shown in [75] and maybe more transparently in the appendix of [35] makes it clear that we are to take 
a "retarded" prescription here in which the derivative on the l.h.s is taken as a finite difference and the fields on 
the r.h.s. arc determined at the previous time step. 
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Here K'{r) = ^, S{r) = ^ and C is a Gaussian noise with = S'^^'S^^^y. Ki{r) and 

S{r) are taken to be = if r = 0. This is justified by the fact that the apparent singularities at 
these points are completely canceled in the JIMWLK Hamiltonian these Langevin equations are 
equivalent to. 

As a consequence of this additional discretization, one needs to carefully check that the results 
are free of any artifacts in the sense that the numerical results can be reliably extrapolated to 
infinite volume, that the spatial continuum limit can be extracted and that "time" discretization 
errors are negligible. Most of these are rather technical issues and I refer to [35] for details. A 
careful study of the field theoretical continuum limit, however, has led to a clearer understand- 
ing of the nature and importance of running coupling corrections for the sclfconsistency of the 
method. This will be reviewed carefully in Sec.lHI 



4.3 Shower operators: the operators that generate soft gluon clouds 

The translation to Langevin form was motivated above with numerical implementations. There 
is in addition a conceptual side effect of this reformulation, whose importance is not obvious at 
this point of the discussion. I will nevertheless start to prepare for the discussion in Sec. 15.81 
by writing an apparently formal solution to the Langevin equation that contains a version of an 
operator dubbed a "shower operator" in [54]. This operator can be interpreted to build up the 
ordered soft emission part of the projectile wavefunction, the soft gluon cloud, although this is not 
fully apparent from the present context. This interpretation will become obvious by comparison 
with the treatment of soft emission in jets in Sec. [3 Let me now come to the formal solution 
itself. To this end I rewrite the r.h.s. of Eq. (|140|l in terms of a "single emission operator" (also 
this terminology will become self explaining in the discussion of amplitudes in Sec. [51): 



dy Uy,3 



i I d ud V 



This can now be integrated using K-ordering: 



(147) 



Uy., 



■"y cxp ■ 



dY d^ud^v 



]Uy„. 



(148) 



The operator 



exp 



{'f 

JYo 



dY d^ud^v 



(149) 



is the shower operator responsible for soft gluon emission. Please note the analogy with Eq. H87(l : 
the only special feature is the complicated nature of what plays the role of the angle variable 
uj. This now contains U at earlier Y and in a nonlocal manner. This is a nice illustration of 
how the JIMWLK equation implements the resummation of the second nonlinearity: the change 
in U turns nonlinear in U itself. One might be tempted to call that "second eikonalization:" 
eikonal factors iterativcly enter the exponents themselves. That these operators indeed have 
strong analogies with path ordered exponentials is obvious from the appearance of Y ordering. 
This leads to useful properties as will be seen in identities like Eq. H2U2|) . 

The relevance of this type of operator lies in the fact that it not only appears in the Langevin 
context, but also in the construction of amplitudes with arbitrary numbers of soft giuons. n-th 
order truncations simply correspond to n-gluon amplitudes. An example at small x is provided 
by the virtual photon wave function entering the DIS cross sections in Eq. Q), now taken to 
all orders in soft gluon emission. In this sense the operator in H149|l creates the "Cold Gluon 
Cloud" [Alex Kovner's favorite rendering of CGC] and clearly summarizes the same physics 
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as the JIMWLK equation. I will now present an alternative derivation of nonlinear evolution 
equations of the JIMWLK type based on soft gluon amplitudes and gluon shower operators. To 
show the generality of the concepts I will do this for a new physics example, a specific type of 
jet observable. 

5 JIMWLK and soft gluon emission: comparing with and 
learning from jet physics 

At first sight it would appear to be fairly far fetched to expect saturation type physics of any kind 
in a jet measurement, in particular if the jets are created in "empty space," for example by an 
e+e~ collision. There one is typically looking at observables that have long been described purely 
perturbatively as multiple emission processes without any reference to density at all. Although 
it is conceivable that in certain corners of phase space -for example by looking at in some sense 
"strongly focussed" jet events- that these particles start to overlap in a way comparable to 
the small Xbj situation, in the sense that one would have to take density effects into account, 
most would have relegated such phenomena into the realm of extremely rare events with little 
experimental relevance. 

This has changed with the identification of so called non-global jet observables by Dasgupta 
and Salam [76, 77] already for e~^e~ collisions. Furthermore, it turns out that for these same 
observables the importance of the nonlinear effects is increased if the jets originate in or travel 
though a dense medium created by other means - for example in a heavy ion collision. It is to 
be expected that many jet observables accessible in the heavy ion context will acquire non-global 
characteristics and become sensitive to medium and density effects. This type of consideration 
will become particularly important in the heavy ion program at the LHC. I will turn to a brief 
discussion of this situation in Sec. 15.81 but will present the underlying theory in the simpler 
e+e~ case. Already in this situation one finds evolution equations with similar structure and 
interpretation of the nonlinear effects as in the small (c.f. Eq. lfT5^ and the BK 

equation (|161|l ). 

To formulate an archetypical observable of this type, consider a e+e~ collision at total in- 
variant energy and look at a 2 jet event (considering 3 and more jet events as a perturbative 
correction for the time being) as depicted in Fig. Typical soft- and coUinear-safe jet observ- 




Figure 12: Two jet event in e'^e~ with prescribed geometry and soft emission into the region away 
from the jets. Cin and Cout denote the regions inside and outside the jets which for intuitiveness are 
characterized by an opening angle 9in. 

ables are for example broadening B or thrust T, both defined as sums over kinematic variables 
of measured hadrons. To be explicit, B ^Yl,h ^'^^o^^ 1 — T = Yl,h ^''^'"'''n — — where {pt)h and rjh 
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are (three dim.) hadron transverse momentum and rapidity w.r.t. the jet axis. Setting V = B 
or 1 — T the corresponding jet shape distributions are then defined as 

where dan{Q) is the n-hadron distribution, (tt{Q) the total cross section and Vh is a single 
hadrons contribution to the observable. The most important perturbative contributions to such 
variables are enhanced by soft and collinear phase space logarithms, a single emission carries at 
most both of them. 

Thrust and broadening belong to a large class of observables in which these enhanced con- 
tributions exponentiate [78] to give a systematic perturbative expansion for the exponent of the 
form^^ 

oo 

lnE(g,y) = ^a3(g2)"(A„ln"+V + B„ln"T/ + ...) . (151) 

The An series is referred to as double logarithmic and i?„ as single logarithmic, the remaining 
contributions are not logarithmically enhanced. The variables considered up to this point are 
referred to as global: all final state hadrons are counted in the sums "Ylih^hi no matter where 
they end up in relation to the physical jet. Non-global observables are more restrictive, they 
measure for instance contribution only from hadrons that end up in the region away from the 
physical jet. To prepare for such a measurement one distinguishes the region containing the 
jets (Cin) from the region safely away from the jets (Cout) as indicated in Fig.^l To provide a 
prototypical example, consider the "inter jet energy fiow" defined by 




(152) 



Here Vh = ojt is the energy of the measured hadron and the sum is restricted to Cout- No event 
with more than a total energy i?out outside the jets is counted. In a perturbative calculation, 
the phase space in energy extends up to ~ (5/2 which will give rise to large logarithms in 
{Q /2)/ Eout ~ E/Eo^t, but there will be no collinear contribution since the measured particles 
are emitted at large angles (from the jet as a source). This implies that in the representation (|151|l 
only the single logarithmic terms contribute. While this is an important difference with the earlier 
global examples, it is not immediately evident, why this observable should bear any similarities 
with the dipole cross sections of the CGC. The careful analysis of Dasgupta and Salam [76,77] 
has uncovered that in distinction to the global case these single logarithmic contributions are 
sensitive to the internal structure of the jet. To the required accuracy the measured soft particles 
can not be treated as if emitted by the initial hard qq pair. Through emission from secondary 
"hard" particles, the observable becomes sensitive to the internal structure of the jet. To retain 
this information one must keep track of the nonlinear multiple emission structure inside the jets. 
Nevertheless, the connection with BK and JIMWLK evolution was not apparent until a beautiful 
paper by Banfi, Marchesini, and Smye (BMS) [72] reformulated the numerical approach of [76,77] 
(in the large Nc limit) for exactly the observable described above in the form of an evolution 
equation. The equation itself is nonlinear in the part describing contributions arising inside the 
jet regions and linear outside. The contributions inside, quite surprisingly for most experts, have 
a striking analogy with the BK equation. For small angle emission in a certain frame, the analogy 

^*This result is nontrivial: For example it states the absence of contributions of the form a" In™ V with 
m = n + 2 . . .2n in the exponent. Only one of the "exponentiated" gluons carries both logarithms. 
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becomes even more striking: the contributions of BMS inside the jet regions become identical to 
the BK equation. That such an analogy is restricted to the in-region, is clear: in the derivation 
of the BK equation one assumes, in fact requires, the presence of non- linear effects everywhere 
and any geometric exclusivcncss of the type imposed by jets was simply never considered. What 
appeared to be more surprising is the relation of jet physics to BFKL dynamics, which dominates 
the BK equation in the small density limit. This link has been established in [79] and indeed 
comes into play in the limit of small angle emission in a judiciously chosen frame. 

Subsequently in [54], I have shown that the analogy goes further: there exists a functional 
evolution equation that generalizes the BMS equation to finite Nc which parallels the JIMWLK 
equation. This should be not too surprising, as BK emerges from JIMWLK by means of a 
factorization assumption at least strongly linked with the large Nc limit. There arc useful lessons 
to be learned for small Xbj physics from both the structure and the derivation of this new jet 
evolution equation. The structure teaches us about possible consequences of exclusive constraints 
on observables. It shows, for instance, how non-emission requirements lead to Sudakov like 
resummations in the restricted areas -there only the first of the rcsummations into eikonal factors 
is at work-, whereas the unrestricted areas remain subject to nonlinear evolution and experience 
both types of nonlinearities. The derivation was built on earlier results on amplitudes [80, 81] 
and thus gives direct access to them. This should provide greater technical flexibility for any 
generalizations that might turn out to be desirable also on the small Xbj side. Not least has the 
presentation in the introduction been heavily influenced by the results derived in this context. 



5.1 The analogy of BMS and BK equations 

To present the analogy, let me give a short account of the BMS equation before I compare with the 
BK equation. I will closely follow the exposition in [72] . The observable considered by BMS has 
been given already in H152|l . The process is dominated by iterative soft emission from eikonalized 
gluons that are not deflected from their original trajectory. The natural kinematic variables are 
thus the directions of the emitters {a,b,k,p,p in what follows) and the corresponding energies. 
Phase space integrals will be written in terms of solid angles and energies. 

A Mcllin transformation represents the cross section in terms of transition probabilities 

Gab{E, i?out) 

/dv e'^^°"^ 
— - — : Gpp{E,iy ) ~ Gpp(£^, i?out) , (153) 
Zmi' 

where p and p are the directions of the original q and q respectively and a saddle point approxi- 
mation has been taken. 

Gab{E, Eont)i or G ab{E) for short where the explicit i?out dependence is not needed, satisfies 
the BMS equation 

EdE Gab{E) - j ^ asWabik) [u{k) Gak{E) ■ Gkb{E) - Gab{E)\ , (154) 
where (Xg := asNc/ir and 

(h\ = (Pa-Pfc) ^ 1 - cos Ogb 

'"'""^ ' [pak)[kpb) (l-COS^afc)(l-COS0fc,) ^ 

is the soft emission kernel. Conventions are adapted to lightlike momenta and the energies, 
ijjp,ujq have been factored out in {pq) compared to the 4-vector product: p.q = LOpU!q{pq) = 

UJpLOqil - COS 6 pq). 
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The jet geometry is encoded in the definition of u{k), which emerges from a Melhn factoriza- 
tion of the energy constraint in Eq. H152() : 



/ \ f dv p'^^out 

e^Eout- '^'j = / 2 v n"(^«) ' "(^■) ^ + e"''"0out(fc) , (156) 



with the and Gout having support inside and outside the jet regions respectively. 

The strategy employed by BMS in the derivation of Eq. (| 154(1 is surprisingly straightforward. 
The evolution equation was derived from the knowledge of the structure of the real emission part 
alone. Virtual corrections were reintroduced using the requirement of real-virtual cancellation in 
the infrared. 

The real emission part of G needed for this argument can be written as 



G(r')(£;,Kut)=i+E/ n 

n=V^ i=l ^ 



as — ^u{h) eiE-Lu,)}Wniph ■ ■ ■ KP) (157) 



where the phase space of soft gluons is cut by E. Wn is the large Nc factorized version of the 
transition probability of a hard qq color singlet into ((7g)hard5"ofti 

Wn{pk,...Kp) = ^} . (158) 

The line of argument then starts by taking a logarithmic derivative EOe of Eq. p57|l which 
uniquely fixes the quadratic term in Eq. (|154|l . This term describes the creation of a real gluon 
in the final state with its energy at the phase space boundary E. This term is affected by 
the jet-geometry as signalled by the factor u{k). The virtual corrections are then taken into 
account by subtraction of a linear term (no real giuon in the final state, virtual corrections only 
renormalize the existing part components of the wave function) with a coefficient that ensures 
infrared finiteness through real virtual cancellation. Virtual corrections occur globally, both 
inside and outside the jet region, thus the factor u{k) is absent. This yields the second term in 
Eq. (HSU). 

That this indeed ensures infrared finiteness is best understood after rewriting l(154|) in the 
form 

EOe Gab{E) = -EdE R^a^E) ■ Gab{E) 

(159) 

asWab{k)u{k) [Gak{E) ■ Gkb{E) - Gab{E)] . 



An 

?(o) 



where R\i^ {E) is called the single log Sudakov radiator for bremsstrahlung emission 



R'^abiE) ^ jj^J ^-''(^) [1 - "(^■)] = ^-^ab, (160) 

BMS show that to good accuracy one may take u{k) to restrict the phase space integrals to 
the outside region for the Sudakov radiator term and to the inside region for the remainder of 

Eq. (iisnj- 

With integration in the Sudakov term restricted to the out region, one recognizes that there 
is no danger of encountering any ill effects from the poles of the kernel in this term. For the 
other term one needs to prove (see [72]) that Gaa{E) = 1. Then the quadratic and the linear 
terms cancel where the kernels diverge, rendering the expression finite. 
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For more details on the ingredients as well as the physics of this equation see [72,76,77,80,81] 
and references therein. 

It is now straightforward to compare the BMS equation, Eq. ()154fl , to the BK equation in its 
incarnation in terms of SY;xy'- 

dySy-xy — ^^2^ J d^Z iCxzy {Sy-xzSy-zy — Sy-xy) ■ (161) 

Ignoring for the moment the Sudakov radiator term in Eq. I|159|l . or setting u{k) 1 for the time 
being, the similarity of Eqns. (|159|l and (|161|l is striking. A one to one relationship of structures 
emerges if one maps directions onto transverse coordinates, the kernels onto each other and tries 
to view Gab{E) as the average of some ti{UlUb)/Nc (where a and b represent the directions of 
the leading hard particles): 

GabiE) ^ {tiiUlUb)/N,)E . (162) 

Note that this is fully compatible with the requirement that Gaa{E) ~ 1, needed to ensure real 
virtual cancellation in the second term of Eq. (|159|1 . While this appears to be quite intuitive 
and fully in line with the physical interpretation, it is not immediately clear how precisely to 
perform that latter part of the translation. One would ask, for instance, how in detail to arrive 
at a consistent definition of these U s and the averaging process suggested by Eq. H162(l in the 
light of the tree like branching process underlying the physics of the BMS equation. Additional 
questions would be how to reconcile their role in probabilities with that in amplitudes and many 
more in the same vein. 

That the interpretation proposed in Eq. Hlt)2|) should be possible suggests itself even more 
strongly, if one uses the observation made by BMS that in the small-angle-emission limit in a 
carefully chosen frame where the measure turns flat and the kernels agree completely: 

Wabik) -> K.^^ . (163) 

This is the region where Marchesini and Mueller [79] have established the link between jet and 
BFKL-dynamics. 



5.2 From JIMWLK to BK: generalizing BMS by analogy 

To understand what the corresponding nonfactorized equation might look like, recall the chain 
of argument leading from the JIMWLK equation Eq. (|81|l via Eq. (I91|l to the pre-factorized from 
of the BK equation (|94|l . From this short calculation it is pretty obvious that the linear and 
nonlinear terms arise in essence from the first and second term in Eq.JSSJ: 

HnuwLK =-l^ ^^^v [*V:*V^ + iVl^l + (*V^zV^ + zV^zV^)] . 

Note that this is in line with the interpretation of these terms as virtual and real contributions 
respectively. To arrive at different coefficients in linear and nonlinear terms as in Ea. H154|l one 
simply has to allow for different Kernels respectively. 

To write down the desired Fokker-Planck Hamiltonian in the jet case, one needs 

1. to use the correspondences of variables and kernels noted in the above. In this step it 
proves useful, if I slightly redefine the functional aspect of the invariant vector fields to 
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produce 5 functions adapted to the structure imposed by the sohd angle measure. To this 
end I will write 



tW;Up = f'UpSip ^ q) (164) 
where S{p ~ q) contains Jacobian factors such that 

^~S{p-Q)fi<l)-fip) ■ (165) 

This will simplify expressions considerably. 

2. The different coefficient in real and virtual parts are implemented by allowing for different 
integration kernels /C''^ i ~ 1,2 for the two terms. Clearly both analogues of /C'*) will be 
proportional to Wpq (fc) with the proportionality factor (fc) carrying the jet geometry via 
some u{k) dependence. 

The result of this exercise is a Fokker-Planck Hamiltonian that has been dubbed Hng (for "non- 
global") in [54], here I refer to it generically by iJjct- I will present it in two forms, to parallel the 
two versions of the BMS equation, H154(l and (|159|l . To this end I will introduce an additional 
function f^'^\k) that will serve as the coefficient of the generalization of the Sudakov radiator 
term. One has the following forms of the Hamiltonian: 



^jct := - g Wu. (k) [/(I) (fc) {^IzVl + iVltV:) + (fc) [Uk] + 



(166a) 



in analogy with p54() , and 



+ /^'Hfc)(*V?.*VS+*V?.*K+ [;:^fc]"''(*K*Vt + »VS*Vt))] (166b) 

to parallel (|159|l . In both cases u,v and k are integrated over according to an "integration 
convention" that uses as its measure for any of the momenta. Retracing the steps leading 
from the definition of the JIMWLK equation to (|94|l by eye should make it obvious that indeed, 
the separation of terms in Eq. (|166b|l is such that the first line generates the Sudakov radiator 
and the second the nonlinear evolution inside the jet cones. 

A comparison of the result with BMS allows one to determine only the leading Nc part of 
the /('^ Already the l/N^ corrections are not controlled by the matching. To this accuracy the 
/ are Nc independent: 

/(2)(fc) ^u{k) (167a) 
/(i)(fc) =1 (167b) 
/(i)(fc)=(l-u(fc)) . (167c) 

Eqns. H166|l . 1167|) together with a Fokker-Planck equation of the form 

EOeZeP] = ^H;,tZE[U] (168) 

^^Note that the calculations are somewhat simplified by the fact that Wpp{k) = 0. The relation between kernels 
on the functional side to those on the BMS side, as given by Eq. 1921 for JIMWLK/BK pair, becomes a matter 
of simple proportionality. 
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define the finite Nc generaUzation of the BMS equation given in [54]. One might argue that 
this has to be the correct result on naturalness reasons -any Nc dependence in these equations 
would make it virtually impossible that they should be able to incorporate all Nc corrections in 
a consistent manner- but an explicit construction of this equation from the amplitude level is 
much more instructive and will be reproduced below. 

The most important feature of this equation is its infrared finiteness which directly relates 
to the real virtual cancellations being correctly encoded in Eq. I)166b|) . Considering the Sudakov 
term in iJjet one finds that just as in the BMS limit, the phase space integration over k is restricted 
to the outside region where no hard particles are to be found. Therefore there is no divergence 
from the poles in Wuv{k). The second term may be recast in analogy to Eq. H82(l and the same 
finiteness argument used there clearly carries over to this situation. By looking at the BMS limit 
one sees that this same cancellation of terms corresponds to a real-virtual cancellation in the 
conventional sense. On the functional side, real and virtual contributions arise from different 
terms. Contributions with an additional factor Uk correspond to real emission, the rest to virtual 
corrections. 

The idea to interpret jet transition probabilities as correlators of path-ordered exponentials 
in singlet projections properly takes into account gauge invariance, just the same way as the 
formulation of the JIMWLK equation does. The equation suggested therefore would appear to 
consistently encode the physics it is intended to cover with a set of very strong constraints. 



5.3 An equivalent Langevin description 

Just as in the small Xhj situation, a translation to a Langevin type description is possible. The 
steps are almost identical to the ones presented in Sec. 01 just slightly complicated by the added 
exclusiveness, the distinction into in and out regions. 

Prerequisites are (termwise) positive (semi-) definiteness of the Fokker-Planck Hamiltonian 
in question. Only then can one introduce a bounded noise integral - if needed separately for 
each positive term. If the Hamiltonian can be (termwise) separated into two conjugate factors 
in analogy to the structure of Eq. 1)82(1 , any correlation effects can be included in the Langevin 
equation. The noise can then be taken to be (termwise) uncorrelated. 

Fortunately all of the prerequisites listed above are also met for Hj^t- In order to achieve 
termwise factorization, I first separate off the /^-^^ terms and further treat the two contributions 
proportional to /^^^ separately. All of these are positive due to the definition of u{k): 

u{k) = ei„(fc) + e-''"eout(fc) (169) 

ensures < u{k) < 1 and the same then for the coefficient functions in Eq. (|166b|l . 

To parallel the treatment in Sec l4.1l and l4.2l one would write the Hamiltonian on the JIMWLK 
level as 

^jct -^*v;ix:t*v^ (170) 

where u,v are integrated over with ^^^n^ 



ab _ / aiifc 



-1 ab 



wUk) f^'Hk){l + UlU,) + f^'Hk){l-UlUk){l-UlU,) . (171) 



Factorization of x according to x — is achieved by defining a three component structure 
representing the three separately positive terms alluded to before: 
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Correspondingly one has to introduce independent white noise for all of the components, 

1(^(1))^;.^ (^(l'))fe^;.^(^(2))^;.|^ (173) 



The only thing left to cope with is the measure and the (5-functions in the noise correlator, which 
need to be such that 



U J '^^Pk ){J J^^ql -I )/ 



ab 
Apq 



Now the measure is 



4ir 

((S^);^'^(S^)f)=47r5(c 



and one needs the correlators to read 



(174) 



(175) 



With these preparations and the identification Y ^ ln£'/£'outOiie has a Langevin equation that 
reads 



(176) 



In comparison with the JIMWLK case all sigma-terms vanish because Wpq{k) satisfies Wpp{k) = 0. 

Note again that any continuum notation is deceptive. The equation to consider is really 
a finite difference equation and upon iteration will allow for an interpretation of subsequent, 
ordered soft gluon emission. To this end it is instructive to rewrite the Langevin equation in 
terms of functional derivatives: 



d^q dVlk 



d^^^k pabifiTrj ^ "b;/j 
~^^Vk \-^y;p\ -Y:k 



^£^^^^\UY,p\ E%^W'lJ,,^}[UY■,p] 



(177) 



The operator in the last version contains two terms in its third component, the component 
responsible for the in-region. Using a somewhat more compact notation for the Y dependence 
they read: 



fdflg 


f dflk 


1^ J 


1 4tt 



i{£ 



dV,„ dVtk las p 



Air An 



TT p.q 



(178) 



Taking into account the fact that the Langevin description is defined only in a y-discretized 
sense, one might be tempted to conclude that the terms containing factors Uq correspond to real 
gluon emission, while the others would generate virtual corrections, not only in this, but also the 
other components. 

A formal solution to (|176|l is 



Uy^p = P.exp{. / dr' / -^^v^] , 



(179) 



which again is to be interpreted in a discrete sense. To write this expression I have also adapted 
the definition of the functional derivatives to include a (5-function in Y: iVypUyq — —UYpt6{p— 
q)5{Y-Y'). 
The operator 

i-Y 



(180) 
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is the analog of the JIMWLK shower operator (| 149(1 and is not only relevant to the Langevin 
formulation itself. It implements subsequent real emission with virtual corrections correctly 
taken into account. This idea will be substantiated below in a three step process: The first step 
consists of explicitly constructing the real emission part of these amplitudes and comparing it 
to known results [80,81]. The second step is the use of this information to recover the evolution 
equation (|168|) -this step reinstates the virtual corrections- and the third step is the equivalence 
of Fokker-Planck and Langevin treatment which closes the line of argument. 

5.4 Amplitudes in the strongly ordered domain 

In order to find a suitable starting point to derive results beyond the large Nc limit, one has to 
go rather far back and start with a general discussion of soft gluon amplitudes along the lines 
given already in [80,81]. 

To help organize the argument, I will start with the definition of a generating functional for 
the (tree level) amplitudes of (real) soft gluon emission from a qq pair in the strongly ordered 
region as already done in [54]"^° 

t^fc„ < t^fe„-l < . . . t^fci < t^p < UJq (181) 

written as 

n=0 '' 

The amplitudes for n-gluon emission from a initial qq pair arc denoted by Ai^qpi^-^"""!^ ), where the 
momenta fci . . . kn and color indices z, j, ai . . . a„ are explicitly listed, a corresponding set of 
Lorentz-indices //i ...//„ is suppressed. (For diagrammatic representations, see Eq. 1)184(1 .') They 
are isolated by n-fold (hmctional) differentiation w.r.t. at S = 0. At this stage, S is just 
an external source, the relationship to the noise of the Langevin description will become clear 
shortly. 

The ^(qpfc\''' fc" ) are known to satisfy an iterative structure in which the color singlet 
qqg'^^^ -aToplitvLde follows by induction from the corresponding amplitude with n soft gluons. 
As a result, the 0-gluon term in Eq. ((182(1 fully determines the whole functional, although the 
explicit construction of the n-th order expression with the tools available to date becomes more 
and more cumbersome with growing n. The rules of the game have been explicitly demonstrated 
in [81] and require no recourse to the 1/A^c limit. 

The important observation now is that there exists [54] a relatively simple, closed operator 
form for this generating functional. In fact it consists of an even more general operator acting 
on the zeroth order "hard seed" term in the sense that this operator can be applied to any "hard 
seed" to generate the full soft gluon cloud accompanying it, the full shower operator already 
presented in Eq. 1(180(1 . This will be explained in this section using the existing information on 
tree level amplitudes and hence the discussion leads to the real emission part of this operator as 
given in Eq. ((198(1 . 

For this purpose, a slight generalization of the definition ((182(1 is operationally somewhat 
more useful, in that it allows me to specify an explicit functional form that casts the iteration 
step from qqg^^^^ to qqg^^l of [81] as a simple functional operation. This operation will involve 
nothing more complicated than the functional differentiation rules given in Eq. 1(85(1 and a suitably 

■^"Compared to the treatment there S has been replaced by iB. to make the correspondence to the shower 
operators immediate. 
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defined n = term to start the process of generating the amplitudes. The starting point is the 
definition 



n=0 



P^du^], ■ (183) 



Here I have simply added a factor of U in the appropriate representation to each leg of the 
diagrams contained in the definition of the ordered amplitudes. Diagrammatically the first few 
contributions read: 






(184) 



Obviously Eq. H182() is just the special case at U = 1: A^J^lE] = A^J^[1,E]. The rationale behind 
this definition lies in the very nature of soft gluon emission underlying the construction of the 
amplitudes. Per definition, soft gluon emission does not change the direction of the parent (the 
emitter). In that role as a parent the entity (be it a quark or a gluon) is therefore written as 
a path ordered exponential along the direction of some momentum. Emission of a softer gluon, 
that in later steps will also serve as a parent will then necessarily add an adjoint factor Uk and 
an cikonal emission vertex J/^t" (with J^^ := ^), where the generator is to be taken in the 
representation of the emitting object. As all these objects are encoded as eikonal lines U, this 
will automatically occur if one writes it via invariant vector fields iV" just as in the rewrite of 
the Langevin equation at the end of Sec. 15.31 In functional form, one ends up with the single 
emission operator 

47r 47r ' ^ ' 

Note how closely this resembles the real emission part of the Langevin equation (|176|l . Eq. (|183|l 
should then be completely determined by the n = term in the sum which describes a qq- 
pair without additional soft gluons: higher orders should just follow from repeated application 
of H185(l . This needs now to be checked against what is known about the soft gluon amplitudes. 
First one needs to give an initial condition for this iteration, the bare qq term. This reads 




M2{q,p)[UlUq] 



(186) 



and corresponds to a zero order amplitude of the form Ai^^^) = M2{q,p)6^^^^ 
The above prescription then claims that the qqgsoft term is given by 

dfli dQk 



47r 47r 



(187) 



To demonstrate once how to use the machinery, I will go through the steps explicitly. First insert 
[A^°^Yp3^[U,E] from fTM|l . This yields 



[AW];^J[/,S] = {»y'. 



dfli dQk 



(188) 



^^For Mn I have adopted the notation of [80,81] 
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All that is left to do is to use the differentiation rules H85|l (with the adaptation to phase 
space ()164|l taken into account) and differentiate the two factors and Uq. The variation 
eliminates the I integral and one is left with 




+ 




-ill 




where 



Now one isolates the amplitude via a variation in S and reads off 



AC'" 



pqk > 



Af2(p,9)A^,(fc)r]y =: M3(p,A;,(?)r] 



(189) 



(190) 



(191) 



with M3{q,k,p) again as in the notation of [80,81]. Eq. (|191|l in fact is identical to Eq. (40) 
of [81] and an additional gluon emission, i.e. iteration with Eq. (|185|l . leads to Eq. (42) of [81] 
for the amplitude with 2 soft gluons: 




Uk2'i-Ek2 




Ukii'S.ki 



(192) 



(The color algebra for the explicit comparison is given in [54], App. A, case n = 1.) The pattern 
simply follows the rules of functional differentiation and a comparison of analytic expressions for 
the amplitudes as given by the iteration rules of [81] was performed explicitly using symbolic 
algebra tools up to n = 6, where the task starts to become time consuming. To compare beyond 
finite orders, a few additional tools are required. 

In the above I have implemented the ordering "by hand," simply by requesting the hierar- 
chy (|181|l for the ujki ■ To write the all orders expression advertised above, I need one more 
ingredient: the addition of an energy integral to the emission operator H185() and an appropriate 

defined such that 



extension of the definition of iV" to iV|^ 



(193) 



where S{p — q) is adapted to the measure used in the momentum and energy integrations. The 
natural choice for the latter is to replace the solid angle integrations in (|185|l by full fledged phase 
space integrations: 



dQk 
in 



(27r)4 ^ ' (27r)2 47r ^ ' 



(194) 
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d{p — q) is then defined as the appropriate S function on the forward light cone, i.e. with Jacobian 
factors such that 

[dcf>k~6{p-k)f{k)^f{p). (195) 



(Note again the similarity of the roles played by ^-ordering in the Langevin description and 
energy ordering in this case.) With these definitions it is straightforward to verify that a replace- 
ment of the solid angle integrations in the above by phase space integrations leads to no change in 
the previous results, if one implements energy ordering for subsequent definitions with ordering 
0-functions, just as one would for the contour parameter in case of path ordered exponential. To 
leading logarithmic accuracy this will not change the result for the evolution equation. 

By now it should appear to be the natural choice, if I write for the generating functional 
(again suppressing the w-label for compactness) 

A%[U, 5] =P., exp J g 5t^''*V?} [A(°)];^Jt/, S] , (196) 

where P^^^ implements the strong ordering of Eq. H181|l . I should actually use a notation that 
displays the dependence on the functional form supplied in the zero emission part and for instance 
write Aj/"jC/,S,Af2(p,(z)]. 

In analogy with path ordered exponentials, this implies that the exponential series is to be 
interpreted as 

P^, expji j gJ^^Ui" H^-'-iVr} 

(197) 

where the 1/n! are absent due to the explicit ordering of the subsequent soft emissions. For 
compactness [54] let U[S, [/] denote the operator 

U[S, U] := P., exp { J g J^^^fj-b st^^iVf} . (198) 

The explicit mention of the JJ-field in the notation will be needed later when different types of 
eikonal fields appear. U[S, U] has a meaning independent of the bare (n = 0) term it is used to 
act on. It is called the (real) shower operator and arises wherever one is forming that ordered 
soft gluon cloud around a hard seed. 

With these tools a direct comparison with the iterative scheme of [81] at all orders becomes 
possible. In essence all that is left to confirm is that the iteration in the shower operator leads to 
the same result as the induction step from n to n + 1 gluons in [81]. This has been done in [54]. 

It is already apparent that the structure contained in the above is very similar to the informa- 
tion contained in the Langevin equation. This has been used to derive the RG equation stated 
earlier [54]. The main observation in this regard is that, as a consequence of the strong ordering, 
any in which k is at (or near) the phase space boundary, there will be only a contribution 

from the first emitted gluon (n = 1). As a consequence, for such variations one gets (h.b. stands 

•^^Also the results below for the structure of the multiple ordered soft emission can be obtained both ways - to 
have a formulation with all ingredients (including the ordering) written explicitly is merely a convenience from 
this perspective. What the formalism does, is to allow to give a closed expression for the generating functional, 
Eq. GTO}. 
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for hard boundary) 



h.b 



Note that 



(199) 



(200) 



is the qqg counterpart to Eq. p96|l . In short, Eq. 1)199(1 states, that under variatfon at the phase 
space boundary, one relates the generating functional for a given tower of amplitudes to that of 
a tower with an additional hard gluon. For the above example it relates the amplitudes with a 
hard qq to those with hard qqg content: 



A;/JC/, S, Ahip, q)] = Af^,[U, S, A^^^{k)M2{p, q)] 



h.b. 



(201) 



or diagrammatically 



h.b. 





(202) 



These two items, the ordered nature leading to ((199|l under "hard" variations, and the fact that 
in this situation the number of hard legs of the soft emission amplitudes is increased by one, as 
exemplified in (|202|l . will be the core observations behind the evolution equation I am aiming at. 



5.5 Transition probabilities from H averages 

As a first step to understand how to translate the above into expressions for transition proba- 
bilities, let me take the situation of a (real) qqg^^f^ soft emission amplitude in which one must 
obtain the standard result 

n„+2(0 n„+2(i') 

xtr(t'"i t'""^"'^ f'l), (203) 

with gluon momenta in the final state labelled by qi, . . . Indeed, the same expression is 
obtained, if one considers the S averaged product of functionals 

([A(")];^J[/,S][A(")];^J[/,S]t)„. (204) 

if the average is defined with a Gaussian distribution for S in which only strongly ordered modes 
gi, . . . , g„ have support. I write 

(. . .)„H := det(M„)^ j D[~] . . . e'H rf*pd0,H;A/„-/^,H, ^205) 

where d(j)k denotes the phase space integral for momentum k and I have suppressed discrete 
indices. M defines the S correlator, which in the fully exclusive case considered above would 
read 

n 
i=l 
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In this case, because of the strong ordering of momenta in the amphtudes, the 2 integral will 
strictly pair ojf the one available momentum ki in the amplitudes falling into the same range as 
Qi with the latter. 

As such, this is only true if I have ordering in the emission vertices, either done by hand 
in the 2-d version, or via P^j^ in the 2+1-d formulation. Only then, unwanted cross terms are 
excluded. The U factors just cancel trivially. 

To make contact with the non-global observables of the BMS setting, one needs to depart 

n 

from the fixed n situation and replace the ^ S{q — qi), which selects a given set of final state 

1=1 

momenta in the case where the number of final state gluons is specified, by a cutoff criterion 
that restricts real emission by geometry and energy. By writing 

M^f = {E^'^E^n^ = S^'g^'^Eip - q) 6{E - LO,)u{q) , (207) 

one can indeed go beyond a final state with a fixed number of gluons and consider 

(A;yc/,s]Aj,yc/,s]t)H , (208) 

which will provide the finite A^c-generalization of g'^^^^^^ of BMS. 

In fact any degree of inclusiveness may be imposed by modifying the restrictions on the 
allowed modes appearing on the r.h.s. of Eq. H2U7|) . 

This realization should make it clear how to include virtual corrections right from the outset, 
by a fairly simple modification of the above. I will refrain from doing so here and instead keep a 
closer parallel with the strategy employed by [72,80,81] and use the expressions for real emission 
to deduce the structure of the evolution equation. 



5.6 Evolution equations for soft semi-inclusive quantities 

For semi- inclusive quantities like the non-global observables of [72], in which one sums over a 
given, limited phase space volume of soft gluons, it is natural to ask for the dependence on that 
phase space boundary. Unlike a direct calculation of the average in Eq. H208() this should lead 
to a tractable result that exhibits new structure. In fact, one recovers the generalization to the 
BMS equation suggested above. 

To arrive at an evolution equation, I will simply take a derivative w.r.t. the phase space 
boundary of the emitted gluons in the semi-inclusive probability (|208|l . This is in direct corre- 
spondence with the procedure used in BMS. Here it is essential that one includes the phase space 
boundary into the definition of M . Then the result is most transparently displayed using the 
following simple relationship, based on functional differentiation and Legendre transformation in 
the case of a Gaussian action ( "free theory" ) : 



, _ .t*m:t 



{W[E])^ =N I D[E] W[E] e-^^ = w[-jj]e~'^ ^ ^ ^ ^ ~" W[Eo 



J=o 



Ho=0 

(209) 



for any functional VK[S]. Here notation has been condensed even further, with all momenta, 
integration signs and measures suppressed. '^^ Legendre machinery is used for the last equality 
sign. The "classical field," Eq = iMJ', vanishes at J' = 0. 



•^■^For compactness I will do so throughout this section. Whore indices and momenta are shown it is with the 
understanding that there will be integration conventions for repeated momenta with d<I>fe as the measure. 
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The canonical example for W['E] here of course is 



{W[E])^ ^ {A^^\U, S]Ajf jC/, S]t)^ , (210) 

which, by construction, is the finite Nc generalization of G^^^^^\ 

For this example one immediately realizes that the exponential of the second order differential 
operator on the right hand side will perform precisely the "sewing" implemented by the Gaussian 
weight in Eq. (|208|l . as long as M is diagonal in energies as in H207() . In this case, the energy 
ordering ensures that 

1 S ..6 ... ... ^ _ ^ nrf^A^l ^it 



Let me now carry out the derivative with respect to the phase space boundary in analogy 
with the derivation of the BMS equation sketched in Sec. 15.11 With M = M [E) a function 
of the phase space boundary, one immediately finds a general expression for the logarithmic 
-E-derivative of the above expectation value that reads 

EdE {Wn)^{E) = ^-e-^^- ^ °EdEM{E)-^ W[Eo\ ■ (212) 

The main point about this seemingly trivial exercise is its use in conjunction with what is 
already known about functional derivatives of the amplitudes of interest: Eqns. (|199() . (|201() . 
Note that OeM (E) will force the variations in the factor taken down from the exponential to be 
at the phase space boundary. Choosing W[E] = Al^glU, E]A'^Jg[U, S]^, the square of the generating 
functional of qq g^^f^. amplitudes, this location at the boundary will lead to the appearance of the 
square of the generating functional of qqg g^^f^ on the right hand side as follows from H2U1|) , (|202|) . 

This marks an important difference compared to the BMS case at infinite Nc'. the equa- 
tion does not close. Instead higher and higher hard correlators will enter: Clearly, an evolution 
equation for the qqg g^^^. amplitudes will then couple in turn to the corresponding object for 
qqg^ ff"oft amplitudes and so forth. One is faced with an infinite hierarchy. All quantities encoun- 
tered in this hierarchy can be defined in complete analogy to the examples presented in H196|l 
and H2()()(l . One simply uses the desired combination of L/- factors as the "bare" term in the tower 
of amplitudes one is interested in. 

This implies that firstly, it is not sufficient to consider only the evolution equation for qqg^^^^ 
amplitudes alone. Instead one has to capture the complete infinite coupled hierarchy of equations 
in one. Secondly, it indicates that such a step is possible, since the objects of interest have already 
been identified and indeed can be written down in a general form. The evolution equations then 
follow from ifH^ . 

I therefore proceed to define the (tree level) generating functional based on m hard particles, 
a general "antenna pattern," by 

a;v;;X[c/,s]a^V.vX[^'S]^ = u[s,c/] {uv^)^^} ® ...®{uv^)^^l u[s,y]t|^^^. (213) 

The factor (t/F^)^!^ ®...® (C/V^^)^!l in this expression generalizes the expression for the bare 
particles of the qq case (the Ap^^[[/, 5]Ap-'g[C/, S]^^ from above) from 2 to m hard legs, and the 

U[S, J7] ... U[S, implement the showering. I have given only q and q factors as any 

gluon can be written in terms of these as {UV^f^ = 2tr(i"([/yt)t t''([/yt)^,). 
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One can now study the evolution of the expectation value of a given object of the type (|213|l 
or directly study a generating functional for all such objects in one go. The latter is given by 

\N[j\j] :=(U[S,C/] e»2tr/c/yt+,2tr(c/y+)t, [j[E,V]^) _\^_^ (214) 

with integrals in the source exponents understood. Appropriate variations with respect to j and 
at j = will then select a given tower of in hard g^^f^ probabilities entering the hierarchy. 
Since, besides the phase space constraints contained in M, it is the particle content of the bare 
terms that define the amplitudes, a notation that emphasizes this is needed. To this end I will 
write 



S 



=■■{ iuv^)l\^ ® . . . ® (c/W)m yj^'^ (215) 

where the UV^ average in the second line is to be taken with a yet unknown weight that re- 
produces the l.h.s.. The idea is to interpret the soft gluon average as an average over UV^ 
configurations. This is completely legal, if somewhat formal as long as the only definition of the 
weight of the average is through (|2ir)|l . 

Nevertheless, it is this reinterpretation that will allow to make contact with the Fokker-Planck 
formulation of Sec. 15.21 and that is all that is really needed. 

For m = 2 and an adapted choice for the color structure one encounters 

{tT{{UV^)l{UV^),)/N,y^^, G'-r'\E) (216) 

and thus provides the precise definition of what was still a bit vague in p62|) . Eq. H215|l is of 
course much more general than that: it defines the real emission part for a general antenna 
pattern, as these will all be needed for the evolution equations. Inclusion of virtual contributions 
will be discussed below. Note that the eikonal lines appearing in l|215() have the interpretation 
of a product UV^ of contributions from both the amplitude and the conjugate amplitude. 
With (pn?^ one then has 



W[j\j] =(U[S, U] ([/yt)(t) ^ . . . ® (t/yt)(t) u[S, V]^)__ 



u=v 



u=v 



_^gi2trjt(7yt+i2tr(;7yt)+j^''cal (217) 

This demonstrates that ,j] is in complete analogy with (the real emission part of) Z[J, 
in the derivation of the JIMWLK equation of Sec. 13.11 

To find the evolution equation for this (meta-) functional and with it the infinite hierarchy of 
evolution equations alluded to above, all that is left to do, is to put together Eqns. H212|l . H211() . 
and H199(l to get the dependence on phase space boundaries for arbitrary [/^-correlators (bare 
m-jet probabilities). This yields the r.h.s. for the real emission contribution: 

e- '°" ^ '^" U[So,£/] 

{g JilUf^^n^ } [EdEM{E)]1l':^^, {g ^'..V^"' ^l^, } e«2tr,^c/yt+.2tr(c/yt)t,- 

U[So,y]^ . (218) 
u=v 
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To arrive at a meaningful answer one still has to supplement virtual corrections. They are 
correctly incorporated as usual, by simply subtracting the corresponding term proportional to 
the original transition probability. This was done in the derivation of the BMS equation and 
it also applies to the full hierarchy. For concreteness, and without loss of generality, I will 
demonstrate how to achieve this for the case of m hard particles shown in Eq. H215|) . 

To understand the constraints on the virtual corrections, let me first study the real emission 
contribution in some more detail. Omitting, for the moment, the operator that implements the 
soft gluon shower and the expectation value in the expression for the real emission part. 



5_ 



u=v 



-u[So, [/]... u[So,y 

one has to consider 

{gJt,Uf^n][E^EMiE)]u'S'''g,,,{gJ[:lv::'\y^^^^^ ([/0)(t) ^ . . . ^ (f/yt)(t) (219) 

and evaluate the variations. Inserting the BMS choice (|207|l for M, it is easy to collect all the 
factors and to make the k integration explicit (I keep the integration convention for / and V for 
readability) : 

^E5{E-uJ,)^^wa'{k){Uv')t^V'iJ^^V'y. (t/yt)(t) ® . . . « (c/yt)(t) 
'^^'^ '''wn'{k){UV^)ft^tj,^^'' (C/yt)(t) ® . . . ® ([/yt)(t) . (220) 



In this expression ui and will be made hard when the functional derivatives hit the U and V. 
The energy part of the corresponding phase space integrals will be eliminated at the same time. 
In summary one will end up with 3 "hard" momenta k,l,l' . 

Carrying out the differentiations, a sum over pairings will emerge, with different color struc- 
tures depending on whether one hits q or q lines. Concentrating on the individual terms in 
this sum, it is obvious that one has to deal with 4 different combinations, corresponding to the 
pairings qq, qq, qq, and qq completely in parallel with Eqns. of the JIMWLK case. 

At this point it is, however, more efficient to directly use the unifying properties of the left 
and right invariant vector fields without exploring these structures explicitly. This may be done 
by noting that Eq. (|219(l allows for a very elegant rewrite: 

{UV^)g^ <»...®{UV^Y^l . (221) 

Eq. H221|) follows from a direct correspondence of iV^j^^t^^ with «V^^ acting on the amplitude, 
and with «Vy^ acting on the complex conjugate factor. In Eq. 1221() it is possible to 

anticipate that the energies will all be made hard by functional differentiation on the hard factors 
{UV^)p- and thus to go back to the original definition of the variations without the energy S- 
function. This allows to use an integration convention that employs the solid angle integrations 
only and a symmetric treatment of k,p, q in this equation. Reinstating the shower-operators and 
the expectation value, one arrives at the real emission part of the evolution equation 

{^Wp,ik)u{k)iUV^)t\^V'l^y,^^^V\^y,^^ + ^V^^^,t),^V^^^t) J 

{UV%\^®...®iUV^)gl)^^^iE) . (222) 
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Virtual corrections do not depend on the physical process at hand: While the real emissions 
are confined to the inside regions and hence carry a factor u{k), virtual corrections appear 
everywhere and therefore have no such factor. Moreover the factor {UV^)l'' will be replaced 
by S°'^ and the vector fields will have to act twice within cither the amplitude or its complex 
conjugate, leading to zV^j^^t^ iV°j^^^-^ and an analogous barred contribution instead of the 
mixed ones in Eq. H222|l . For symmetry reasons, one therefore expects the virtual corrections to 
read 



(223) 

where the overall sign and normalization is fixed by real virtual cancellation: only with this choice 
do they add up to an infrared finite Fokkcr-Planck Hamiltonian, which happens to coincide with 
the expression conjectured in Sec. 15.11 Indeed, the operator appearing in the sum of (|222|l 
and H223() is nothing but — i?jet of Eq. H166a|l with the completely determined. In particular, 
there is no room for any Nc dependence of the coefficients from this argument. 
To summarize, it has been shown that 

EdE (Je*2tr,*C/Vt+,2tr(C/yt)t,\^^^^^^^ ^ _^ ^^^^ g.2trjtC/yt+«2tr(C/yt)tj ^^^(E) . (224) 

The only argument still missing is about how to reverse the step leading from H81|l to (|90|l . 

This follows from the fact that the identity holds for arbitrary j and , i.e. the fact that the 
above is completely independent of the type of correlator considered. The result is, as advertised 
in Sec. 15.21 the evolution equation 

EdE Ze[UV^] = -i/jet Ze[UV^] (225) 

for the weight Ze\UV^ used to define the averages (. . .)uvi{E). This completes the derivation 
of Eq. mia . 



5.7 The two types of nonlinearities; Nc and correlator factorization 

Now that the interpretation of Gpq{E) in terms of link operators suggested in (|162|l is understood, 
it seems worthwhile to recapitulate the role of the two types of nonlinearities discussed repeatedly 
already in the JIMWLK case: eikonalization of gluon fields into Wilson lines which enter the 
definitions 



Gpg ■.=tr{{UV^)liUV^)g)/N, (226a) 
1 S_ 



Gp,{E) :^{Gpg)uv^{E) = — ,,,, , ,,,, , W^[/,j] , (226b) 



as the first stage and nonlinearities in the Wilson lines themselves contained in the Fokker-Planck 
Hamiltonians (|166b(l or as "second eikonalization" in the shower operators H18()|l as the second 
stage. 

Jet observables provide the perfect example that the second stage is not always visible; 
in global observables u{k) suppresses the inside contributions completely and only the outside 
terms survive. This completely suppresses the nonlinear terms. How this happens is again best 
illustrated with the evolution equation for the two point function. Starting from (|225|l . one 
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immediately finds 

EdE{triiUV^)l{UV^),)/N,)uv^iE) 



(227) 

The first term on the r.h.s., carrying the factor u(fc), originates from the real contribution, the 
second from the virtual one. By H93fl this is equivalent to 

EdEGp,{E) = J ^^Wpg{k)(^u{k){GpkGkgN, - ^) - 2GfGpg){E) (228a) 

^asWpg{k)(^u{k){GpkGkg - Gpg) - ^(1 - uik))GpgyE) . (228b) 

First and second terms contribute inside and outside the jets, the former is completely suppressed 
in global jet obscrvablcs. In this case, all traces of the second nonlincarity arc gone: since the 
equation is linear, nothing couples to higher n-point functions. The counterpart to the Balitsky 
hierarchy fully decouples. What is left is simple Sudakov exponentiation which corresponds to 
the fact that G still corresponds to a correlator of Wilson lines. 

On the other hand, this formulation has a direct link between large Nc limit and factorization 
of correlators which leads to the analogue of the BK equation that triggered these explorations, 
the BMS equation. To see this explicitly, use (|2()9|l to represent the averaging involved in 

Gp,iE) = {Gp,)uvf (E) = {A%[U, ^A%[U, E]% (229) 

and then repeat the argument of [81], Eq. (50). It is worth emphasizing that the real- virtual 
cancellation argument of BMS restated in Sec. applies already to Eq. (|228|l itself, without any 
reference to the factorization argument. Infrared finiteness in the Sudakov term follows directly 
from the argument of BMS, while in the other term it is the definition of G that guarantees the 
necessary cancellations. This provides an illustration of how the more general reasoning above 
will lead to infrared finite results via real- virtual cancellation in any evolution equation for color 
singlet objects contained in the new evolution equation. 



5.8 Connections with CGC physics 

Here I want to discuss two types of issues concerning the relationship of the treatments of 
evolution in the CGC case and that of non-global jet observables. The first concerns what one 
can learn from the above for the interpretation of the JIMWLK equation and the extension of the 
underlying methods to describe new, more exclusive observables. The second type of questions 
concerns the treatment of jets that are created in a dense medium such as in heavy ion collisions 
at RHIC and the LHC. 

Turning to the first of these, let me begin with the observation that the two types of terms in 
the JIMWLK Hamiltonian indeed have the interpretation of virtual correction and real emission 
terms. Any restriction on (real) soft gluon emission imposed by a choice of variables must 
come along with a modification of the phase space available on real emission and introduce 
the analogue of the factor u(k) appearing in the jet case. This analogy is important with 
respect to a generalization away from purely central collisions in a way that remains under 
control perturbativcly. The problem is conceptual and affects both JIMWLK and BK equations, 
although in the latter case the problem is easier to state. To this end, return to the discussion of 
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IR safety of BK (and JIMWLK) evolution as illustrated with Fig. |5| IR safety was ensured by 
the presence of a saturation scale, which in turn emerges from the existence of a dense gluonic 
medium. This can confidently be expected for central collisions, but will certainly fail as soon 
as one starts to probe the targets periphery. In this situation the gluon density becomes small, 
the evolution equation reduces to the BFKL equation with all its problems of diffusion into the 
infrared. The culprit for this in a practical sense is the long range nature of the emission kernels 
K. in both JIMWLK and BK. As a consequence, if one attempts to prepare an initial condition 
to JIMWLK or BK that has large densities only in a finite spatial range, evolution will not 
respect this structure - the finite size object will not stay finite but grow exponentially [82,83] 
. The situation for the saturated region is more subtle [84], but that does not change the basic 
problem of exponential growth. Any calculation that would lead to finite range of the JIMWLK 
or BK kernels would necessarily have to take confinement into account, an idea somewhat at 
odds with the pcrturbative philosophy taken in this review. The same applies to attempts to 
phenomenologically treat the impact parameter (or b) dependence of real experiments using 
6-profiles [50-52], although the results are intriguing. 

The above treatment suggests a different approach. In the jet case, space is predominantly 
empty, in fact the observables (be it in the non-global or global case) are designed around for- 
bidding radiation above certain thresholds into most of space. The radiation veto then leads to 
Sudakov logarithms that govern the mostly empty regions, while the JIMWLK or BK nonlin- 
earities are limited to the regions that are allowed to become dense by choice of observable. The 
only trace of confinement physics entering this perspective is parton hadron duality, and often 
nonperturbative power corrections, which can in principle be understood. All of these problems 
appear to be entirely separate from the evolution equations. A highly attractive alternative to 
the DIS-type treatment existent in the literature would thus be the definition of semi inclusive 
observables in the small Xbj case of the type of diffractive and rapidity gap events that treat the 
empty regions in a manner similar to the jet case. 

While this remains a task for the future, the above developments have made much of the in- 
terpretational lore that dates back to Gribov [85] to emerge transparently in the actual technical 
treatment. While the JIMWLK equation is an equation for the cross section, the objects entering 
the formal solutions of the corresponding Langevin equations, the "shower-operators" -Eq. (|149(l 
for JIMWLK and Eq. 1)180(1 for the jet case- give access to the underlying amplitudes. These 
amplitudes result from applying the shower operators to any type of hard "seed." The proce- 
dure simply adds fully developed clouds of soft gluons to the hard objects used as a "seed" as 
exemplified by DIS at small Xbj , where the amplitude generated by acting on a hard qq pair may 
be thought of as the soft gluon part of the photon wave function. This, strictly speaking, is the 
justification for the description of the structure of the virtual photon wave function in Sec. 11.21 
A derivation of the JIMWLK equation or an evolution equation that describes diffractive events 
with rapidity gaps^^ from this perspective, based on amplitudes and with a full understanding 
of the operatorial content would be highly instructive. 

On the other hand there are direct applications to jet measurements if the jets are created 
in medium. This is the generic case for jet measurements at the LHC. Let me again consider a 
measurement of the inter-jet energy fiow H152() but with the jets created in medium as depicted 
in Fig. ^] The first obvious modification is that only a subset of the soft gluons that would be 
measured without the medium will reach the detector, the others are absorbed by the medium. 



See [86] for an early BK based attempt. 
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Figure 13: In medium analog of the two jet event depicted in Fig. 1121 Left: Soft gluons in the inter jet 
region Cout now may be totally absorbed by the medium or be created with more energy than allowed 
in the out region but then softened until they fall below i?out- Right: If the medium is large enough one 
of the jet cones may even be completely be absorbed. 

Instead of (|152() . the best one can hope to measure is 




where E'loss is the net energy loss into the medium. Phase space logarithms now are of the type 
\n{E / {Eout ^ E\oss)) and generically larger than in the previous case. Note that one can not lower 
E'out below E'loss once the medium is present. Also the evolution equations would change as new 
contributions arise. Consider for example a gluon generated inside the medium with an energy 
somewhat above Eout which loses enough energy to fall below i^out before it reaches the detector. 
Also this mechanism enhances the importance of the nonlinear effects. [Technically this requires 
a redefinition of the soft emission kernels and thus of the expressions entering H158() .] 

Upon considering to transplant other global jet observables like thrust or broadening from 
the e+e" case to the in medium situation sketched in Fig. El one realizes that global sums like 
in 1)150(1 are not experimentally accessible as part of the contributions will be absorbed in the 
medium. This would appear to indicate that also other observables which are usually considered 
global will acquire non-global characteristics. This should lead to nonlinear evolution equations 
also in such cases. 

The issues raised here pose interesting and important questions to be dealt with in future 
work - for the remainder of this review I will return to physics results from the existing evolution 
equations. 



6 Numerical results at fixed coupling: JIMWLK and BK 

6.1 Numerical results from JIMWLK simulations 

While there are many numerical studies of the CGC based on the BK equation [61-69] there ex- 
ists only one set of simulations of the full JIMWLK equation by Rummukainen and Weigert [35] . 
Fig. 1141 shows the evolution of the dipole operator obtained there, averaged over 8 independent 
trajectories on a 256^ lattice, from interval agY = ln(l/a::bj) = 0. . .2.5.'^^ The main quali- 
tative observation here is that the dipole operator soon approaches a specific functional scaling 

''^In [35] one uses volumes 30'^-512-^, with 8—20 independent trajectories (with independent initial conditions) 
for each volume. This amount of statistics proved to be sufficient; trajectory-by-trajectory fluctuations in the 
dipole operators arc quite small (i.e. already one trajectory gives a good estimate of the final result). 
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form, where the shape is preserved but the length scale is shortened under the evolution as dis- 
cussed in Sec. 13.41 This is clearly visible on the logarithmic scale plot on the right: the Gaussian 
initial condition quickly settles towards the scaling solution which evolves by moving towards left 
while approximately preserving the shape. Such had been seen earlier in numerical studies of the 
BK case, non-factorizing contributions do not seem to change that feature. With scaling estab- 



256" R = 22 



256 =22 




Figure 14: Evolution of the dipole operator from the JIMWLK equation on an 256^ lattice, shown on a 
linear (left) and logarithmic (right) distance scale. The general pattern as predicted from the existence 
of the fixed point and illustrated in Fig. |7| shows up clearly. On the right one sees a scaling form emerge: 
a stable shape of the curve that merely shifts to the left unaltered. [Plots taken from [35]] 



lished, one can now measure the evolution of the saturation length scale Rs- There are several 
inequivalcnt ways to measure Rs from the dipole function; perhaps the most straightforward and 
robust method is to measure the distance where the dipole function reaches some specified value. 
One defines Rs{Y) to be the distance where the dipole function N{{x — y)^) = tr(l — U^Uz) /Nc 
reaches a value c: 

NiRl) := c (231) 

where c should not be much smaller than 1. The precise definition of course is c-dependent - 
different choices will lead to a rescaling of the value of i?s - a common choice is c = 1/2. From 
Rs{Y) one obtains the instantaneous evolution speed 

A = -dyl-aRsiY) . (232) 

Note that contrary to Rs this quantity is unique and independent of the convention for c in 
Eq. H231|) inside the scaling regime. Outside this uniqueness is lost. It was by measuring A 
that [35] discovered a problem with the continuum limit in both fixed coupling JIMWLK and 
BK simulations. 

While the infinite volume limit constitutes no problem -this is simply the saturation scale 
at work to protect against contributions from the IR- the continuum extrapolation turns out 
to be hard. This is shown in Fig. ^1 which compares the fixed coupling BK result with the 
values taken from the JIMWLK simulations. Any attempt at a continuum extrapolation from 
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• 0.1 0.2 0.3 0.4 

a/R ~aQ 

Figure 15: Approach to the continuum limit of A/a against a in units of Rs- The dashed line shows 
the continuum value obtained from a simulation of the BK equation. A continuum extrapolation from 
this alone (continuous line) appears unreliable. 

these data would appear to be inherently unstable. It is simply the necessity of having a very 
large hierarchy of relevant scales a <^ Rs < L that interferes with a continuum extrapolation. 
This has been found not to be specific to the JIMWLK equation but to equally affect the BK 
equation [35] . Comparison with BK has allowed to extend the systcmatics of the continuum limit 
more carefully using a much larger range of volumes and initial conditions with reduced errors and 
will be discussed in Sec. 16.31 after reviewing the size of violations of correlator factorization (|97|l 
that distinguishes JIMWLK from BK. 

6.2 Size of factorization violations 

A feature that allows a direct comparison of JIMWLK and BK simulations is the smallness of 
factorization violations in the simulations of [35]. 

One might expect that the BK equation to be a fairly good approximation of the JIMWLK 
equation, in particular if the initial conditions are uncorrelated and correspond to factorized 
correlators for which 

{tY{UlUMUlU^)) - MUlU,)){tr{UlU,)) (233a) 

or 

jtrjUlUMUlU^)) - (tr(C/t£/.))(tr(£/t£/.)) 
(tr(C/^C/.))(tr([/JC/,)) 

vanish. This would appear to be essentially true for generic separations as shown in Fig. 1161 
There, for concreteness, points are chosen so that {x — z) \\ ei, {z — z) \\ 62 and |a3 — z| = \z — z\. 
The left shows (|233a() and one should recall that the natural magnitude for these correlators is 
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~ 1. The right shows a plot of (|233b|l based on the same simulation results with one leg kept at 
length 10 in lattice units and the length of the other leg on the horizontal axis. Evolution would 
appear to erase scaling violations in this plot. Thus, given the initial conditions used in [35], 
it would appear that the correlators cancel with 1-2 % accuracy, with the violations staying 
roughly at the same size throughout the Y interval covered. 




r= 1x1 = 13.1 



Figure 16: Factorization violations for a 4-field correlator. The geometry is explained in the text. 
Unnormalized correlator 1233a|l on the left (from [35]). Typical error size is indicated by dashed lines for 
a curve in the middle. Normalized correlator Il233bll on the right. The dots indicate the value of Rs{Y) for 
the individual curves. Factorization violations are generically small and indicate that evolution does not 
induce strong violations of factorization. Consequently, the BK equation will be a good approximation 
as long as the initial condition does not contain strongly nonfactorizing contributions. 

Recently, however, lancu and Mueller [87] have argued that evolution builds up correlations 
that violate factorization. They argue that factorization violations are caused by rare configu- 
rations, build up during evolution and have the effect of reducing the evolution rate. Although 
this deserves further study, explicit comparison of shapes for the dipole correlator from both 
JIMWLK and BK simulations at same size lattices also give only very small differences. In 
particular, as shown in Fig. ^1 one does not see strong deviations between evolution rates of 
JIMWLK and BK simulations on comparable lattices. This seems to indicate that these simula- 
tions in general were not sensitive to the correlations referred to in [87] . For the present purposes 
it would appear to be safe to conclude that for the type of initial conditions explored to date, 
the BK equation should give a good approximation of the salient features of the evolution. This 
is the basis on which BK simulations have been used to refine the understanding of the UV 
extrapolation by augmenting Fig. 1151 with results from BK simulations in which one also uses a 
lattice representation for transverse space to parallel the treatment of the JIMWLK equation. 

6.3 Continuum limit of BK and JIMWLK at fixed coupling 

In order to study the cutoff-dependence and to facilitate the comparison with JIMWLK simula- 
tions, one best solves the BK equation with the same type of discretization in transverse space 
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that is necessary to implement the former. This is precisely what has been done in [35]. 

Independent results for the continuum have been obtained by another method based on 
Fourier transformation presented in [61,68]. In these approaches the kernel involves only a 1- 
dimcnsional integral for which a logarithmic momentum variable S '-^ In fc^ is used. Consequently, 
it becomes easy to include a very large range of momenta in the calculation; typically, more than 
20 orders of magnitude are used. Thus, this allows to obtain a cutoff-independent continuum 
limit result for the evolution exponent A. The key result obtained with this method was already 
shown in Fig. E| for comparison with the JIMWLK extrapolation. The large range of momenta 
needed to get convergence is already a clear indication of what can be quantitatively seen in the 
discretized coordinate space version I will turn to now. 

Fixed coupling, volume 4096" 

I 1 2.5, 



fixed coupling continuum limit 




a/R ~aQ 



Figure 17: (From [35]) The evolution of N{r) on 
a 4096^ lattice, measured along on-a:xis direction. 
The curves are plotted with interval 5Y = 0.25/a, 
and the total evolution time is 7/ a. 



Figure 18: (From [35]) Continuum extrapolations 
for A using JIMWLK and BK simulations. The UV 
cutoff is plotted in units of the physical scale Rs- 
The continuous line is A/a^ = 2.26 - 2.0{a/ Rs)°-^^ 



Fig.ll7lshows the evolution of N measured for a 4096^ lattice. The similarity with the dipole 
function obtained with the JIMWLK equation is striking, see Fig. E| Again the initial shape 
settles towards the scaling form, which then continues to evolve by shifting towards smaller R 
while preserving the shape. Nevertheless, even with a scaling form reached, the spacing of curves 
and hence A is not constant. This again is due to the growing influence of the UV cutoff. 

The most sensitive quantity turns out to be X{Y), which in this case is free of statistical errors. 
A function of the form A/a^ = ci + C2{a/RsY is an excellent fit to the data and fixes ci = 2.26, 
C2 = 2.0 and v = 0.32. The a ^ result agrees perfectly with the momentum space result 
although this was not included in the fit procedure. This is shown in Fig. ^1 The agreement 
of BK and JIMWLK results at corresponding cutoff scales is remarkable over the whole range. 
Most striking is an evaluation of the active phase space. While on the IR side activity is limited 
to within one order of magnitude of the saturation scale Qs , one needs the cutoff on the UV side 
to be about 6 orders of magnitude larger than Qs to get within 1% of the continuum value, as 
is easily extracted from the power law fit. 

This situation is clearly unnatural: large open phase space of this magnitude usually generates 
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large logarithmic corrections that need to be resummed to get physically reliable results. This 
has been used to argue that the main source for such corrections are due to running coupling 
effects and study the effects which, among other new qualitative features, will bring a large 
reduction of active phase space in the context of the BK equation that erases the potential for 
large logarithmic corrections. 



7 Running coupling effects in the BK equation 

7.1 Running coupling and UV phase space 

If one wants to discuss running coupling effects and how they affect the UV behavior one has first 
to state that JIMWLK and thus also BK have been derived in leading logarithmic approximation 
and thus, contrary to state of the art BFKL calculations, one does not know what sets the scale 
in these evolution equations from a reliable calculation. One is therefore left with a discussion 
of the structure of the diagrams leading to the driving terms of the evolution equations to argue 
a reasonable choice. For simplicity [35] have argued and calculated in the context of the BK 
equation. One should always keep in mind that qualitatively these arguments will have similar 
implications also for the full JIMWLK treatment. 

In the context of both JIMWLK and BK it is common to argue that a good scale choice is 
Qs{Y), since it is there that the r.h.s. of the evolution equation peaks. Since with this choice 
the scale in the coupling changes during evolution, this has often been dubbed running coupling 
(c.f. [55]). In fact, it is quite easy to see, that by a simple reparametrization of the Y axis 
such a scale setting is mapped on a fixed coupling treatment, provided one uses the same initial 
condition. For this reason it might be less controversial to use the term scale setting in this 
context. To be explicit, fixed coupling and scale setting at Qs differ by an additional factor of 
f(V) = "°^^'fxld"'' on the r.h.s. of the BK equation. If one uses Y to parametrize evolution with 
the Qs scale setting and Y' in the fixed coupling case, f{Y) is in fact the Jacobian dY' /dY of 
the reparametrization 



This immediately leads to the conclusion that the phase space structure of the evolution equation 
can not differ between the two and thus, this scale choice will not cure the problems encountered 
above. 

The difficulty can be traced to the fact that the scales in a diagram describing a single 
emission step during evolution for a fixed size parent dipole will generically have nothing to do 
with Qs, if that parent dipole has a size very different from 1/Qs- The situation in the IR and 
UV (relative to Qs) will behave very differently: While the nonlinear effects will generally block 
evolution in the IR, UV scales that dominate the diagrams will become visible and enter the 
coupling. While on this level one can not calculate (even for BK) which of the scales (the sizes 
of parent and daughter dipoles) or combination thereof truly enters the coupling, a scale of the 
order of the parent dipole size appears sensible. Like the authors of [73,74], also [35] choose 
to account for running coupling effects by having the coupling run at the natural scale present 
in the BK-cquation, the size of the parent dipole {x — y)^. Such a treatment is often called 
"DGLAP improved" if the term "running coupling" is already occupied for what is here called 
Qs scale setting. Explicitly in Eq. H98|l . one replaces the constant as by the one loop running 



Y 




(234) 



74 



coupling at the scale \l{x — y)'^\ 

as-^as{l/{x-yf) = — ^ with /3o - (lliV, - 27V/)/3 . (235) 

A in Eq. I|235(l is related to Aqcd by a simple factor of the order of 2tt. The equation then reads 

dyNy.xy = ""'^^^^^'/^'^^^ Jd^Z iC^.y [Ny,^, + Ny,,y - Ny,^y - Ny.^.Ny.^y'^ . (236) 

Such a simple replacement, of course, is not what we expect to arise from an exact calculation in 
which the two loop contributions that generate the running will emerge underneath the integral 
over z. This comes on top of the fact that at next to leading order there will be corrections beyond 
the running coupling effects that are entirely unaccounted for by this substitution. Nevertheless, 
it appears natural to argue that Eq. H235() . when used in the BK equation, captures the leading 
effect associated with the "extraneous" phase space uncovered quantitatively in Sec. 16.31 This 
has been the outcome of a long discussion in the context of the BFKL equation and its next 
to leading order corrections [88-94]. Here one simply takes this as the starting point with the 
attitude that a scale setting technique like the one employed in Eq. (|235|l is but the easiest way 
to get even quantitatively reliable initial results. 

With running coupling implemented, one has to face the presence of the Landau pole in 
Eq. (|235|l . The main point here is, that, while indeed one does not control what is happening 
at that scale, one is justified in providing a prescription for its treatment. With Qs 3> Aqcd 
the integral on the r.h.s. of the BK equation yields extremely small contributions at scales of the 
order of Aqcd, typically of the order of 10^^ times what is seen around Qs- This is the reason 
why one is allowed to ignore whatever happens there and to freeze any correlators such as N 
at these scales at the values provided by the initial condition. Note that the same reasoning is 
already necessary at fixed coupling in some sense. If one turns on running coupling, the physical 
thing to do is to switch off the influence of the divergence at the Landau pole, for instance by 
freezing the coupling below some scale /io- Of course, one has to make sure that the details of 
where one does that will not affect the results. 

Such an argument about the integrals in the BK equation is valid at zero impact parameter or 
infinite transverse target size, where the nonlinearity is present and sizable everywhere. For large 
impact parameter one necessarily enters the region where the nonlinearity is small and there one 
would remain sensitive to the value of /io- That is the central reason for not considering such 
situations. 

That the effects of running coupling slow down evolution has been known from the earliest 
numerical treatments of BK. The study of the UV phase space in [35] explicitly demonstrates 
the reason for this: Fig. 1191 shows an example for the extrapolation of \/ UsiQsiXY) to its 
continuum value at one value of Y (corresponding to as{Qs{yy') — 0.2) and compares it to the 
fixed coupling case shown earlier. This feature persists over a large range in Qs as shown in [35]. 
The size of active phase space becomes particularly clear in Fig. 1201 which shows activity centered 
within an order of magnitude around Qs- The rightmost panel shows remarkable agreement with 
the analytical results of Mueller and Triantafyllopoulos if their free parameter is used as a fit 
parameter. 

Despite the presence of Aqcd through the running coupling the solutions show a remarkable 
degree of scaling, called near scaling in [35]. Interestingly one may use the evolution rate A to 
check if this behavior is reached. Fig.l21lshows that, starting from G-B+W type initial conditions 
of the form 

Ny,^r^\-e~^''Q'Oro?'\ (237) 
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Figure 19: (From [35]) example of a continuum extrapolation of A/Qs(Qs(y)^) for a fixed initial condi- 
tion at a y-value where aa{Qs{Y)^) = 0.2, compared to the fixed coupling case. 



with different values for the characteristic scale the system quickly converges onto an asymptotic 
line. On this line the system shows scaling behavior to very good accuracy as illustrated in 
Fig. 1221 In fact, the universal features of this equation have been studied far into asymptotic 
region and Fig. 1231 summarizes many of the results. Fig. 1231 extends the range shown in Figl21l 
enormously and adds comparisons with the fixed coupling case and the analytic results of Mueller 
and Triantafyllopoulos. Two simple parametrizations are shown to emphasize how strong the 
change compared to fixed coupling is, not only in size but also in functional shape. Let me 
discuss these issues in turn. 

Fig. 1231 shows A plotted against as{Ql) and the fixed coupling situation will emerge in the 
guise of Qs scale setting. This is represented by the steeply rising straight line (solid) which 
simply corresponds to X{Y) = 2.26as{Q^), where the coefficient was extracted in a fixed coupling 
simulation. Why is this the relevant comparison? Recall the discussion of Sec. 13.41 in the fixed 
coupling case scaling necessarily implies that Qg is of the form 



Qs{r) - 



(238) 



with constant c. The value of c is uniquely determined by the shape of the scaling solution in 
transverse space and turns out to be c = 2.26. Qs scale setting on the other hand amounts to the 
reparametrization (|234|l . and one immediately reads off how A changes from the fixed coupling 
case: 

X{Y) - dy InQsiY') = c a^^'^^ — = c «,(g.(y)2) . (239) 

c remains the same constant as in the fixed coupling case and thus to the fixed coupling curve in 
Fig. 123 Eq. H239() implies that the form of the leading term with c = 2.26 simply corresponds to 
Qs scale setting and thus necessarily suffers the same UV problems as the fixed coupling case. 

The Mueller- Triantafyllopoulos result is shown to agree with the asymptotic line of near 
scaling throughout the whole range (dashed). At large Qs both the fixed coupling and the MT 
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Figure 20: (From [35]) UV dependence of the BK-equation on the asymptotic line after the parent 
dipole scale is taken to determine the coupling. Left: A against a in units of lattice size. Middle: A 
against a in units of Rs, the natural units. Right: the same rescaled with the Y dependence scaled out 
using the results of [73,74], I121II and 112211 with c — 1.5 {Y just leads to a reparametrization). Active 
phase space is limited in the UV (and IR) by 5 to 10 x Qs. 

curves match up. This is due to the fact that there the MT result (which was given in Ea. H121|l 
as a function of Y instead of Qs) at large Y has an expansion in powers of as{Ql) in the form 

\{Y) = cas{Qs^f) + smah corrections . (240) 

This can be understood in the following way: One makes the ansatz 

Arrr 

\{Y) := dy ln(Q,(y)/AQCD) = Aq {as{Qs{Yf)\ " a,{Qs{Yf) - 



with po = (llA^c - 2iV/)/3, and finds 
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qcd; 
(241) 



^^QCD 



^^^^^^ =exp{[(n + l)2A„A"(y-Fo)+[ln(%^)]"-]^} . (242) 



Recovering A from this leads to 



A(r) = ^-^ _n-^iM^^_^ ^243) 

\n + l)2X,{f^Y{Y - Fo) + [H^)r'] """^ 

so that the leading term of Eq. I|121(l corresponds to ri = 1 as advertised, with matching numerical 
coefficients. [I have followed the spirit of [73,74] and absorbed the freedom in the initial condition 
at Yq into the corresponding Y of p21|l .] 

This behavior at large Y can also be expected for the numerical treatment through a simple 
extension of the scaling argument of [55]. Starting from the near scaling observation, one repeats 
the steps leading to Eqns. H101() and \i{)2\ and ends up with an expression for \{Y) of the form 

27rA(r) = ^ / ^^a„(l/r2)(Ar(«^Q?) + N{v'QI) N{r^Ql) N{u^QI)N{v^QI)) +... 

(244) 

^^To construct the solution introduce f(Y) := ln(Qs(y)^/AQQjj) and solve by separation of variables in / and 
Y. 
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Figure 21: (From [35]) A = 9i-lnQs(V) for dif- Figure 22: (From [35]) Approximate scaling be- 
ferent initial conditions. All curves starting from havior of AV.^y on the asymptotic line of Fig. 1211 
Gaussian initial conditions stay below an asymp- 
totic line along which we find (approximately) scal- 
ing dipole functions N. 

in which the dots stand for the small scaling violations. The only source for the difference between 
the fixed coupling expression Eq. Hl()2(l and 1)244(1 or the corresponding lines in Fig. lies in 
the r dependence of the running coupling that is now inside the integral. To obtain the strong 
relative reduction in active phase space observed in Fig. 1191 as must vary considerably over the 
range of scales contributing to the integral in the fixed coupling case. Wherever this happens, 
X{Y) will be a nontrivial fimction of as(Qs(^)^)- Conversely, at asymptotically large Qs, where 
the as(l/r^) probed in Eq. (|244|1 turn essentially constant over a range comparable to the phase 
space needed in the fixed coupling case, one should return to a situation as in Eq. H24()() . 

It is essential to remind oneself here that the limit in which the leading term of (|240|1 is 
sufficient is not under control. The price to pay for this behavior at large Y is the reopening of 
phase space to the size encountered already at fixed coupling. 

As in [35], this is illustrated by tabulating (see table ^ a measurement of active phase space 
in terms of the values of a/R needed to find A within a given percentage of its continuum 
value at different as{Qs{Y)^). This is then compared with the fixed couphng result. In the far 



deviation from 


as(O') = .2 




= .15 


Qs(Q? 


= .12 


as(Q? 


= .1 


fixed coupling 


continuum value 


A/a = .845 


A/a = 


.96 


A/a = 


1.052 


A/a = 


1.14 




10% 


a/R = 0.17 


a/R = 


0.15 


a/R = 


0.12 


a/R = 


0.1 


a/R = 0.0011 


5 % 


a/R = 0.1 


a/R = 


0.08 


a/R = 


0.066 


a/R = 


0.05 


a/R = 0.00013 


1 % 


a/R = 0.027 


a/R = 


0.015 


a/R = 


0.007 


a/R = 


0.0012 


a/R = 8.2 • 10" 



Table 1: Phase space needed to extract continuum values for A is reduced for larger values of the 
couphng. The fixed coupling situation reemerges asymptotically for large Qs, i.e. small coupling. 
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asymptotic region one is back to the unsatisfactory situation that important corrections are most 
hkely missed, with the additional drawback, that now one has no real idea of what those might 
be. Fortunately, this would appear to occur only at such extremely large values of Qs that this 
concern is purely academic. 

2, / GeV 

10** lo' 100 10 5 




Figure 23: (From [35]) XiY) as a function oi aaiQsiYY') with running coupling. A perfect overall fit is 
obtained with the MT result Eq. 112211 . The shape of the curve only depends on c, shown is a fit with 
c = 1.5. Simple parametrizations emphasize how drastically the functional shape is changed away from 
a simple c as{Ql). 



The fact that the asymptotic expansion of Mueller and Triantafyllopoulos of \{Y) -Eq. H121(l - 
apparently contains only a single term beyond the questionable Qa scale setting result, and yet 
provides a fit over the whole Qg range with its large deviations from the leading term, may cause 
confusion. The resolution to this simply is hidden in the change of variables Y , used by MT, 
and Qs, as used here. Any attempt to translate Eq. H121(l from Y to Qs or as{Ql) leads to a 
infinite series that will converge only slowly, so that it is in fact wrong to think of this as a single 
correction term from this perspective. 

The discussion thus far has been about the behavior of the dipole within one of two orders 
of magnitude of Qs- Careful studies have also been performed at much shorter distances [61,95] 
in the scaling window, where gluon densities are still small, but scaling and the presence of Qs 
still plays an important role as it leads to behavior of the form H123|l with 



= c 



-In 



1 



RsiYy 



1 



1^0 



Both studies attempt to extract the anomalous dimensions of Eq. H123|l , starting from distances 
typically a few orders of magnitude smaller than Rs- 

They find that for this quantity which only makes sense at such short distances, running 
coupling effects remain important over an enormous range. Residual information about the 
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initial conditions used appears to be more persistent than close to Qs ■ At the same time running 
coupling effects remain particularly important in this range. Both these facts can be understood 
by looking at the BK equation at running coupling and the integral expression for X{Y). Writing, 
for instance. Ea. H236|l in terms of rescaled variables f :— rQs{Y) etc. one obtains 

dyNr = ^ [d^z ICi^iy |a^^ + Ni,-Nr- iV^ivJ ' (^45) 

) — In(r^) J 

^^QCD 

Clearly, with increasing Qs the importance of running coupling effects (in the sense that goes 
beyond Qs scale setting) will be pushed further and further away from the scaling region where 
r = rQs{Y) > 10~^, while at shorter distances, where the anomalous dimensions are measured 
form a scaling fit to H123() . f dependence of the coupling remains important. These features 
are reflected by the numerical results of [61,95] which emphatically confirm the existence of the 
scaling window: they cannot even find any evidence that scaling of the form (|123|l would break 
down at the upper end of the scaling region.''^ 

To summarize the results of this section, one should state that with parent dipolc running, 
the contributions to A arise from within little more than an order of magnitude of the saturation 
scale, as Fig. 1201 shows impressively. [This concerns the region in cutoffs over which there is 
an appreciable change in the value extracted.] This indicates that the dominant phase space 
corrections are indeed taken into account and the basic physics ideas that had initially motivated 
the density resummations are realized. While the IR side of the resummation was working 
already at fixed coupling, the UV side of phase space is only tamed if the coupling runs at a scale 
directly visible in the diagrams comprising the r.h.s. of the evolution equation such as 
the parent dipole size. This is equivalent to incorporating subleading terms in the asymptotic 
expansions H240|) at least to the order done by Mueller and Triantafyllopoulos in H121|) . Only 
with both the nonlinearities and a running coupling that goes beyond Qg scale setting in the 
sense explained above are we justified to claim that the evolution is dominated by distances of 
the order of Rg and thus the coupling involved in typical radiation events during evolution is of 
order as{Qs{Y)'^). [Emission from very small dipoles would simply not be sufficiently suppressed 
with Qs scale setting.] While the result of [73,74] provides a physically motivated explanation of 
the curve, a convenient way to convert this into an explicit function of Qs{Y) does not exist. For 
this purpose the fits shown in Fis'. l23l Drovide a useful, quantitative, rule-of-thumb representation 
for the numerical result. The existence of the scaling window is numerically confirmed: Scaling 
features extend into the low density region and persist throughout the scaling window to distances 
much shorter than the saturation scale. 



7.2 Evolution induces longitudinal correlations: breaking of A s -scaling in DIS 

Already in the introduction it has been argued and illustrated in Fig. |H1 that for uncorrelated 
scattering centers Qs should scale with larger atomic number like A^^^ as shown in Eq. 0. This 
has recurred in the context of the MV model and its generalizations rSec l2.2(l : generically one 
expects that going to larger nuclear targets ought to have an enhancing effect on the importance 
of the nonlinearities encountered. To recall the argument: If one is at small enough Xbj (high 
enough energies) for the projectile to punch straight through the target, going to a larger target 



means interaction with more scattering centers as indicated in Fig. 5(b) Taking these to be 



color charges located in individual nucleons inside the target nucleus the interaction will occur 



■^^At the short distances involved there, this might be a problem with short distance corrections entirely absent 
in the BK equation, but this is a question that goes beyond the existence and relevance of the scaling window 
per se. 



80 



with more and more uncorrelated color charges the larger the target. Viewed from up front this 
amounts to a smaller correlation length in the transverse direction. In terms of Qg this would 
lead to a variant of relation © , 

Qf{Yf ^cAiQPiYf (246) 

where c is some, yet unspecified constant that takes into account nontrivial geometry and other is- 
sues not incorporated in the McLerran-Venugopalan model (see also the discussion below Eq. I)28|l 
and [42]). Note that none of these arguments would give a clear delineation of the range of Y 
values for which this would be appropriate. It is clear, however, that in the context of small 
Xhj evolution such a statement becomes even less trivial: the whole idea of deriving a small Xhj 
evolution equation is based on the premise that the logarithmic corrections determining the evo- 
lution step arc entirely target independent. The only place where target dependence is allowed 
to enter is the initial condition as already indicated in Eq. (jnj. A careful determination of c can 
not be done without phcnomcnological effort, but a discussion of how evolution changes such a 
relationship can be done on the basis of the evolution equations alone. This is what I focus on 
here. 

In simple cases the two notions -evolution and scaling at all Y- are not conflicting with each 
other: indeed a rescaling like in Eq. (|246|l is fully compatible with fixed coupling evolution on 
the asymptotic line, where (in the continuum limit) 

0,(r) ^ e^(^-^°)g,,(ro) (247) 

with a universal value for A, so that a simple rescaling of the properties of the initial condition 
according to 

Q^{Yo)' = cAiQP{Yor (248) 

automatically yields Eq. H246() at all Y declaring it fully compatible with evolution. 

At running coupling, however, this ceases to be the case. The straightforward physical argu- 
ment is that evolution is characterized by q;s(Qs)j which initially is larger for protons than for 
nuclei, because of the difference in saturation scales at the initial condition. In this situation 
evolution will be faster for smaller targets and eventually catch up with that of the larger ones. 

Technically, with running coupling (and again on the asymptotic line), one finds a Y depen- 
dent A that may be parametrized through some function / of as{QsiY)'^) as in the previous 
section: 

A(y) = dy InQsiY) = fiasiQsiYf)) • (249) 

In this generic case 

0,(y) = e.^^o^^'/(-(Q=(^)^)Q,(y„) (250) 

and a rescaling of (5^o)^ will not simply factor out as in Eq. H247|l . Therefore a simple rescaling 
like Eq. (|246(l will conflict with evolution. All phcnomcnological fit functions encountered in 
FigESlshare this property, which is at the heart of the reasoning used by Mueller in [96] to argue 
that evolution will eventually erase A-dependence in Qs- 

Instead of a numerical analysis (c.f. [95]) and purely to illustrate the point, one may trace 
the memory loss with a simple parametrization like Eq. H241|l . equivalent to a y dependence of 
Qs according to Eq. H242(l . Recall that for n = 1/2 this should be even quantitatively correct 
over the physically interesting domain. 

Scahng in the A-dcpendcncc at Yo according to Qf{Yo)^ = cA^QP^Yoy , one uses Eq. H242|l 
to predict the combined Y and A dependence of the saturation scale. The result is in qualitative 
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agreement with [96] : the ratio of saturation scales in the bulk of the region shown in Fig. 
behaves like 



Qf(r)2 



(n + l)2Ao(f)"(r-ro) 



[ln( 



^^QCD 



exp 



(n + l)2Ao(fr(r-yo)+[ln(2^)]"^^ 



(251) 



This is plotted in Fig. |^ over a very large range in Y. Clearly larger nuclei suffer stronger 
initial "erasing." It should be kept in mind that realistically one should not expect real world 
experiments to cover more than a few orders of magnitude in Xhj and thus a small interval in 
Y-Yo. 
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Figure 24: (From [35]) A-dependence for proton saturation scales from 1 to lOOGeV. c parametrizes 
possible shape dependence as explained in the text. 



8 A brief foray into phenomenology 

Phenomenology in the context of the color glass condensate is a very hot topic, in particular in 
connection with present activities at RHIC. As a consequence, many of the topics arc not yet 
settled enough to be included in a review. Nevertheless, this review would be incomplete without 
trying to give a flavor of the ideas and experiments being discussed. 

The key point in all of them is the appearance of Qs and associated scaling in the interpre- 
tation of experimental data. This is done with respect to the Xbj dependence as well as the A 
dependence where available, both in the saturation region and the scaling window. 

The prime example still is the application of saturation ideas to HERA data at Xbj < 10^^. 
This will be the first topic discussed below in Sec. 18.11 

Nucleus nucleus collisions form a much more complicated environment than Icpton nucleus 
collisions and the step from fully developed theory to comparison with data is much larger. 

There is one approach to heavy ion data which is directly build on the technology described 
in this review, and that aims at the initial condition of heavy ion collisions. This will be briefly 
touched upon in Sec. 18.21 
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A different approach based on (fc =) fcT-factorizcd formulae for particle production in terms 
of unintegrated gluon densities as defined in Eq. (|42|l [61,85,97-102]. allows to apply the results 
of the evolution equations to AA' experiments: 

i?^ = ^^^^x l^'dfc^a, (252) 

Despite some obvious shortcomings concerning perturbative and hadronization corrections, a 
treatment based on this formula is currently the best one can do. Keeping this in mind, the 
formula imprints the consequences of any Qs scaling found in unintegrated distributions onto 
transverse momentum spectra in nucleus nucleus collisions. 

A very illustrative example in this context is the Cronin effect and its expected behavior under 
evolution discussed in Sec. 18.31 Particle multiplicities are another example, but indications are 
that this topic will come into its full right only with LHC energies. These and other ideas relevant 
to build a solid understanding on how and exactly in which kinematical domain the CGC shows 
up in experiments will briefly reappear in Sec. 18.41 



8.1 Geometric scaling in HERA and lepton nucleus collisions 

The discovery of geometric scaling in HERA data at Xbj < 10~^ as prompted by the Golec- 
Biernat and Wiisthoff fit [32,34] was the first phenomenological application of saturation ideas in 
small Xbj physics, predating all the theoretical developments around JIMWLK and BK evolution 
equations that demonstrated the scaling features of these evolution equations. It brought with 
it the by far most efficient parametrization of HERA data at Xbj < 10~^. Recall the DIS cross 
section formula from the introduction (with the momentum fraction integral restored): 

1 

atot{Y,Q^) = j d?r j da\ii^^,{ay,Q^)\'' aai^oXeiXy) (253) 



where a is the longitudinal momentum fraction of the q or q, and z their relative transverse 
separation. and Xbj-dependence is clearly separated into wave function and dipole cross 
section respectively. The G-B+W fit parametrizes the latter via an Xbj-dependent saturation 
scale Qs according to H46|l and a normalization factor to carry the dimensions from the impact 
parameter integral: 

/-SA 2\ f 72 7 Mi^^^xU^). G-B+W -r^O^(y)/4^ /c^r a\ 

cTdipoioiY, r-") ^ J d-'b ( — > cro(l - e '^'^'^^^^) . (254) 

At the same time they preempted the fixed coupling scaling form for Qs of Eq. (|103|l . Qs{Y) = 
P^HY -Yo) Q ^^Yq) = {xhiQ/x\jj)^Qs{xo). This fit ansatz is just a particular realization of the scaling 
behavior later found to result from JIMWLK and BK, both at fixed and running coupling. Suc- 
cessful saturation based fits incorporating scaling forms like l|123(l have also been performed in the 
meantime [103]. The most important phenomenological result, however, is independent of the ex- 
act shape of the dipole cross section, it is the simple fact that HERA data show geometric scaling. 
Once the dipole cross section depends on Xbj only via Qs{Y) as crdipo\c{Y, r^) = (Tdipoicir^ Ql{Y)), 
geometric scaling emerges directly. Rescaling the r integral in Ea. H253|l and using the fact that 
the wave function factor scales like f{r'^Q^)/r^, this implies 

'TtotiY,Q^) ^ a,,,(Yo,-^^Ql{Yo)) (255) 
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x-evolution,X = 0.188, j: = 5.5x10" -0.025 BK with running a, ^ = 0.182 




Q" / GeV" Y 



Figure 25: Geometric scaling in HERA data. Shown is the structure function F2, which at small a-bj is 
simply related to the cross section by a rescaling by Q^: a{Y,Q^) ~ F2{Y,Q^) ■ Q^. 

and should allow to map all the cross section data onto a single scaling curve. This is shown in 
Fig. 1251 While the explicit parametrizations shown in Fig. |211 make use of BK evolution results 
from [35], this is not essential, at least for the scaling plot on the right. In fact, successful fits of 
HERA data based on CGC evolution exist [67,103] and they are of a quality comparable to that of 
DGLAP fits. It is not the quality of the fit to the data (leaving aisde comments on negative gluon 
distributions often seen in global DGLAP analyses) that distinguishes the CGC perspective, it 
is the natural emergence of Qs and the scaling features which set it apart. Scaling features have 
since been looked for also in older NMC [104] and E665 [105, 106] data for the ratios of nuclear 
to deuterium structure functions 2F^ /AF.P . Despite the fact that both Xbj and ranges are 
such that CGC saturation arguments remain marginal despite the moderate enhancement via 
A, both [107] and [108] find that scaling in both x^j and A holds surprisingly well. The scaling 
in A was fitted with a fixed power close to 1/3 ignoring the refinements of Sec. 17.21 Due to the 
limited range in Xhj this can not be taken as an experimental confirmation of CGC physics - 
such would require a measurement of Xhj dependence and this is clearly a reason to explore this 
question more closely in experiments at higher energies. This is being actively pursued. One of 
the options is the Electron Ion Collider (EIC) in planning at BNL [109] as a machine dedicated 
to such experiments or via pA collisions at the LHC where higher energies are available at the 
price of a more complicated projectile. 

8.2 Gluon production in the initial stages of heavy ion collisions 

The evolution equations have mainly been discussed in the context of eA collisions, but it should 
be intuitively clear that a new perturbatively large scale Qs, present in the nuclear wave function 
at small Xbj, would also manifest itself in AA' collisions. The original idea is laid out in [23], 
where the incoming nuclei are treated as sources of strong color fields in the spirit of the MV 
model. While the nonlinear nature of the fields does not allow to simply add them, it is still 
possible to add the color currents that create them and then to solve the resulting Yang Mills 
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equations. This will only modify the solutions after they are in causal contact, i.e. in the 
forward lightcone of the collision point. At the edge of this forward cone the fields generated by 
the individual nuclei serve as an initial condition for what happens after the collision for a given 
source configuration. The source configurations themselves are subject to ensemble averaging, 
say with (|29|l (per nucleus), if one stays within the MV model or with the corresponding field 
averages expressed alternatively via 6, a or [/, if one goes beyond this model to incorporate Xbj 
dependence via BK or JIMWLK equations. 

[23] provides the generic framework but only attempts perturbative solutions valid for large 
times, however, the field equations for large densities are truly nonlinear in the very initial stages 
of the collision of interest. In this region numerical methods are required. These have been 
developed and used by [110,111]. More recent work discusses the possibility to create "elliptic 
flow" [112,113] and independent checks have been provided in [114]. The main interest of such 
studies is the creation of high density effects in the very initial stages of heavy ion collisions, 
which would appear to be necessary to explain the very short equilibration times often quoted 
as necessary to explain the collective phenomena supported by data. 

8.3 Evolution erases the Cronin effect 

The last experimental topic to be discussed here that has met with strong interest is the question 
for the energy dependence of what is known as the Cronin effect. I should begin this section 
with the warning that all input and results in terms of the CGC are formulated on a partonic 
level, while the underlying definitions of cross sections (|256|l and actual measurements are done 
on a hadronic level. There are strong indications that there may by important contributions 
from hadronization to the measured quantities. For this reason one might want to talk about 
what is discussed below as a Cronin-like effect. Nevertheless, the main message that small Xbj 
evolution has a specific and rather strong effect on transverse momentum spectra remains an 
exciting result. 

The Cronin effect refers to an enhancement of the number of particles produced in the in- 
termediate transverse momentum range of nucleus nucleus collisions over those in proton proton 
collisions which is commonly attributed to a momentum broadening due to multiple rescatterings 
in the final state. At first sight this would appear to be at odds with the ideas of saturation 
which would lead to the expectation that particle production should be suppressed for momenta 
below Qs and even, as argued by [115] within the entire scaling window. This, however, is but 
less than half the picture. 

The first observation is that the original MV model does contain Cronin type enhancement 
[49, 116, 117]. This is based on a sum rule [43] for the gluon distribution in the MV model that 
equates the total number of gluons in a (large) nucleus with A times that of a nuclcon. At the 
same time the momentum spectrum is distorted by eikonalization: small momenta are indeed 
depleted, while at very large momenta a perturbative tail remains undisturbed. This can 
be readily recovered from (|44|l together with H45|l under the assumption of A^/^ scaling of the 
exponent in either gluon distribution or dipole cross section. (The remaining A^/^ comes from 
the (Pb integral as a global prefactor.) By necessity the "excess" has to accumulate at moderate 
transverse momenta which leads to Cronin enhancement. 

The second observation comes from the fact that the sum rule is broken by evolution. One 
aspect of this is simply the modification of A-dependence of the saturation scale through evolution 
reviewed in Sec. 17.21 The effect on the gluon distribution and thus the Cronin effect is even much 
more pronounced: The Cronin peak is suppressed over a very short range in rapidity. 

^*If applied to the gluon field A, this is a gauge dependent statement! 
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The underlying object measured is a ratio RAB.CD{y,P,b) that compares numbers of pro- 
duced particles per rapidity Y (i.e. Y for present purposes), transverse momentum p and cen- 
trality b in the collision of two nuclei A, B with the same in the collision of another pair C, D as 
a function of transverse momentum, both normalized to the number of collisions Nco\\{b)'. 



Ra B;C d(Y,p, b) 



(^con(fc))c+u dNA+B/dYd^pd^b 
{N,on{b))A+B dNc+o/dYd^pd^b 



(256) 



For practical measurements (|256(l is too differential and so simplified variants are used. On the 
one hand; there is the ratio Ra b hi which the multiplicities of (|256() are integrated over impact 
parameter 6, and C, D are both taken to be protons for reference. This is the standard way to 
show the Cronin effect. On the other hand, there is Rep {c,p stand for central and peripheral 
respectively) in which (A, B) = (C, D) and the multiplicities in numerator and denominator are 
integrated over central and peripheral regions respectively. [The same applies to the average 
number of collisions used to normalize.] Both quantities should show similar trends, loosely 
speaking, the peripheral region looks like a collision with a smaller nucleus. 

i?d Au had been studied carefully in the context of the BK equations in [69] . For this the- 
oretical study initial conditions of the form H44|l together with (|45() were used to provide an 
initial condition that shows Cronin enhancement. The result is shown in Fig. 1261 and shows the 
fast disappearing Cronin enhancement. Again evolution is considerably slower in the running 
coupling case. 
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Figure 26: (Courtesy of the authors of [69]) Suppression of the Cronin peak through BK evolution both 
at fixed and running couphng. Running coupling evolution is slower than the fixed coupling case as seen 
already for generic quantities in Sec.|3 



Both i?d Au and i?cp have been measured by the BRAHMS experiment at RHIC [118] and 
are shown in Fig. 1271 top and bottom strips respectively. In both of these measurements it is 
the depletion at forward rapidities that is the interesting effect: Other approaches [119,120] that 
concentrate on multiple scattering to reshuffle particles in the momentum distributions typically 
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Figure 27: Brahms data [118], Top strip: i?d Au as a function of p — pr for different pseudorapidities rj 
(here to be identified with Y). Bottom strip: The same for Rep in d + Au. 



do not take into account the nonlinearity of the small Xhj evolution equation which leads to 
the effect shown in Fig. 1211 Instead the Cronin enhancement tends to increase, reflecting the 
increasing number of scattering centers at small Xbj ■ 

Experimental work and theoretical discussion on this topic is still ongoing, see for exam- 
ple [121] for a recent fit of available data. 

8.4 Many ideas abound 

The list of ideas given above is by no means exhaustive. The effect of Qs on particle multiplic- 
ities in pA and AA' collisions has been discussed as well, initially for RHIC, where they have 
been essentially ruled out for the central rapidity region. This is in keeping with the tentative 
interpretation of BRAHMS data discussed in Sec. 18.31 but judging from the same data, one 
would expect this to be highly relevant at the LHC at central rapidities for which predictions 
exist [122]. The applications typically combine methods to calculate particle yields with classical 
Glauber formulae for multiple scattering situations -used in order to understand the reaction 
geometry- with a factorized ansatz for inclusive production in terms of the unintegrated gluon 
densities of the participants as discussed in Eq. 1)252(1 . The latter are then parametrized using 
information about the CGC, typically with some parametrizations for the energy dependence of 
the saturation scale. This should give first ideas about the effects to be expected, but clearly 
could profit from a closer integration with the evolution equations and the concepts of Sec. 18.21 

Recently two particle and jet azimuthal corrections have been suggested [123, 124] to further 
study the onset of parton saturation suggested by the BRAHMS data, ideas to measure the 
transverse size of the area populated by gluons have been formulated [125] and too many others 
to give a complete list of. 

It is clear, however, from the added number of theoretical complications in heavy ion exper- 
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iments that a lepton nucleus experiment at high energies would be an important tool to get a 
clean baseline for all the applications to these more complicated experiments. This is the goal of 
the EIC project currently in the planning stage at BNL. 

9 Conclusions 

The concepts to describe saturation effects in small Xhj have come a long ways since the early 
discussions sparked by the Gribov-Levin-Ryskin (GLR) equation [126]. Stages on the way are 
represented by the Mueller-Qiu [16] equation (aiming at the same double logarithmic limit as 
the GLR equation), Mueller's original dipole model [17, 18, 127, 128] designed to give evolution 
towards small Xhj and the McLerran-Venugopalan model meant to describe density effects at 
small but fixed Xbj in large nuclei. 

JIMWLK and BK evolution equations have provided a new framework, derived from small 
coupling QCD, that allows to take into account gluon induced density effects at small Xbj- 

The theoretical treatment contains two types of nonlinearities, eikonalization of "soft" gluon 
fields into path ordered exponentials or Wilson lines as a first stage and a second stage that even 
exponentiates these objects -inside of shower operators- to describe emission of "hard" gluons. 

The first stage is sufficient to obtain bounded results for cross sections as evidenced by the 
dipole formula for the DIS example Q or Sudakov form factors in jet observables. Such expres- 
sions, since bounded, are always efficiently parametrized as exponentials as in Eqns. (j^SJ, without 
any reference to a specific averaging procedure for the gluon fields - interpretation in terms of a 
saturation scale as in Eq. H44(l is most natural in the small Xbj context. 

The second stage exponentiation via the shower operators [54] implements x^j dependence 
of objects like the photon wavefunction in DIS. Via translation through Langevin formulations 
[13,75] and JIMWLK equation [7-15], this generates an infinite hierarchy of coupled equations for 
cross sections, the Balitsky hierarchy, which reduces to the BK equation, if correlators factorize 
into two point functions. 

In whichever guise this is presented, the Xhj dependence is characterized by the emergence of 
scaling: universal behavior that ties the Xbj dependence to dependence on the saturation scale. 

The evolution equations impose this behavior quite quickly: Near the saturation scale details 
about initial conditions are efficiently erased, as is evidenced by the quick convergence onto the 
near scaling curve shown by the plot of the evolution rate X{Y) = dy ^nQs{Y) in Fig. |231 Once 
near scaling is reached, the scaling form only changes very slowly - scale breaking is present but 
weak [107]. This feature also manifests itself in other quantities such as, for example, anomalous 
dimensions as measured from BK evolution by [61,95]. A subtlety that was manifest -although 
unappreciated- already in numerical simulations at fixed coupling is the importance of running 
coupling to tame UV phase space. Only with running coupling implemented is the evolution 
activity centered around the physical scale Qg. 

The generic structures of the phase space diagram Fig.0]with saturation, scaling and extended 
scaling region have all been seen numerically. The perhaps simplest physical interpretation is 
this: Starting from a dipole function that is saturated at large dipole sizes and perturbatively 
small at small dipole sizes, evolution is driven by contributions from small sizes. This follows the 
BFKL pattern of exponential growth. Creation of large objects that would fall into the already 
saturated domain is prohibited by the nonlinearities, corresponding contributions are absorbed 
by the medium. Due to the growth at small distances, densities rise also there and the boundary 
where the nonlinearities start to act readjust. This happens selfconsistently in JIMWLK and 
BK. In [73, 74] the picture is taken literally and an even quantitatively comparable result is 
achieved by combining BFKL evolution with an absorptive boundary which is selfconsistently 
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adjusted to limit growth. This type of analysis leads to scaling forms for the dipolc cross section 
valid at small distances, and predicts the size of the scaling window via the corresponding BFKL 
diffusion radius. 

Scaling features as found in these evolution equations had been seen earlier in HERA data at 
small Xhj and called geometric scaling. This has encouraged phenomenological applications to 
various observables like particle multiplicities or the Cronin effect. Besides its Xhj dependence, the 
A dependence of Qs as implied by the MV model and somewhat modified by evolution [35,95,96] 
is an interesting tool to tie together experimental results in Icpton proton, lepton nucleus, and 
nucleus nucleus experiments. While qualitative features in existing experiments -HERA and 
RHIC featuring most prominently on this list- appear to be compatible with state of the art 
CGC phenomenology, matters arc far from settled. Experiments at smaller Xbj, like the LHC, 
and simpler environments, like the EIC with its electron projectiles, are needed to fully clarify 
the situation. Current experimental activities at RHIC and phenomenological work to match 
and assess the enormous amount of new results, for example by suggesting new measurements 
and phenomenological consequences, are very encouraging. 

At the same time it would help to close the gap between reliable theory and phenomenological 
applications. Our understanding of small Xbj evolution equations at this point is limited to the 
central part of Icpton nucleus collisions. The technological tools and insights gained through 
comparison with soft emission phenomena in non-global jet observables should allow to formulate 
new, IR safe observables together with the associated evolution equations also in the small Xbj 
case. This field has matured from its origins, but new developments are just around the corner. 
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